Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 207,483 A→C 15.9% V786V (GTA→GTC dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913207,4830AC15.9% 246.3 / 39.5 106V786V (GTA→GTCdnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (48/40);  new base C (6/11);  total (54/52)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-02

Only 100 of 109 total aligned reads displayed.

TTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGT  >  NC_000913/207384‑207582
                                                                                                   |                                                                                                   
ttACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTa                                                                                                     <  1:2439791/100‑1 (MQ=255)
   cGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTaa                                                                                                    >  2:641930/1‑98 (MQ=255)
    ggATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGa                                                                                                 <  2:1905309/100‑1 (MQ=255)
    ggATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGa                                                                                                 >  1:2845684/1‑100 (MQ=255)
     gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                                <  1:195096/100‑1 (MQ=255)
     gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                                >  1:2723022/1‑100 (MQ=255)
     gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAc                                                                                                >  1:862418/1‑100 (MQ=255)
          aaCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgccg                                                                                           >  1:477556/1‑100 (MQ=255)
          aaCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgcc                                                                                            >  1:1255230/1‑99 (MQ=255)
           aCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGCCTGAAAGCGCACTATCCTGCGGAGTTTATAACGGCGGTAATGACCACCGa                                                                                          >  1:1232853/1‑100 (MQ=255)
            cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGTGGTCATTAccgccgcc                                                                                         <  2:1041168/100‑3 (MQ=255)
            cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                         <  1:2031048/100‑3 (MQ=255)
            cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                         <  1:2491973/100‑3 (MQ=255)
            cAAATAGCACTCTGCGGCCTATGGTGTAATGTCATAGCAAACGTTATTGCTTAAAGTGCATTACCCTGCGGAGTTTATGGCGGCGGTCATTAccgacgca                                                                                         <  2:2035386/100‑3 (MQ=255)
             aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcca                                                                                        <  1:1177579/100‑4 (MQ=255)
             aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                         <  2:2153007/99‑1 (MQ=255)
              aaTCGCACTCTGCGGCCTATGTTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                         <  1:1809262/98‑3 (MQ=255)
              aaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                       <  2:423825/100‑5 (MQ=255)
              aaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGata                                                                                        <  2:2524924/99‑1 (MQ=255)
               aTCGTACTCTGCGGCCTATGCTTTTTTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgtcgccata                                                                                      <  2:1637916/100‑6 (MQ=255)
               aTCGCACTCTGGGGTCTATGCTTTGTTGTCATATCAAACGTTATGGCTGAAAGCGCATTATCCTGCTGATTTTATGGCGGCGGTAATGAccgccg                                                                                           <  2:2569921/95‑1 (MQ=255)
               aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATg                                                                                      >  2:2404413/1‑100 (MQ=255)
               aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATg                                                                                      <  1:6145/100‑1 (MQ=255)
               aTCGCACCCTGGGGCCTATGCTTTTGTGTCATATCAAACGTTATGGCTGAAAGCGCACTACCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                       <  2:2945883/99‑5 (MQ=255)
                tCGCACTCTGCTGCTTATTCTTTGGTTGCATATCAAACGTTATGGCTGAAAGCGCACTATCCAGCGGAGTTGATGGCGGCGGTAATGACCGCCGATATgg                                                                                     <  1:1882381/100‑1 (MQ=255)
                tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATgg                                                                                     <  1:2601359/100‑1 (MQ=255)
                 cGCACTCTGCGGCCTATGCCTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGa                                                                                    <  1:2678367/100‑1 (MQ=255)
                  gCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATAGCGGCGGTAATGACCGCCGATATGGac                                                                                   >  1:1329739/1‑100 (MQ=255)
                   cACTCTGCGGCCTATGCTTTGGTGTTATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                  <  1:663434/100‑1 (MQ=255)
                   cACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                  <  1:693356/100‑1 (MQ=255)
                    aCTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccataaactc                                                                                 <  2:2191897/100‑11 (MQ=38)
                     ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaac                                                                                >  1:1530242/1‑100 (MQ=255)
                     ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaac                                                                                <  1:2421591/100‑1 (MQ=255)
                       ctGCGGCTTATGCTTTGGTGTCATATCAAACTTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccataaac                                                                                   <  1:1176104/95‑9 (MQ=37)
                       ctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                              >  1:1946425/1‑100 (MQ=255)
                       ctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                              >  1:2166904/1‑100 (MQ=255)
                         gCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCg                                                                            <  1:426224/100‑1 (MQ=255)
                            gCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCAGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCgag                                                                          >  2:2490885/1‑99 (MQ=255)
                             ccTATGCTTTGGTGTCGTATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGaa                                                                        >  1:1006159/1‑100 (MQ=255)
                                aTGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCTGAGTTTATGGCGGCGGAAATGACCGCCGATATGGACAACACCGAGAAgg                                                                      >  1:2694658/1‑99 (MQ=255)
                                     ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTGATGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                                <  2:1610412/100‑1 (MQ=255)
                                     ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                                <  2:1277056/100‑1 (MQ=255)
                                      tGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTggg                                                               <  2:2166881/100‑1 (MQ=255)
                                       ggTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATAACCGCCGATATGGACAACAACGAGAAGGTGGTGGGt                                                              >  1:1275797/1‑100 (MQ=255)
                                        gtgtGATATCAAGCTTCATGGCTGAAAGCGCACTATCGTGCTGAGTTTAGGGCGGCGTTAGTGACCGCCGATATGGACAACACCGAGAAGGTGGTggg                                                               <  2:2824462/98‑1 (MQ=255)
                                         tgtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCt                                                            >  1:2620430/1‑100 (MQ=255)
                                          gtCATATGAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtg                                                           >  1:1593817/1‑100 (MQ=255)
                                          gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGAAGATATGGACAACACCGAGAAGGTGGTGGGTCtg                                                           >  1:1752399/1‑100 (MQ=255)
                                          gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGAGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtg                                                           >  1:1038663/1‑100 (MQ=255)
                                            cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACAGCGAGAAGGTGGTGGGtat                                                            >  2:1386325/1‑95 (MQ=255)
                                            cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggt                                                         >  2:271561/1‑100 (MQ=255)
                                              tatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                       <  2:103431/100‑1 (MQ=255)
                                              tatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                       >  1:1544814/1‑100 (MQ=255)
                                              tatCAAACGTTATGGCTGAAAGCGCACTAGCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGtggtgg                                                                >  2:664631/1‑91 (MQ=255)
                                               atCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGa                                                      >  1:1298521/1‑100 (MQ=255)
                                               atCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCAATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGa                                                      >  1:85889/1‑100 (MQ=255)
                                                tCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATAGTGACAACACCGAGAAGGTAGTGGGTGTTGTGGAt                                                     >  1:1313975/1‑100 (MQ=255)
                                                   aaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAg                                                  <  1:2138835/100‑1 (MQ=255)
                                                        cgatcttGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTg                                             <  2:1732717/94‑1 (MQ=255)
                                                         aTGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTgg                                            >  1:1338246/1‑100 (MQ=255)
                                                          tGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGc                                           >  2:22984/1‑100 (MQ=255)
                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGTGTCTGGTGGATGAGTGCTGGCg                                          <  2:2967187/100‑1 (MQ=255)
                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGTTGAGTGCTAGCg                                          >  2:290374/1‑100 (MQ=255)
                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCg                                          <  2:1329177/100‑1 (MQ=255)
                                                              tGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTgagatcg                                       >  1:1732737/1‑94 (MQ=255)
                                                               gAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATg                                      >  1:401884/1‑100 (MQ=255)
                                                                aaaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgg                                     >  1:454727/1‑100 (MQ=255)
                                                                 aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                    <  1:1661202/100‑1 (MQ=255)
                                                                 aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                    <  1:2496760/100‑1 (MQ=255)
                                                                   gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                  <  2:3008100/100‑1 (MQ=255)
                                                                    ctcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCt                                 <  2:53382/98‑1 (MQ=255)
                                                                     gcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTg                                >  2:1354612/1‑100 (MQ=255)
                                                                     gcACTATCCTGCGGAGATTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGCTGAGAGCTGGCGGACGGGGCTg                                >  2:1483154/1‑100 (MQ=255)
                                                                      cACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAAAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCtag                               >  1:53382/1‑98 (MQ=255)
                                                                          aTCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAAt                           <  1:2217008/100‑1 (MQ=255)
                                                                           tCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATc                          <  1:1098289/100‑1 (MQ=255)
                                                                                cGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGcc                     <  1:2035900/100‑1 (MQ=255)
                                                                                   aGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacc                  >  1:1656456/1‑100 (MQ=255)
                                                                                    gTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGGTGAAAATTCTGcca                    >  2:3039016/1‑97 (MQ=255)
                                                                                    gTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacca                 >  1:2679084/1‑100 (MQ=255)
                                                                                     tttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAg                <  2:401881/100‑1 (MQ=255)
                                                                                     tttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAg                <  1:271564/100‑1 (MQ=255)
                                                                                      ttATTGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGa               <  2:636159/100‑1 (MQ=255)
                                                                                      ttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGa               >  1:173649/1‑100 (MQ=255)
                                                                                      ccatatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGTTGGCGGATGGGGCTGAAAATCCTGCCACCAGa               >  1:2027730/7‑100 (MQ=255)
                                                                                       tATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat              <  2:1635002/100‑1 (MQ=255)
                                                                                       tATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat              >  1:1676656/1‑100 (MQ=255)
                                                                                       catatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat              >  2:2703353/6‑100 (MQ=255)
                                                                                       catatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGat              >  1:3048510/6‑100 (MQ=255)
                                                                                        atatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTTGCGGATGGGGCTGAAAATCCTGCCACCAGata             >  1:2153350/5‑100 (MQ=255)
                                                                                        atatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGata             >  1:1946293/5‑100 (MQ=255)
                                                                                         tatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTAGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGatat            >  1:1513139/4‑100 (MQ=255)
                                                                                          ggcggcGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAtata           >  1:2843775/1‑100 (MQ=255)
                                                                                                ggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTcc     <  1:995033/100‑1 (MQ=255)
                                                                                                ggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTACTAGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTcc     >  1:767713/1‑100 (MQ=255)
                                                                                                  tAATGACTGCTGATATGGATAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGTCATCAGatat            <  2:2607747/91‑1 (MQ=255)
                                                                                                  tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg   <  2:2749885/100‑1 (MQ=255)
                                                                                                  tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg   <  2:66568/100‑1 (MQ=255)
                                                                                                  tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTAGATGAGTGCTAACGAATGGGACTGAAAATCCTGCCACCAGATATAAACTcc     >  1:851682/1‑98 (MQ=255)
                                                                                                   aaTGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGt  <  1:664635/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGT  >  NC_000913/207384‑207582

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: