Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,189,524 Δ1 bp 5.2% coding (4175/4224 nt) rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,189,5240G. 5.2% 1137.7 / 55.8 308coding (4175/4224 nt)rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (150/142);  new base . (10/6);  total (160/148)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 309 total aligned reads displayed.

TGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGA  >  NC_000913/4189380‑4189673
                                                                                                                                                  |                                                                                                                                                     
tGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTgg                                                                                                                                                  >  1:1495393/1‑150 (MQ=255)
     ggacTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCTTATGAGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGaa                                                                                                                                             >  1:1442006/4‑150 (MQ=255)
      gCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAAc                                                                                                                                            >  1:3631434/1‑150 (MQ=255)
         tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgc                                                                                                                                         >  1:2492461/1‑150 (MQ=255)
           aaaGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgcgg                                                                                                                                       >  1:3417731/1‑148 (MQ=255)
              gagaACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAAc                                                                                                                                    >  2:2440914/1‑150 (MQ=255)
               agaACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACg                                                                                                                                   <  1:3856846/150‑1 (MQ=255)
               agaACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACg                                                                                                                                   >  1:1282608/1‑150 (MQ=255)
                 aaCGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCa                                                                                                                                 <  2:603412/150‑1 (MQ=255)
                    gTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGt                                                                                                                              >  1:2865442/1‑150 (MQ=255)
                     ttATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTc                                                                                                                             <  1:620398/150‑1 (MQ=255)
                      tATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCt                                                                                                                            <  1:1117352/150‑1 (MQ=255)
                          gTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                                                                                                        >  2:2543575/1‑150 (MQ=255)
                            gggTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCgg                                                                                                                      <  1:788380/150‑1 (MQ=255)
                              gtcgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGtt                                                                                                                    >  2:3936039/1‑150 (MQ=255)
                                tgggtcgtctgatcccggcaggtaccggttaACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCt                                                                                                                  <  1:4093506‑M2/119‑1 (MQ=255)
                                cgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCt                                                                                                                  <  2:859668/150‑1 (MQ=255)
                                 gtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTg                                                                                                                 <  1:2207745/150‑1 (MQ=255)
                                  tcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                                                                                                <  1:120244/150‑1 (MQ=255)
                                  tcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                                                                                                <  1:988460/150‑1 (MQ=255)
                                    tGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATa                                                                                                              >  2:2132656/1‑150 (MQ=255)
                                    tGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATa                                                                                                              >  1:3492261/1‑150 (MQ=255)
                                      aTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAAc                                                                                                            >  1:1065563/1‑150 (MQ=255)
                                        cccGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGa                                                                                                          >  1:1533675/1‑150 (MQ=255)
                                         ccGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAg                                                                                                         >  2:3088904/1‑150 (MQ=255)
                                          cGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGt                                                                                                        >  2:1432921/1‑150 (MQ=255)
                                              aGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATc                                                                                                    <  1:2739075/150‑1 (MQ=255)
                                                gTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGt                                                                                                  <  1:1157829/150‑1 (MQ=255)
                                                    cGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAAt                                                                                              >  2:3418449/1‑150 (MQ=255)
                                                    cGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAAt                                                                                              >  2:1495495/1‑150 (MQ=255)
                                                      gTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATcc                                                                                            <  1:2047224/150‑1 (MQ=255)
                                                       ttACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCg                                                                                           >  1:3665658/1‑150 (MQ=255)
                                                           gcgTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCaaa                                                                                       >  2:2887212/1‑150 (MQ=255)
                                                             gTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCACCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATaa                                                                                    <  1:3945635/150‑1 (MQ=255)
                                                                ccaccaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACg                                                                                  <  1:3753354/150‑1 (MQ=255)
                                                                 caccaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt                                                                                 >  1:1424431/1‑150 (MQ=255)
                                                                 caccaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt                                                                                 >  1:1980725/1‑150 (MQ=255)
                                                                  accaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTa                                                                                <  2:2334510/150‑1 (MQ=255)
                                                                  accaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTa                                                                                >  2:759594/1‑150 (MQ=255)
                                                                  accaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTa                                                                                >  1:2592307/1‑150 (MQ=255)
                                                                  accaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTa                                                                                <  2:2383561/150‑1 (MQ=255)
                                                                   ctaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTaa                                                                               <  1:497499/148‑1 (MQ=255)
                                                                    caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCA‑‑‑‑GGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAcc                                                                          >  1:1508826/1‑150 (MQ=255)
                                                                    caGGATCGTACGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCA‑‑‑‑GGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAcc                                                                          >  2:3201615/1‑150 (MQ=255)
                                                                      ggATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTaaaaa                                                                            <  1:1559759/150‑1 (MQ=255)
                                                                         tCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAccc                                                                         <  1:3610159/150‑1 (MQ=255)
                                                                          cGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCGGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCg                                                                        <  1:271552/150‑1 (MQ=255)
                                                                          cGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCGGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCg                                                                        <  1:1327139/150‑1 (MQ=255)
                                                                           gTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAccc                                                                         <  2:677323/150‑1 (MQ=255)
                                                                           gTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAccc                                                                         <  2:1753811/150‑1 (MQ=255)
                                                                             aTGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCAGGTAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCtt                                                                     <  1:1021967/150‑1 (MQ=255)
                                                                             aTGCGTCGCAGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCtt                                                                     >  1:3378281/1‑150 (MQ=255)
                                                                              tGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTc                                                                    >  1:2706840/1‑150 (MQ=255)
                                                                               gCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcg                                                                   >  1:3682295/1‑150 (MQ=255)
                                                                                 gTCGCCGTGCTGCGGGTTCAGCTCCGGGTGCACCGCT‑‑GGTGACTGCAGATGACGCATCTGCGTGCCTGGGAGAACTTCTGAACTCAGGTCTGGGCTGTTCTGATAACGTGTAATCGTTAATCCGCAAATAAAGTAAAAACCCGCTTcggc                                                                 <  2:2358592/150‑1 (MQ=255)
                                                                                  tCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggcg                                                                <  1:1693543/150‑1 (MQ=255)
                                                                                       gTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttt                                                           <  1:3002756/150‑1 (MQ=255)
                                                                                       gTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttt                                                           <  1:71162/150‑1 (MQ=255)
                                                                                          ctgcGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGtttttt                                                        >  1:2842350/1‑150 (MQ=255)
                                                                                           tgcGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAATCCCGCTTCGGCGGGttttttt                                                       >  2:1181429/1‑150 (MQ=255)
                                                                                           tgcGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttttttt                                                       >  2:2802638/1‑150 (MQ=255)
                                                                                           tgcGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttttttt                                                       >  1:2338132/1‑150 (MQ=255)
                                                                                           tgcGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCACCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTa                                                      >  2:3358614/1‑150 (MQ=255)
                                                                                            gcGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTACTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTa                                                      <  2:3485250/150‑1 (MQ=255)
                                                                                               ggTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAATACCCGCTTCGGCGGGTTTTTTTATgg                                                   >  1:1795694/1‑150 (MQ=255)
                                                                                                gTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATggg                                                  >  2:1325854/1‑150 (MQ=255)
                                                                                                 tGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgggg                                                 <  2:523305/150‑1 (MQ=255)
                                                                                                     gCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAg                                             >  1:566651/1‑150 (MQ=255)
                                                                                                     gCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAg                                             >  2:1510350/1‑150 (MQ=255)
                                                                                                     gCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAg                                             >  1:456064/1‑150 (MQ=255)
                                                                                                      cTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGt                                            >  2:3964438/1‑150 (MQ=255)
                                                                                                         cGGCTGCACCGCA‑‑GGTGACTGCAGACGACGCATCTGCCCGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGGACACACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa                                         <  2:200657/150‑1 (MQ=255)
                                                                                                          ggCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATCACGTAAAAACCCGCTTCGGCGGGTTTTTTTAGGGGGGGAGTTTAg                                        >  2:2666655/1‑150 (MQ=255)
                                                                                                           gCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAgg                                       >  1:3321144/1‑150 (MQ=255)
                                                                                                               cACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGaaa                                   <  2:1632345/150‑1 (MQ=255)
                                                                                                               cACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGaaa                                   <  2:2213813/150‑1 (MQ=255)
                                                                                                                     a‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAt                             <  2:4201709/150‑1 (MQ=255)
                                                                                                                         gTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAttt                           >  1:2258772/1‑150 (MQ=255)
                                                                                                                              tGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAg                      >  2:776953/1‑150 (MQ=255)
                                                                                                                               gCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGa                     >  1:4156559/1‑150 (MQ=255)
                                                                                                                                cAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGaa                    <  1:1380376/150‑1 (MQ=255)
                                                                                                                                cAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGGAATCGGTAATCCGCAAATACCGTAAAAACCCGCTTCGGCGGGTGTGTTTAGGGGGGGAGTGGAGGGAAAGAGCATTTGTCAGaa                    >  2:1053865/1‑150 (MQ=255)
                                                                                                                                 agaagaCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATCACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGGTTAGGGAAAGAGCATTTGTCAGAat                   >  2:3752410/1‑150 (MQ=255)
                                                                                                                                  gaagaCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCGTCGGCGGGTTTTTGTCTGGGGGGAGGGGAGGGAAAGAGCAGTTGTCAGactc                  >  2:1571802/1‑147 (MQ=255)
                                                                                                                                     gaCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATAttt               <  1:303620/150‑1 (MQ=255)
                                                                                                                                     gaCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATAttt               <  1:2650084/150‑1 (MQ=255)
                                                                                                                                     gaCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATAttt               <  2:3717786/150‑1 (MQ=255)
                                                                                                                                      aCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTa              >  2:2609213/1‑150 (MQ=255)
                                                                                                                                       cGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa             <  1:3616194/150‑1 (MQ=255)
                                                                                                                                       cGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa             <  1:2543495/150‑1 (MQ=255)
                                                                                                                                        gCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAg            >  1:3273544/1‑150 (MQ=255)
                                                                                                                                        gCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTGTAGGGAAAGAGCATTTGTCAGAATATTTAAg            >  1:1074770/1‑150 (MQ=255)
                                                                                                                                           tCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGaa         <  2:3612044/150‑1 (MQ=255)
                                                                                                                                             tGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAAtt       >  2:2572321/1‑150 (MQ=255)
                                                                                                                                               ccAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc     >  1:4046264/1‑150 (MQ=255)
                                                                                                                                               ccAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc     >  1:3994016/1‑150 (MQ=255)
                                                                                                                                               ccAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc     >  1:2143553/1‑150 (MQ=255)
                                                                                                                                               ccAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc     >  1:3145822/1‑150 (MQ=255)
                                                                                                                                               ccAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGGTTAGGGAAAGAGCATTTGTCAGACTATTTAAGGAATTTc     >  2:3007141/1‑150 (MQ=255)
                                                                                                                                                cAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCt    <  2:430438/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGA  >  NC_000913/4189380‑4189673

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