Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,189,448 C→T 22.7% Q1367* (CAG→TAG)  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,4480CT22.7% 254.9 / 72.4 132Q1367* (CAG→TAG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (48/54);  new base T (17/13);  total (65/67)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

Only 100 of 139 total aligned reads displayed.

CCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAA  >  NC_000913/4189302‑4189596
                                                                                                                                                  |                                                                                                                                                      
  ttCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGAt                                                                                                                                                   >  2:1504730/1‑150 (MQ=255)
    cATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCg                                                                                                                                                 <  1:749874/150‑1 (MQ=255)
     aTCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGt                                                                                                                                                >  1:1279877/1‑150 (MQ=255)
      tgtcCGCGGCATCGTTCCAGGAGACCACTCGCGGGCTGACCGAAGCAGCCGTTGCGGGCATACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTa                                                                                                                                               <  1:372976/148‑1 (MQ=255)
      tctcCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTa                                                                                                                                               <  1:1540721/150‑1 (MQ=255)
      tctcCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTa                                                                                                                                               >  2:72862/1‑150 (MQ=255)
      tctcCGCGGCATCGTTCCAGGAGACCACTAGCGGGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGTTCGTCTGATCCCGGCAGGTACCTGTTACGCGTACCACTTGGATCGTa                                                                                                                                               <  2:724441/150‑1 (MQ=255)
      tctcCGCGGCAGCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTCTCGTGGGTCGTCGGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTa                                                                                                                                               <  2:328350/150‑1 (MQ=255)
      actcCGCGGCATCGTGCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTTTCGTGGGTCGTCTGATCCCGGCAGGGACCGGTTACGCGTACCACCAGGATCGTa                                                                                                                                               <  2:355510/149‑1 (MQ=255)
       ctcCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGAGCGTAt                                                                                                                                              >  2:1321916/1‑150 (MQ=255)
        tcCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTCATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATg                                                                                                                                             >  2:591940/1‑150 (MQ=255)
          cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCg                                                                                                                                           <  1:398425/150‑1 (MQ=255)
          cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTCGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATACTATGCg                                                                                                                                           >  2:495068/1‑150 (MQ=255)
           gcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGt                                                                                                                                          >  2:835537/1‑150 (MQ=255)
            cgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGAAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTc                                                                                                                                         >  1:1436513/1‑150 (MQ=255)
             tccgcggcatcgttccaggagaccactcgcgtgctgaccgaagcagccgttgcgggcaaacgcgacgaactgcgcggcctgaaagagaacgttatcgtggttcgtctgatcccggcaggtaccggttaACCACCAGGATCGTATGCGTCg                                                                                                                                        >  2:1242940‑M2/129‑150 (MQ=255)
               cATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGcc                                                                                                                                      >  1:583228/1‑150 (MQ=255)
                    ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGtgct                                                                                                                                 <  1:555118/150‑1 (MQ=255)
                     tCCCGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGtgctg                                                                                                                                <  2:1892976/150‑1 (MQ=255)
                      ccAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTAGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGtgctgc                                                                                                                               >  2:986450/1‑150 (MQ=255)
                       cAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCg                                                                                                                              >  1:804087/1‑150 (MQ=255)
                       cAGGAGACCACTCGCGTGCTGACAGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCg                                                                                                                              >  1:508930/1‑150 (MQ=255)
                         ggAGACCACTCGCGTGCTGACCGAAGCAGCCGATGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCggg                                                                                                                            <  2:1204118/150‑1 (MQ=255)
                            gaCCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGTAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGa                                                                                                                         <  1:978189/150‑1 (MQ=255)
                             aggagaccactcgcgtgctgaccgaagcagccgttgcgggcaaacgcgacgaactgcgcggcctgaaagagaacgttatcgtgggtcgtctgatcccggcaggtaccggttaACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGaa                                                                                                                        >  1:872620‑M2/113‑150 (MQ=255)
                                aCTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCt                                                                                                                     >  2:596705/1‑150 (MQ=255)
                                 cTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTc                                                                                                                    >  1:1871871/1‑150 (MQ=255)
                                   cgcgTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCg                                                                                                                  >  2:949282/1‑150 (MQ=255)
                                      gTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCTGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCt                                                                                                               <  2:1702635/150‑1 (MQ=255)
                                      gTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCt                                                                                                               <  2:444768/150‑1 (MQ=255)
                                       tGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTTAAGCTCCGGCTg                                                                                                              >  1:1286683/1‑150 (MQ=255)
                                       tGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTg                                                                                                              >  2:395183/1‑150 (MQ=255)
                                       ggCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTg                                                                                                              <  2:606431/149‑1 (MQ=255)
                                        gCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGc                                                                                                             >  2:1872107/1‑150 (MQ=255)
                                        gCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGc                                                                                                             <  1:1061471/150‑1 (MQ=255)
                                        gCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGc                                                                                                             >  1:1386827/1‑150 (MQ=255)
                                         cTGCCCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCa                                                                                                            <  1:1826173/150‑1 (MQ=255)
                                           gACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGGGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCCCCAGGATCGTATGCGTCGCCGTGCTGCGGGTGCAGCTCCGGCTGCAcc                                                                                                          >  1:113182/1‑150 (MQ=255)
                                             ccGCAGCAGCCGCTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGc                                                                                                        <  2:1783818/150‑1 (MQ=255)
                                              cGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCa                                                                                                       >  2:954002/1‑150 (MQ=255)
                                              cGAAGCAGCCGTTGCGGGAAAACGCGACGCACTGCGCGGCCTGAAAGAGAACGTTCTCGTGGGTCGGCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCa                                                                                                       <  1:1581709/150‑1 (MQ=255)
                                                aaGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGGCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑gg                                                                                                   <  2:1244791/150‑1 (MQ=255)
                                                 cgcagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGt                                                                                                  <  1:1105820/149‑1 (MQ=255)
                                                 agcagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGt                                                                                                  >  1:1338658/1‑150 (MQ=255)
                                                  gcagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTg                                                                                                 >  1:1334108/1‑150 (MQ=255)
                                                   cagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGa                                                                                                >  2:1771287/1‑150 (MQ=255)
                                                    agcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGTGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGAGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGAc                                                                                               >  2:402653/1‑150 (MQ=255)
                                                     gcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCCCCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACt                                                                                              >  1:624693/1‑150 (MQ=255)
                                                     gcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGCCt                                                                                              >  2:1406014/1‑150 (MQ=255)
                                                     gcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACt                                                                                              >  2:1806558/1‑150 (MQ=255)
                                                     gcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACt                                                                                              >  2:1413918/1‑150 (MQ=255)
                                                             gggCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCagaaga                                                                                      >  2:1465429/1‑150 (MQ=255)
                                                              ggCAAACGCGCCGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCGCCGGCTGCACCGCA‑‑GGTGACTGCagaagcc                                                                                     >  2:1789438/1‑148 (MQ=255)
                                                              ggCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGGGTGACTGCagaag                                                                                       >  1:1786917/1‑150 (MQ=255)
                                                                cAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGCTCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCTCCGCA‑‑GGTGCCTGCAGAAGACGc                                                                                   >  1:31037/1‑150 (MQ=255)
                                                                 aaaCGCGACGACCTGCGCGGCCTGAACGAGACCGTTATCGGGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCa                                                                                  <  1:1172984/150‑1 (MQ=255)
                                                                      cgacgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGc                                                                             <  1:1737546/150‑1 (MQ=255)
                                                                      cgacgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGAACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGCCGCATCTGc                                                                             >  1:768979/1‑150 (MQ=255)
                                                                      cgacgaACTGCGCGGACTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGc                                                                             >  1:1387535/1‑150 (MQ=255)
                                                                         cgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAg                                                                          >  2:1106659/1‑150 (MQ=255)
                                                                            aCTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGCCGCATCTGCCAGCCt                                                                       >  1:250899/1‑150 (MQ=255)
                                                                             cTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTg                                                                      <  2:249910/150‑1 (MQ=255)
                                                                              tGCGCGGCCTGAAAGAGAACTTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGGGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCAAGCCTgg                                                                     <  2:488743/150‑1 (MQ=255)
                                                                              tGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTgg                                                                     <  1:609169/150‑1 (MQ=255)
                                                                                cgcgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCa                                                                   <  2:181699/150‑1 (MQ=255)
                                                                                cgcgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCa                                                                   <  2:1096058/150‑1 (MQ=255)
                                                                                  cgGCCTGACAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGGTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCAGGCAGa                                                                 <  1:491407/150‑1 (MQ=255)
                                                                                    gCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAAc                                                               >  2:1494533/1‑150 (MQ=255)
                                                                                        gAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGTA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgct                                                           <  2:1336027/150‑1 (MQ=255)
                                                                                         aaaGAGAACGTTATCGTGGGTCATCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgctg                                                          <  2:1890182/150‑1 (MQ=255)
                                                                                         aaaGAGAACGTTAGCGTGGGTCGCCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGGCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgctg                                                          <  2:379214/150‑1 (MQ=255)
                                                                                                aCGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCGTGGCAGAACTGCTGACCGCAg                                                   >  1:610596/1‑150 (MQ=255)
                                                                                                 cGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGGGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAgg                                                  <  1:180426/150‑1 (MQ=255)
                                                                                                 cGTTATCGGGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAgg                                                  <  2:822535/150‑1 (MQ=255)
                                                                                                           gggcggcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGGGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGt                                        <  2:1383328/144‑1 (MQ=255)
                                                                                                           ggTCGTCTGCTCCCGGCAGGTACCGGTAACGCGTACCACCAGGATCGTATGCGTCGCCGTGCGGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGt                                        <  1:696598/150‑1 (MQ=255)
                                                                                                           ggTCGGCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGGCGCCGTGCTGCGGGTGAAGCTCCGGCGGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGt                                        <  1:246751/150‑1 (MQ=255)
                                                                                                             tcgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTc                                      >  2:883012/1‑150 (MQ=255)
                                                                                                              cgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGGCGCCGTGCTGCGGGGGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCt                                     <  1:1005276/150‑1 (MQ=255)
                                                                                                               gtgtGATCCCGGCAGGTACCGGTTTCGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTg                                    <  2:448620/147‑1 (MQ=255)
                                                                                                                tcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                   >  2:987616/1‑150 (MQ=255)
                                                                                                                   gATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATaa                                >  1:1238037/1‑150 (MQ=255)
                                                                                                                    aTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAAc                               >  2:1145911/1‑150 (MQ=255)
                                                                                                                        gggCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGGCGCCGTGCTGCGGGTGAAGCGCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGt                           <  1:495068/149‑1 (MQ=255)
                                                                                                                        cGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGt                           >  1:973191/1‑150 (MQ=255)
                                                                                                                               tACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGtt                    >  2:1230434/1‑150 (MQ=255)
                                                                                                                                    gTTACGCGTACCACTAGGATCGTATGCGTCGCAGTGCTGCGGGTGAAGCTCAGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCCGCCTGGCAGCACTGCTGCACGCAGGTATGGGCGGTTCTGATAACGAGTAATCGTTTATcc               >  1:1443000/1‑150 (MQ=255)
                                                                                                                                    gTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGTA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATcc               >  1:3904/1‑150 (MQ=255)
                                                                                                                                    gTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATcc               <  2:789604/150‑1 (MQ=255)
                                                                                                                                    gTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACCGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATcc               <  1:231024/150‑1 (MQ=255)
                                                                                                                                    gTTACGCGTACCACCAGGATCGTATGCGGCGCCGTGCTGCGGGTGAAGCTCCGGCTGCAGCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATcc               <  1:1247545/150‑1 (MQ=255)
                                                                                                                                      tACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGAATGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGc             <  2:1279877/150‑1 (MQ=255)
                                                                                                                                       aCGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCa            <  1:460287/150‑1 (MQ=255)
                                                                                                                                       aCGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCa            <  2:94148/150‑1 (MQ=255)
                                                                                                                                       aCGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCa            <  1:1102560/150‑1 (MQ=255)
                                                                                                                                        cgcgTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCaa           >  1:1307094/1‑150 (MQ=255)
                                                                                                                                              ccaccaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCGTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATTACGAGTAATCGTTAATCCGCAAATAACg     >  1:1399964/1‑150 (MQ=255)
                                                                                                                                              ccaccaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACg     <  2:804087/150‑1 (MQ=255)
                                                                                                                                                accaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACg     >  2:1652357/1‑150 (MQ=255)
                                                                                                                                                 ccaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTaa  >  1:1600854/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCA‑‑GGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAA  >  NC_000913/4189302‑4189596

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