Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,677,833 A→G 5.9% I173T (ATC→ACC) ‡ yhjG ← AsmA family protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,677,8330AG 5.9% 462.8 / 13.4 152I173T (ATC→ACC) ‡yhjGAsmA family protein YhjG
Reads supporting (aligned to +/- strand):  ref base A (72/71);  new base G (6/3);  total (78/74)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 156 total aligned reads displayed.

CCTGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCG  >  NC_000913/3677684‑3677982
                                                                                                                                                     |                                                                                                                                                     
ccTGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGa                                                                                                                                                       >  1:1367790/1‑150 (MQ=255)
  tGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCAAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAAACTCCAGATACGCTTTGCTTACTTTGTCATCAATGGCGaac                                                                                                                                                     >  1:2111930/1‑148 (MQ=255)
    cgcCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGGTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCg                                                                                                                                                   >  1:2079376/1‑150 (MQ=255)
     gcCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGc                                                                                                                                                  <  2:670511/150‑1 (MQ=255)
      ccTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGcc                                                                                                                                                 >  1:1823757/1‑150 (MQ=255)
         tCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCtt                                                                                                                                              <  2:471644/150‑1 (MQ=255)
         tCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCtt                                                                                                                                              <  1:1539529/150‑1 (MQ=255)
            ggCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGCTCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGAt                                                                                                                                           >  1:2001856/1‑150 (MQ=255)
            ggCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGAt                                                                                                                                           >  2:116360/1‑150 (MQ=255)
            ggCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAAGGGCGGATCGCCATTGAt                                                                                                                                           <  1:1161031/150‑1 (MQ=255)
             gCCTAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACGTTTCGATCCAGTAACTTCGGTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATc                                                                                                                                          <  1:1695344/150‑1 (MQ=255)
             gCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCATTAACTTGGCTGAACGGAAGCGGCTTGCCTAACGGAGCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATc                                                                                                                                          <  2:1880978/150‑1 (MQ=255)
             gCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATc                                                                                                                                          <  1:2403/150‑1 (MQ=255)
              ccAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCg                                                                                                                                         >  1:167558/1‑150 (MQ=255)
              ccAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCg                                                                                                                                         >  1:1120677/1‑150 (MQ=255)
              ccAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGAACAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCg                                                                                                                                         >  2:1986639/1‑150 (MQ=255)
               cAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGa                                                                                                                                        >  1:183850/1‑150 (MQ=255)
               cAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAAGGGCGGATCGCCATTGAtgac                                                                                                                                        <  2:2423426/150‑4 (MQ=255)
                 aaaaCGTAATCGCCCACCTTTTCTTTCTCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaa                                                                                                                                      <  2:1602875/150‑1 (MQ=255)
                  aaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaaa                                                                                                                                     <  2:345304/150‑1 (MQ=255)
                  aaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaaa                                                                                                                                     <  2:388376/150‑1 (MQ=255)
                  aaaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGaaaa                                                                                                                                     <  1:1067216/150‑1 (MQ=255)
                   aaCGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAg                                                                                                                                    <  1:223110/150‑1 (MQ=255)
                      ttAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGTACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAat                                                                                                                                 <  2:1285990/149‑1 (MQ=255)
                      gTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGTATCAACAAATTTCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAat                                                                                                                                 <  2:1824254/150‑1 (MQ=255)
                      gTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGAACCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAat                                                                                                                                 >  1:1668004/1‑150 (MQ=255)
                       tAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAata                                                                                                                                <  1:2398591/150‑1 (MQ=255)
                        aaTCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAatat                                                                                                                               >  1:1915190/1‑150 (MQ=255)
                          tCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAAtatta                                                                                                                             <  1:796910/150‑1 (MQ=255)
                          tCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAAtatta                                                                                                                             <  2:2226181/150‑1 (MQ=255)
                               cACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAg                                                                                                                        >  2:306420/1‑150 (MQ=255)
                               cACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATACTATCCAg                                                                                                                        >  2:2092994/1‑150 (MQ=255)
                                 ccTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGcc                                                                                                                      <  2:1584819/150‑1 (MQ=255)
                                  cTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCg                                                                                                                     >  1:756940/1‑150 (MQ=255)
                                  cTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCg                                                                                                                     >  1:2436514/1‑150 (MQ=255)
                                   ttttCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGa                                                                                                                    <  2:1642144/150‑1 (MQ=255)
                                    tttCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGaa                                                                                                                   <  2:853072/150‑1 (MQ=255)
                                      tCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAAc                                                                                                                 <  1:34035/150‑1 (MQ=255)
                                        tttATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGa                                                                                                               >  2:2311273/1‑150 (MQ=255)
                                           aTCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAACTCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCa                                                                                                            <  2:1561678/150‑1 (MQ=255)
                                              cGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGc                                                                                                         >  2:35614/1‑150 (MQ=255)
                                              cGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGc                                                                                                         <  1:2000597/150‑1 (MQ=255)
                                              cGCTTTAAATTTAGATCCAGAAACTACGCTAAACGTCAGAAGATTGACTAAGAGATAAACAAAAATCTCCAGATAAGAATAGCTAAATTTCTCATCAATGGACATCCGCCCTTGAACGAAAAGAATATTATCCAGCCGAAACGCCCATGc                                                                                                         >  2:1160715/1‑150 (MQ=255)
                                               gCTTTACCTTTCGCTCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCTAAAAGAATATTATCCAGCCGAAACGACCATGcc                                                                                                        <  2:682868/150‑1 (MQ=255)
                                               gCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGcc                                                                                                        <  2:1395063/150‑1 (MQ=255)
                                                       tttCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCtt                                                                                                <  2:923501/150‑1 (MQ=255)
                                                       tttCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCtt                                                                                                >  1:369930/1‑150 (MQ=255)
                                                          cGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGc                                                                                             <  2:266688/150‑1 (MQ=255)
                                                          cGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGc                                                                                             >  2:2457595/1‑150 (MQ=255)
                                                               cAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcg                                                                                        <  1:116360/150‑1 (MQ=255)
                                                                aGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcgc                                                                                       >  2:1031622/1‑150 (MQ=255)
                                                                 gTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTcgcg                                                                                      <  1:752117/150‑1 (MQ=255)
                                                                       tCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTg                                                                                <  2:1367790/150‑1 (MQ=255)
                                                                           tGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTat                                                                            >  1:1191194/1‑150 (MQ=255)
                                                                            gAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCACATCCGCTTTGCTAACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTatc                                                                           >  2:2168395/1‑150 (MQ=255)
                                                                              aCGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTatcat                                                                         >  1:484541/1‑150 (MQ=255)
                                                                               cGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTatcatc                                                                        <  2:421252/150‑1 (MQ=255)
                                                                                ggCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCATACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCg                                                                       >  2:159727/1‑150 (MQ=255)
                                                                                 gCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGGATTCGCGTCTTTGTTATCATCGt                                                                      >  1:1749607/1‑150 (MQ=255)
                                                                                 gCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGt                                                                      >  2:2005173/1‑150 (MQ=255)
                                                                                   aGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTg                                                                    <  1:2270272/150‑1 (MQ=255)
                                                                                        cTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAg                                                               <  2:753818/150‑1 (MQ=255)
                                                                                        cTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAg                                                               <  2:1353441/150‑1 (MQ=255)
                                                                                           gCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGAtt                                                            >  1:813661/1‑150 (MQ=255)
                                                                                           gCCTAAGGGAACAACAAAAATCTCCAGATCCGCATTGCTTACATTGTCATCAATGGCGATCCGCCCATGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTCTTAACATCGATGCCAAGAtt                                                            >  1:1280525/1‑150 (MQ=255)
                                                                                             cTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTaa                                                          >  1:1687008/1‑150 (MQ=255)
                                                                                               aaGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAAc                                                        <  1:1994692/150‑1 (MQ=255)
                                                                                                 cggATCAACAAAAATCTCAAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATACGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGt                                                      >  2:2119406/2‑150 (MQ=255)
                                                                                                    aTCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCa                                                   <  1:1201927/150‑1 (MQ=255)
                                                                                                     tCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAg                                                  <  2:645510/150‑1 (MQ=255)
                                                                                                        acaaAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTa                                               >  1:2057089/1‑150 (MQ=255)
                                                                                                                ctcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCtt                                       >  2:1975204/1‑150 (MQ=255)
                                                                                                                 tcCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCttt                                      >  2:2252425/1‑150 (MQ=255)
                                                                                                                   cAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTc                                    >  2:2080630/1‑150 (MQ=255)
                                                                                                                      aTCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCaga                                 <  1:306420/150‑1 (MQ=255)
                                                                                                                              tGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGaa                         <  1:1500715/150‑1 (MQ=255)
                                                                                                                                cTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATc                       >  2:1203305/1‑150 (MQ=255)
                                                                                                                                 ttACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCa                      <  2:1864788/150‑1 (MQ=255)
                                                                                                                                  tACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAg                     >  2:327802/1‑150 (MQ=255)
                                                                                                                                  tACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAg                     >  2:2337129/1‑150 (MQ=255)
                                                                                                                                   aCTTTCTCATCAATGGCGATCCGCCCTTGATCGAAAAGAAAAATATCCACCCGAAAAGAACATACAGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAgg                    >  2:1652624/1‑150 (MQ=255)
                                                                                                                                    cttttcgatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAACCGACCATGCCGCCGGCTTTGGATTAGCGTCTTTGTTATCAGCGTGGGCATGCTTAAACGGCCAGTTATTGTTCTTTTCGGCGAGGCGAATCAGgc                   >  2:849176/8‑150 (MQ=255)
                                                                                                                                     tttGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGgcg                  <  1:2199028/150‑1 (MQ=255)
                                                                                                                                        tcgatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGgcgcgc               >  1:308547/4‑150 (MQ=255)
                                                                                                                                          gatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt             >  1:189612/2‑150 (MQ=255)
                                                                                                                                          gatcaAGGGCGGATCGCCATTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt             >  1:552835/2‑150 (MQ=255)
                                                                                                                                          catcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt             >  1:632929/1‑150 (MQ=255)
                                                                                                                                          catcaATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGt             >  1:2339728/1‑150 (MQ=255)
                                                                                                                                              aaTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTc                                    <  2:1049137/123‑1 (MQ=255)
                                                                                                                                              aaTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTc                                    >  1:1049137/1‑123 (MQ=255)
                                                                                                                                              aaTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCggg         >  1:1874355/1‑150 (MQ=255)
                                                                                                                                              aaTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCggg         <  2:1324070/150‑1 (MQ=255)
                                                                                                                                               aTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGc        <  2:1823757/150‑1 (MQ=255)
                                                                                                                                               aTGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGAATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGc        >  2:1792460/1‑150 (MQ=255)
                                                                                                                                               aTGGCGATCCGCCCGTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGGTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGc        <  2:1315078/150‑1 (MQ=255)
                                                                                                                                                  gCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCttt     >  1:1793798/1‑150 (MQ=255)
                                                                                                                                                  gCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCttt     >  1:277689/1‑150 (MQ=255)
                                                                                                                                                    gATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTc   >  1:1111007/1‑150 (MQ=255)
                                                                                                                                                    gATCCGCCCTCGAGCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCTTTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTc   <  1:999329/150‑1 (MQ=255)
                                                                                                                                                     aTCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCg  >  1:1595966/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTGCGCCTTCAGGCCAAAAACGTAATCGCCCACCTTTTCTTTATCCGCTTTACCTTTCGATCCAGTAACTTCGCTGAACGGCAGCGGCTTGCCTAAGGGATCAACAAAAATCTCCAGATCCGCTTTGCTTACTTTGTCATCAATGGCGATCCGCCCTTGATCGAAAAGAATATTATCCAGCCGAAACGACCATGCCGACGGCTTTGCATTCGCGTCTTTGTTATCATCGTTGGCAAGATTAAACGTCCAGTTATTGTTCTTTTCAGAGAGGCGAATCAGGCGCGCGTCGGGCTTTTCG  >  NC_000913/3677684‑3677982

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