Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,000 C→G 5.1% P219P (CCG→CCC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0000C→G 5.1% 70.1 ‑2.5 ‑0.1798P219P (CCG→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (27/15):  ref base (371/383):  total (400/398)
Fisher's exact test strand distribution p-value = 8.01e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.40e-01

Only 100 of 1658 total aligned reads displayed.

TCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCG  >  NC_014034/221900‑222100
                                                                                                    |                                                                                                     
tgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                                       >  1:803581‑M1/57‑101 (MQ=255)
      cggtcgggtccgggggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                                                                                 >  1:3353720‑M1/51‑101 (MQ=255)
         gcggatccgtgggatcggtggggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATGGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                                                              >  1:4019863‑M1/48‑101 (MQ=255)
             gtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                                                                          <  2:1177301‑M1/58‑1 (MQ=255)
              tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                                         <  2:3778276/101‑1 (MQ=255)
               ccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                                                        <  2:223599‑M1/60‑1 (MQ=255)
                  tgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                     <  1:2322739‑M1/63‑1 (MQ=255)
                   gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                    <  1:4417150‑M1/64‑1 (MQ=255)
                        cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                               >  1:4146503‑M1/33‑101 (MQ=255)
                         ggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                              <  2:3700935‑M1/70‑1 (MQ=255)
                          aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                                                                             <  2:2936354/101‑1 (MQ=255)
                           tcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTccggg                                                                            <  2:1239703‑M1/72‑3 (MQ=255)
                             gggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                                                                          <  1:42604‑M1/74‑1 (MQ=255)
                              ggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGa                                                                         <  2:2378565‑M1/75‑1 (MQ=255)
                               atccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                        >  1:1734704‑M1/26‑101 (MQ=255)
                                tccgttggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                       >  1:1487105‑M1/25‑101 (MQ=255)
                                tccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                       >  2:2863924‑M1/25‑101 (MQ=255)
                                tccgtcggagcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                       >  2:3123421‑M1/25‑101 (MQ=255)
                                    tgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                   >  1:652221‑M1/21‑101 (MQ=255)
                                    tgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                   <  1:409390‑M1/81‑1 (MQ=255)
                                          cggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTC.GGATcc                                                            <  1:3546489‑M1/87‑1 (MQ=255)
                                           ggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCg                                                           <  1:3211912‑M1/88‑1 (MQ=255)
                                           gcgcgggcccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCg                                                           <  2:4091726/91‑1 (MQ=11)
                                           cGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGCTCGGTCGGGTCCGTGGGATCGGTC.GGGTCCg                                                           <  2:1862255/101‑1 (MQ=11)
                                            gtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTC.GGATCCGt                                                          <  1:3459143‑M1/89‑1 (MQ=255)
                                             tcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTC.GGGTccggg                                                         <  2:1942607‑M1/90‑3 (MQ=255)
                                             tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTg                                                         >  1:1807487/1‑101 (MQ=255)
                                             tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTg                                                         >  2:4035154/1‑101 (MQ=11)
                                              cgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTC.GGATCCGTgg                                                        <  1:2579821‑M1/91‑1 (MQ=255)
                                               gtgtccttgggatcggtcgggtccgtggggtccgtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTggg                                                       <  2:86152/92‑1 (MQ=11)
                                               ggatccgtgggatcggtcgggtccgtgggatcggtggggtcggtgGGAACGGCGGGGTGCGTGGGATCGGTCGGGGCCGTGGGATCGGTG.CGATCTGTggg                                                       >  1:2103484/10‑101 (MQ=0)
                                                gatccgtgggatgggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGa                                                      >  1:2417908/9‑101 (MQ=11)
                                                 gtccgtgGGGTCGGTCGGGTCCGTGGTATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGAt                                                     >  1:2509608‑M1/8‑101 (MQ=255)
                                                 gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTC.GGATCCGTGGGAt                                                     >  2:1657793‑M1/8‑101 (MQ=255)
                                                  tccttggtatcggtcgggtccgtgggatcggtcgggtccgtgggGTCCGTCGGAACCGTGGGATCGGGGGGGCCCGTGGGGTCGGTT.GGATCCGGGGGact                                                    >  2:927021/9‑99 (MQ=14)
                                                  tccttggcgtcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGAtc                                                    <  2:4270357/92‑1 (MQ=11)
                                                  tccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTC.GGATCCGTGGGAtc                                                    >  1:881033‑M1/7‑101 (MQ=255)
                                                  tccggggggtcggtcgggtccgggggatcggtcgcgtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGAtc                                                    <  1:3827873/92‑1 (MQ=11)
                                                   gcgggggatcggtcgggtccgtgggatcggtcgggtccgtcGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGAtcg                                                   <  2:2635600/96‑1 (MQ=11)
                                                   ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTGGGAtcg                                                   >  1:3742973‑M1/6‑101 (MQ=255)
                                                    cttgCGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGAtcgg                                                  >  2:3260368‑M1/5‑101 (MQ=255)
                                                    cgtggtatcggtcgggtccgtgggatcggtcgagtccgtgggGGCCGTCGGATGCGTTGGATCGGTCGGAGCCCTGGGATCAGTC.CGATCCGCGGGAtcgt                                                  >  1:307578/7‑100 (MQ=2)
                                                     gtgggatcgggcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTGGGAtcggt                                                 <  1:3959890/98‑1 (MQ=11)
                                                     gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGAtcggt                                                 >  1:1661257/1‑101 (MQ=255)
                                                      tgggatcggtcgggtccgtcggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTC.GGGGCCGTGGGGtcggtc                                                >  2:3871707/3‑101 (MQ=11)
                                                      tgGGATCGGTCGGGTCCGTGGGAACGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGAtcggtc                                                >  1:1912451‑M1/3‑101 (MQ=255)
                                                       gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTGGGAtcggtcg                                               <  2:2462359‑M1/100‑1 (MQ=255)
                                                        ggggcggtgggatccgtggggtcggtcggggccgtgggatCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGAtcggtcgg                                              <  2:1784105/97‑1 (MQ=11)
                                                        ggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTGGGATCGGTC.GGAGCCGTGGGAtcggtggg                                              >  2:3305057/1‑101 (MQ=11)
                                                         ggtcggccgggtccgttggagtggggggatccgtgggaTCGGTTGGGTCCGTAGTATCGGTCGGGTCCGTCGGATCGGTC.GGATCCGTGGGATCGGTCggg                                             <  2:54551/99‑1 (MQ=2)
                                                         gatcggtcgggtccgtgggatcggtcgggtccgtggGATCGGTCGGATCCGTGCGATCGGTCGGGTTCGTGGGTTCCGTC.TGATCCGTggggcccgtgggg                                             >  2:621675/1‑91 (MQ=11)
                                                            cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGATcc                                          <  2:454772/101‑1 (MQ=11)
                                                            cggtccggtccgtgggatcgttcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTC.CGATCCGTGGGATCGGTCGGATcc                                          <  2:1441491/101‑1 (MQ=14)
                                                             ggtcggatccgtgggatcggtcggatccgtaggatcGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTC.GGATCCGGGGGATCGGTCGGGTCCg                                         >  2:2300378/1‑101 (MQ=11)
                                                              gtcgggtccgtggggtccgtcggatccgtgggatcgGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTTCGTGGGGTCGGTCGGGTCCGt                                        >  2:151624/1‑101 (MQ=11)
                                                              gtcgggtccgtgggatcggtggggtccgtgggatcgGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGGTCCGTCGGGATCCG.GGGGTTGGGCGGGTcccg                                        >  1:477624/1‑99 (MQ=11)
                                                              gtcgggtccgtgggatcggtcgggtccgtgggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGC.GGGTCCGTGGGGTACGTCGGGTCCGt                                        >  2:2107222/1‑101 (MQ=11)
                                                               tcgtgtccttgggatcggtcgggtccgtggggtccgTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGtg                                       <  2:135086/101‑2 (MQ=14)
                                                                cgggtgggttggatcggtcggggccgtggggtccgtCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGtgg                                      <  2:2772329/101‑3 (MQ=11)
                                                                cgggtccgtgggatcgttcgggtccgtgggatcggcCGGGTCCGGGGGATCGGTCGGGTCCGTGGGGTCCGTC.GGATCCGTGGGATCGGTCGGATCCGtgg                                      >  1:3898669/1‑99 (MQ=14)
                                                                 ggctccgttggatcggtcggagccctgtgatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.TGATCCGTGGGATCGGTCGGGTCCGTGgg                                     <  2:693373/98‑1 (MQ=14)
                                                                 ggatccgtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTC.GGATCCGTGGGAGCGGTCGGATCCgaggg                                     >  2:3464358/4‑97 (MQ=14)
                                                                  gatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTC.GGCTCCGTCGGGTCGGGTGGCTCCggggga                                    >  1:2691443/3‑96 (MQ=11)
                                                                   gtccgttgggtccgtctggtccgtgggctcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTGgggt                                   <  2:1319800/101‑3 (MQ=11)
                                                                   atcccggggatcggtcggatccggaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTCGGAt                                   <  2:1022703/95‑1 (MQ=11)
                                                                     tcgtggggtccgtcgggtccgtggggtcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTCGGATCGGTGGGATCCGTGGGAtcg                                 <  2:4356456/100‑1 (MQ=11)
                                                                     ccgtgggatcggtcgggtgcgtgggatcgggcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                 >  1:4349939‑M2/1‑94 (MQ=255)
                                                                     ccgtgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                 <  1:1417161‑M2/101‑8 (MQ=255)
                                                                     ccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGGTCGGTG.GGATCCCTGGCATAGGTCGGGTCCCTGGGAtcg                                 >  1:201337/1‑101 (MQ=14)
                                                                     ccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTC.GGATCCGTGGGAGCGGTGGGATCCGTGGGAtcg                                 >  2:2015548/1‑101 (MQ=11)
                                                                     ccgtcggatcggtcgggtccgtgggagcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTgggatcg                                 >  1:56075‑M2/1‑94 (MQ=255)
                                                                      cgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTC.GGGTCCGGGGGATCGGTCGGGTCCGTGGGagggg                                >  2:785057/1‑97 (MQ=11)
                                                                        tgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGATCCGTgggatcggtc                              >  2:2717741‑M2/1‑91 (MQ=255)
                                                                           gatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggg                           <  1:1189808‑M2/101‑14 (MQ=255)
                                                                            atcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTCGGATCGGTCGGGTCCGGGGGATCGgtcgggc                          >  1:1050274/1‑100 (MQ=11)
                                                                              cggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTG.GGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                        <  2:3983283/101‑1 (MQ=11)
                                                                                gtcgggtgcgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgt                      >  1:1561928‑M2/1‑83 (MQ=255)
                                                                                 tcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTG.GGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtc                     <  1:1284586/101‑1 (MQ=11)
                                                                                  cgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGGgggatcggtcgggaccgtgg                    <  2:1090701‑M2/101‑21 (MQ=255)
                                                                                   ggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTC.GGGTCCGTCGGATCGGTCGGATCCGTGGGATCGgtcggggccgtcgg                   >  2:3042264/4‑101 (MQ=11)
                                                                                    tatccgtggggtcggtccgATCGGTGGGGTCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggga                  <  2:1571360/99‑2 (MQ=11)
                                                                                    ggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggg                  <  1:2855068‑M2/101‑23 (MQ=255)
                                                                                     gtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggt                 <  1:2283620‑M2/101‑24 (MQ=255)
                                                                                     atccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggat                 <  1:1001929/100‑3 (MQ=255)
                                                                                      tccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGGGGGATCGGTC.GGATCCGTAGGATCGGTCGGGTCCGTAGGATGGgtcgggtccgtggggtc                >  1:721288/1‑101 (MQ=14)
                                                                                          taggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTg            >  1:970207/3‑101 (MQ=14)
                                                                                            ggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgg          <  1:3661202‑M2/101‑31 (MQ=255)
                                                                                             ggtccgtcggatCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcggg         <  1:3286315‑M2/99‑32 (MQ=255)
                                                                                             gatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTG.GGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggg         <  2:2555514‑M2/101‑32 (MQ=255)
                                                                                              gtccaccggATCCGGGGGGTCCGTCGGGTCCGTGGGATCAGTC.GGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAt        <  2:3187121/95‑1 (MQ=11)
                                                                                               tcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTTGGATCCGTGAGGTCCGTGGGaac       >  1:2966198/1‑99 (MQ=11)
                                                                                               tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGCCCGTGGGAtc       >  1:2203039/1‑101 (MQ=14)
                                                                                                cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtcc      <  1:2760080‑M2/101‑35 (MQ=255)
                                                                                                cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGGCCGTGGGAtcg      >  2:4164590/1‑101 (MQ=11)
                                                                                                ccgtcggTTCCGTGGGTTCGGTCGGATGCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtcc      <  2:3395885/99‑2 (MQ=11)
                                                                                                 ggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggtcggtcgggtccg     <  2:3930193‑M2/101‑36 (MQ=255)
                                                                                                 ggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccg     >  2:3585713/1‑99 (MQ=11)
                                                                                                   tcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggatccgtg   >  1:1112378‑M2/1‑64 (MQ=255)
                                                                                                   gcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtc   <  1:3154567/100‑1 (MQ=14)
                                                                                                    cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtgg  >  1:1344907/1‑99 (MQ=11)
                                                                                                    |                                                                                                     
TCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTC.GGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCG  >  NC_014034/221900‑222100

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG