Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,155,475 A→C 38.1% M14L (ATC→CTT)  RCAP_rcc02000 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,155,4750A→C38.1% 0.0 ‑0.6 ‑1.8295M14L (ATC→CTC) RCAP_rcc02000hypothetical protein
Reads supporting (aligned to +/- strand):  new base (66/46):  ref base (100/82):  total (167/128)
Fisher's exact test strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.68e-02

Only 100 of 304 total aligned reads displayed.

GCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTCCTGAGCCG  >  NC_014034/2155376‑2155575
                                                                                                    |                                                                                                    
gCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCAt                                                                                                     >  2:852211/1‑101 (MQ=255)
  ccTCGGGCCCCGTCCCCGTCCCCCCCCGCCGGGCGG.AAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCt                                                                                                   <  2:336875/101‑1 (MQ=255)
      gggCGCCGTCACCGTCACCGCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGAGCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTgggg                                                                                               >  2:4471971/1‑101 (MQ=255)
       ggCGCCGTCACCGTCACCGCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTggggc                                                                                              >  1:1148220/1‑101 (MQ=255)
            cgtcaccgtcacCGCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTggggcggggc                                                                                         >  2:2453817/1‑101 (MQ=255)
               caccgtcaccgCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGGCGt                                                                                      >  2:3212821/1‑101 (MQ=255)
                cccgcccccgcccgcGGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGATCTCGTCGATGCCATCTGGGGCGGGGCCGTc                                                                                     <  2:3548526/94‑1 (MQ=255)
                   gtcaccgCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacc                                                                                  >  2:2883947/1‑101 (MQ=255)
                    tccccgccggcCGGGCGG.AAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacca                                                                                 <  2:1460468/98‑1 (MQ=255)
                    cccccccccgcgGGGCGG.AAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacca                                                                                 <  2:3763986/95‑1 (MQ=255)
                    cccccccccgccGGGCGG.AAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcacca                                                                                 <  2:2694474/95‑1 (MQ=255)
                     caccgCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTcaccac                                                                                >  2:1967055/1‑101 (MQ=255)
                      accgCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACg                                                                               >  2:37024/1‑101 (MQ=255)
                       ccgcccggCGGGCGGAAAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACg                                                                               <  2:4211672/101‑1 (MQ=255)
                         ccccgcgGGGCGG.AAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATg                                                                            <  1:216993/100‑1 (MQ=255)
                           ccgcggGGCGG.AAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggc                                                                          <  1:2003607/101‑1 (MQ=255)
                           ccgcggGGCGG.AAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggc                                                                          <  2:1667257/101‑1 (MQ=255)
                            ccgcgGGCAG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGCCGATGCCATCGTGGGCGGGGCCGTCACCACGATggcg                                                                         <  2:4239053/98‑1 (MQ=255)
                             gccgGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcgg                                                                        >  1:1930573/1‑101 (MQ=255)
                              ccgGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcggc                                                                       >  2:2563993/1‑101 (MQ=255)
                              ccgGGCGA.AAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcggc                                                                       <  1:3523029/101‑1 (MQ=255)
                                 ggcgg.cAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATggcggcggc                                                                    >  2:2650484/1‑101 (MQ=255)
                                      aaaaGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTc                                                                <  2:3904860/99‑1 (MQ=255)
                                         cgCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGATCTCGTCGATGCCATCTGGGGTGGGGCCGTCACGACCAGGGCGGCGGGGCTCGcg                                                             >  2:1164835/2‑100 (MQ=255)
                                         aGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGcc                                                             <  1:3843560/101‑1 (MQ=255)
                                         aGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGcc                                                             >  2:3092071/1‑101 (MQ=255)
                                           tctcgAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCgg                                                           <  1:662955/98‑1 (MQ=255)
                                             gcgAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggc                                                         >  1:4255047/1‑101 (MQ=255)
                                             gcgAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggc                                                         <  2:1620904/101‑1 (MQ=255)
                                              cgAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCggaga                                                        >  2:662893/1‑98 (MQ=255)
                                              cgAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcg                                                        >  2:4129139/1‑101 (MQ=255)
                                               gAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcgg                                                       >  1:4063633/1‑101 (MQ=255)
                                                aaaTGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGggcggc                                                      <  2:2311288/101‑1 (MQ=255)
                                                 aaGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCt                                                     >  1:453798/1‑101 (MQ=255)
                                                   gggCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGa                                                   <  2:2362406/101‑1 (MQ=255)
                                                   ggagaaaaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGa                                                   <  1:2792180/91‑1 (MQ=255)
                                                    ggCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAt                                                  <  2:2123091/101‑1 (MQ=255)
                                                      cAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGAtgt                                                <  1:3315153/101‑1 (MQ=255)
                                                        aaaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCCGGCCGGGCGCCTGATTTgg                                              >  2:3688282/6‑101 (MQ=255)
                                                         gATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGc                                             >  1:2108490/1‑101 (MQ=255)
                                                         aaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGc                                             <  1:1334903/97‑1 (MQ=255)
                                                         aaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGc                                             <  2:3005858/97‑1 (MQ=255)
                                                         aaaaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGc                                             <  2:1288123/97‑1 (MQ=255)
                                                           aaaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCAt                                           >  1:2732377/3‑101 (MQ=255)
                                                            cATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCAtt                                          <  2:1723737/101‑1 (MQ=255)
                                                            cATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCAtt                                          <  2:3963653/101‑1 (MQ=255)
                                                            cATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCAtt                                          <  1:4472266/101‑1 (MQ=255)
                                                            aaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGTCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCAtt                                          <  2:567837/100‑1 (MQ=255)
                                                            aaTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCAtt                                          >  1:1674501/2‑101 (MQ=255)
                                                             aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                                         >  1:3131167/1‑101 (MQ=255)
                                                             aTGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTc                                         >  1:1523980/1‑101 (MQ=255)
                                                              tGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCg                                        >  2:4016974/1‑101 (MQ=255)
                                                              tGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCg                                        >  2:496776/1‑101 (MQ=255)
                                                              tGACCGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCg                                        >  1:4343139/1‑101 (MQ=255)
                                                                 ccGATACCAGAGGCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGAcc                                     <  1:3279275/100‑1 (MQ=255)
                                                                   gATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGa                                   <  1:1021790/101‑1 (MQ=255)
                                                                    ataccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAg                                  >  2:2768216/9‑101 (MQ=255)
                                                                    aTGCAAGGGGCATCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAg                                  >  2:1134603/1‑101 (MQ=255)
                                                                      gCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAggg                                >  2:1866100/1‑100 (MQ=255)
                                                                      accagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAggg                                >  1:3739457/7‑100 (MQ=255)
                                                                       ccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTg                               <  2:3126863/96‑1 (MQ=255)
                                                                       ccagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTg                               <  2:4365096/96‑1 (MQ=255)
                                                                       cAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTg                               >  1:3617644/1‑101 (MQ=255)
                                                                        cagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGc                              >  2:3510633/5‑101 (MQ=255)
                                                                        cagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGc                              <  2:2967561/97‑1 (MQ=255)
                                                                        cagaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGc                              <  1:7978/97‑1 (MQ=255)
                                                                         agaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCg                             >  1:1236447/4‑101 (MQ=255)
                                                                         agaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCg                             >  1:3678146/4‑101 (MQ=255)
                                                                          ggggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGGGGTGCgg                            >  2:1560454/1‑101 (MQ=255)
                                                                          ggggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCgg                            <  2:110903/101‑1 (MQ=255)
                                                                          ggggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCgg                            >  1:1227476/1‑101 (MQ=255)
                                                                          gaggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCgg                            <  1:2429926/99‑1 (MQ=255)
                                                                           aggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCggg                           <  2:905341/100‑1 (MQ=255)
                                                                            ggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTTGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgc                          >  1:565205/1‑101 (MQ=255)
                                                                            ggCCTCGTCGAGGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgc                          >  1:3815100/1‑101 (MQ=255)
                                                                            ggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGgc                          >  2:1791048/1‑101 (MQ=255)
                                                                            ggCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGAGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGgc                          >  1:3925617/1‑101 (MQ=255)
                                                                                  gtcgagGTCATCGATGCCCTTTGTGGTGGCGCCGTCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcg                    <  1:4013444/98‑4 (MQ=255)
                                                                                  gtcgagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcg                    >  2:2245803/4‑98 (MQ=255)
                                                                                  gtcgagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGgcgcggcg                    <  2:2393872/98‑4 (MQ=255)
                                                                                   tAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAg                   >  1:2272694/1‑101 (MQ=255)
                                                                                    agagCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAgc                  <  2:1224332/101‑1 (MQ=255)
                                                                                     gagGTCATCGATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGGGCGGgcgcggcggcg                 >  2:2030702/1‑95 (MQ=255)
                                                                                     gagCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAgcg                 <  2:455465/101‑1 (MQ=255)
                                                                                        cTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGcgccg              <  2:3845137/101‑1 (MQ=255)
                                                                                         tcgtcgATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGcgccgc             >  1:2628629/1‑101 (MQ=255)
                                                                                         tcatcgATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGGCGCGGCGGcgccgg             >  1:1709558/4‑100 (MQ=255)
                                                                                           atcgATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGGCGCGGCGGCGCCGGtt           <  1:2572169/100‑1 (MQ=255)
                                                                                            acaaTGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTGGACCGAGGTGCGGGCGCGGCGGCGCCGGttt          <  1:4221510/98‑2 (MQ=255)
                                                                                             cgATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTcc         >  1:2899505/1‑101 (MQ=255)
                                                                                             cgATGCCATCTGGGGCGGGGCCGTCACAACGATGGCGGCGGCGCTAGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTcc         >  1:4591997/1‑101 (MQ=255)
                                                                                              gATGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGGCGCGGCGGCGCCGGttttt        >  1:475624/1‑98 (MQ=255)
                                                                                              gATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTCCt        <  2:4593576/101‑1 (MQ=255)
                                                                                               aTGCCCTTTGTGGTGGCGCCATCACGACGATGGCGGCGGCGCTGGCCGGGCGCCTGATGTGGCATTCGACCGAGGTGCGGGCGCGGCGGCGCCGGtttttc       >  1:2302082/1‑97 (MQ=255)
                                                                                                 gCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTCCTGAg     >  2:1380571/1‑101 (MQ=255)
                                                                                                 gCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTCCTCAg     >  2:546143/1‑101 (MQ=255)
                                                                                                  ccATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTCCTGAGc    >  2:363691/1‑101 (MQ=255)
                                                                                                   cATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTCCTGAGcc   <  1:2987210/101‑1 (MQ=255)
                                                                                                   cATCTGGGGCGCGGCCGTAACCGCGATGGCGCCGGCTCACGCCGTGCGTCTGATGTGGCTTTCGGCCGAGGTGCGGGCGCGGAAGCGCGGCTTCCTTAGcc   <  2:770787/101‑1 (MQ=255)
                                                                                                    aTCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTCCTGAGCcg  <  1:3092315/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GCCCTCGGGCGCCGTCACCGTCACCGCCCGCCGGGCGG.CAAGCCGCGAAAGGGCAGGATCATGAGCGATGCAAGGGGCCTCGTAGAGCTCGTCGATGCCATCTGGGGCGGGGCCGTCACCACGATGGCGGCGGCGCTCGCCGGGCGGCTGATGTGGCATTCGACCGAGGTGCGGGCGCGCAAGCGCCGCTTCCTGAGCCG  >  NC_014034/2155376‑2155575

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG