Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,613,961 G→T 7.9% E99* (GAA→TAA) ‡ nikA → nickel/heme ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,613,9610GT 7.9% 526.8 / 25.5 145E99* (GAA→TAA) ‡nikAnickel/heme ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (74/59);  new base T (9/3);  total (83/62)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

TCGCCCAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTA  >  NC_000913/3613815‑3614109
                                                                                                                                                  |                                                                                                                                                    
tCGCCCAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAAc                                                                                                                                                   >  1:2665832/1‑150 (MQ=255)
   cccAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGt                                                                                                                                                <  2:2562120/150‑1 (MQ=255)
   cccAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGt                                                                                                                                                <  1:2597458/150‑1 (MQ=255)
      agagCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTTCACCGttgg                                                                                                                                             <  1:1441478/150‑3 (MQ=255)
      agagCAGGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATGCCGTGGCTGGCAAACAGCTGGACTCTTTCAGAAGAGGGTAAAACCTGGACCTTCACCCTGCGTGTTGACGTGAAATTCTCCAACGGTTCACCGttgg                                                                                                                                             <  1:169691/150‑3 (MQ=255)
          cATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTGCGATGc                                                                                                                                         >  1:742812/1‑150 (MQ=255)
                ttATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAggc                                                                                                                                   >  1:870947/1‑150 (MQ=255)
                 tATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAggcg                                                                                                                                  >  1:671251/1‑150 (MQ=255)
                 tATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAggcg                                                                                                                                  >  1:650171/1‑150 (MQ=255)
                 tATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAggcg                                                                                                                                  >  1:1371272/1‑150 (MQ=255)
                         aTTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGa                                                                                                                          >  2:2498160/1‑150 (MQ=255)
                          ttGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGaa                                                                                                                         <  1:2005246/150‑1 (MQ=255)
                            ggTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGaaaa                                                                                                                       <  1:1147944/150‑1 (MQ=255)
                              tGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACt                                                                                                                     >  2:459401/1‑150 (MQ=255)
                              tGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACt                                                                                                                     >  2:1178871/1‑150 (MQ=255)
                                    atCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACATGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGACAACTTCcgcg                                                                                                               >  1:2264849/1‑150 (MQ=255)
                                      cAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCa                                                                                                             <  2:1804906/150‑1 (MQ=255)
                                         gCAGACGTTTCGGTGATCCCGTGGCTGGCAAAATGCTGGACTCATTCAGAAGTTGGTAAAACCTGGACCTTCACCCTGCTTGATGACGTTAAATTCTCCAACGGTTAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                          <  1:2074751/150‑1 (MQ=255)
                                         gCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                          <  2:1463435/150‑1 (MQ=255)
                                         gCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                          >  2:778708/1‑150 (MQ=255)
                                         gCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                          <  2:1009524/150‑1 (MQ=255)
                                         gCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                          >  2:2737383/1‑150 (MQ=255)
                                            gACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTc                                                                                                       <  2:1034951/150‑1 (MQ=255)
                                            gACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTc                                                                                                       <  2:594430/150‑1 (MQ=255)
                                                ggtCGGTGATCCCGTGGCTGGCAAACAGCTGGACTCCTTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATa                                                                                                   <  2:1697136/148‑1 (MQ=255)
                                                gTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATa                                                                                                   <  2:1593040/150‑1 (MQ=255)
                                                gTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATa                                                                                                   <  2:1697140/150‑1 (MQ=255)
                                                   cGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGGGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAAcc                                                                                                <  1:2501122/150‑1 (MQ=255)
                                                   cGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATACAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAAcc                                                                                                <  1:863538/150‑1 (MQ=255)
                                                    ggTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCg                                                                                               >  2:1555444/1‑150 (MQ=255)
                                                      tGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTc                                                                                             >  1:692933/1‑150 (MQ=255)
                                                         tCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAAc                                                                                          <  2:1950146/150‑1 (MQ=255)
                                                              tggctggcAAAAAGCTGGACTCACTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCAc                                                                                     >  2:2714531/1‑150 (MQ=255)
                                                              tggctggcAAAAAGCTGGACTCACTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCAc                                                                                     >  2:264227/1‑150 (MQ=255)
                                                                      aaaaaGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCctggct                                                                             >  2:413036/1‑150 (MQ=255)
                                                                          aGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAg                                                                         <  1:561931/150‑1 (MQ=255)
                                                                             tGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTg                                                                      <  1:2167303/150‑1 (MQ=255)
                                                                              ggACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTgg                                                                     >  1:2457247/1‑150 (MQ=255)
                                                                               gACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGc                                                                    >  2:731766/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                 >  2:1384988/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                 >  2:206532/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                 >  2:2565672/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                 >  1:1862675/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                 >  1:466713/1‑150 (MQ=255)
                                                                                     ttCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                 <  2:1162337/147‑1 (MQ=255)
                                                                                     ttCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCa                                                              >  2:2445401/1‑150 (MQ=255)
                                                                                     ttCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCa                                                              >  2:1028351/1‑150 (MQ=255)
                                                                                       cAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGa                                                            <  1:2302877/150‑1 (MQ=255)
                                                                                       cAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGa                                                            <  2:1138562/150‑1 (MQ=255)
                                                                                          aagaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAttg                                                         >  2:422571/1‑150 (MQ=255)
                                                                                          aagaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAgtg                                                         >  2:1565832/1‑150 (MQ=255)
                                                                                            gaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAttgtt                                                       <  1:2270360/150‑1 (MQ=255)
                                                                                            gaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAttgtt                                                       <  2:2039423/150‑1 (MQ=255)
                                                                                            gaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAttgtt                                                       <  2:914138/150‑1 (MQ=255)
                                                                                               ggTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGAt                                                    <  1:602277/150‑1 (MQ=255)
                                                                                               ggTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGAt                                                    <  1:1094520/150‑1 (MQ=255)
                                                                                               ggTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGAt                                                    <  1:2348056/150‑1 (MQ=255)
                                                                                                 tAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGt                                                  >  2:2382063/1‑150 (MQ=255)
                                                                                                   aaaCCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTa                                                >  1:2475725/1‑150 (MQ=255)
                                                                                                       cTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGc                                            >  2:558067/1‑150 (MQ=255)
                                                                                                       cTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGc                                            >  1:2564918/1‑150 (MQ=255)
                                                                                                        tGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCa                                           >  1:2105360/1‑150 (MQ=255)
                                                                                                              ttCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGt                                     <  2:1373237/150‑1 (MQ=255)
                                                                                                              ttCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGt                                     <  2:1373229/150‑1 (MQ=255)
                                                                                                               tCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTa                                    <  2:315024/150‑1 (MQ=255)
                                                                                                                  cccTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTaaaa                                 <  2:1253095/150‑1 (MQ=255)
                                                                                                                    cTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACa                               <  2:573669/150‑1 (MQ=255)
                                                                                                                    cTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACa                               <  1:1555444/150‑1 (MQ=255)
                                                                                                                      gCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACaga                             <  1:886884/150‑1 (MQ=255)
                                                                                                                      gCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACaga                             <  2:870947/150‑1 (MQ=255)
                                                                                                                      gCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACaga                             <  1:744715/150‑1 (MQ=255)
                                                                                                                      gCGTGATGACGTGAAATTCTCCAACGGTGAACCGGTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACaga                             <  2:2347381/150‑1 (MQ=255)
                                                                                                                           atgaCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGc                        >  1:1665851/1‑150 (MQ=255)
                                                                                                                              aCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCaaa                     >  1:21594/1‑150 (MQ=255)
                                                                                                                              aCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCaaa                     >  1:2601037/1‑150 (MQ=255)
                                                                                                                              aCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCaaa                     >  2:639332/1‑150 (MQ=255)
                                                                                                                              aCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCaaa                     >  1:1145608/1‑150 (MQ=255)
                                                                                                                                      ttCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGGAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTg             >  2:1836728/1‑150 (MQ=255)
                                                                                                                                      ttCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTg             >  2:950192/1‑150 (MQ=255)
                                                                                                                                           cAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAg        >  1:886614/1‑150 (MQ=255)
                                                                                                                                           cAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAg        >  1:30118/1‑150 (MQ=255)
                                                                                                                                           cAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAg        >  1:2203823/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAAGGTCACGCCTGGCTGGAGCTGGCAAACCAGAGTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc       >  1:311379/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc       >  1:879322/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc       >  1:1232154/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc       >  1:2679058/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc       >  1:1985577/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc       >  1:2187593/1‑150 (MQ=255)
                                                                                                                                            aaCGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc       <  1:110206/150‑1 (MQ=255)
                                                                                                                                             aCGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgcg      <  1:1095892/150‑1 (MQ=255)
                                                                                                                                             aCGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgcg      <  1:1274657/150‑1 (MQ=255)
                                                                                                                                               ggTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGcc    >  2:1467467/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCct   >  1:1981555/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCct   >  1:1800648/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCct   >  2:303898/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCct   >  2:1026146/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAAATACCCTGAAAAGCGCct   >  2:1499278/1‑150 (MQ=255)
                                                                                                                                                 tGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCcta  >  2:523130/1‑150 (MQ=255)
                                                                                                                                                 tGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCcta  >  1:2556043/1‑150 (MQ=255)
                                                                                                                                                 tGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCcta  >  1:65675/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGCCCAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTA  >  NC_000913/3613815‑3614109

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