Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 632,070 T→G 7.6% intergenic (+71/‑112) cstA → / → ybdD carbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913632,0700TG 7.6% 552.8 / 33.0 172intergenic (+71/‑112)cstA/ybdDcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein
Reads supporting (aligned to +/- strand):  ref base T (76/83);  new base G (5/8);  total (81/91)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

Only 100 of 172 total aligned reads displayed.

GCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATT  >  NC_000913/631925‑632215
                                                                                                                                                  |                                                                                                                                                 
gCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTcccc                                                                                                                                               <  2:391167/150‑1 (MQ=255)
gCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTcccc                                                                                                                                               <  2:231575/150‑1 (MQ=255)
  gaaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCgg                                                                                                                                             <  1:1701103/150‑1 (MQ=255)
    aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTg                                                                                                                                           <  2:774865/150‑1 (MQ=255)
     agaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTgg                                                                                                                                          >  1:2185697/1‑150 (MQ=255)
       aaaCGCCGTATGAACCCATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAACGGCGCACACTACAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTgggg                                                                                                                                        <  2:2474267/150‑1 (MQ=255)
         aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                                                                                                                      >  1:452813/1‑150 (MQ=255)
         aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                                                                                                                      >  1:490038/1‑150 (MQ=255)
         aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGACTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                                                                                                                      >  1:452867/1‑150 (MQ=255)
          cGCCGTATTCACCATTGCGGGAAAAGGCCTTGGAGATCGGGGCGCAGGCAAAAGGCGCACCATAATGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGa                                                                                                                                     <  2:1772583/150‑1 (MQ=255)
           gCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGTCGAATGTTTGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                                                                                                                                    <  2:2469008/150‑1 (MQ=255)
            ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                                                                                                                                   <  1:1116499/150‑1 (MQ=255)
            ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                                                                                                                                   >  1:2475528/1‑150 (MQ=255)
            ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCGCCCCACCGGGGAGAgg                                                                                                                                   <  1:1059934/150‑1 (MQ=255)
            ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCGCCCCACCGGGGAGAgg                                                                                                                                   <  2:2139216/150‑1 (MQ=255)
            ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCGCCCCACCGGGGAGAgg                                                                                                                                   <  1:1786925/150‑1 (MQ=255)
               tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTCAAGTCAGAGTGAGGGGGGGATGGTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggggc                                                                                                                                >  2:1812548/1‑149 (MQ=255)
               tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                                                                                                                                >  1:102960/1‑150 (MQ=255)
                  gAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCGCCCCACCGGGGAGAgggttagg                                                                                                                             <  2:2125004/150‑6 (MQ=255)
                   aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                                                                                            >  1:1796957/1‑150 (MQ=255)
                   aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                                                                                            <  2:631858/150‑1 (MQ=255)
                   aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                                                                                            <  2:1609078/150‑1 (MQ=255)
                      cAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg                                                                                                                         >  1:916361/1‑150 (MQ=255)
                      cAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg                                                                                                                         >  2:480519/1‑150 (MQ=255)
                                 aaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt                                                                                                              <  2:1545882/150‑1 (MQ=255)
                                     tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCa                                                                                                          <  1:1925282/150‑1 (MQ=255)
                                     tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCCAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCa                                                                                                          <  1:950042/150‑1 (MQ=255)
                                        aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCAttt                                                                                                       >  1:1186939/1‑150 (MQ=255)
                                         ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTg                                                                                                      >  2:544750/1‑150 (MQ=255)
                                          gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGt                                                                                                     >  2:1453238/1‑150 (MQ=255)
                                          gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGt                                                                                                     >  2:2225514/1‑150 (MQ=255)
                                                   gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGAt                                                                                            >  1:109564/1‑150 (MQ=255)
                                                   gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGAt                                                                                            >  2:1377882/1‑150 (MQ=255)
                                                   gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGAt                                                                                            >  2:1377887/1‑150 (MQ=255)
                                                   gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGAt                                                                                            >  2:2415027/1‑150 (MQ=255)
                                                     gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATaa                                                                                          <  2:741105/150‑1 (MQ=255)
                                                     gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATaa                                                                                          <  2:741125/150‑1 (MQ=255)
                                                      cAGGCAAAAGGCGCACACTAACGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAg                                                                                         <  2:1836030/150‑1 (MQ=255)
                                                      cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAg                                                                                         <  1:1079330/150‑1 (MQ=255)
                                                       aGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAgg                                                                                        <  1:81564/150‑1 (MQ=255)
                                                        ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                       >  1:1286313/1‑150 (MQ=255)
                                                        ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                       >  1:1446253/1‑150 (MQ=255)
                                                        ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                       >  1:992820/1‑150 (MQ=255)
                                                        ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                       >  1:1286184/1‑150 (MQ=255)
                                                        ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                       >  2:871346/1‑150 (MQ=255)
                                                             aaGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCa                                                                                  >  2:1988536/1‑150 (MQ=255)
                                                                gcgcACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGCCTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATCAGGCGTTCAcgc                                                                               >  1:600805/1‑150 (MQ=255)
                                                                      acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                         <  1:2400879/150‑1 (MQ=255)
                                                                         aaaGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                      <  2:1186939/150‑1 (MQ=255)
                                                                          caGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTACCCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                     <  2:589530/149‑1 (MQ=255)
                                                                          aaGTCAGAGTGAGGGGGGGATGTTGGCGATTGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                     <  1:109496/150‑1 (MQ=38)
                                                                          aaGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                     <  2:1835079/150‑1 (MQ=39)
                                                                             tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                  >  1:2252582/1‑150 (MQ=255)
                                                                             tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                  >  1:416006/1‑150 (MQ=255)
                                                                             tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                  >  1:2189305/1‑150 (MQ=255)
                                                                                 agTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCttt                                                              >  1:923767/1‑150 (MQ=255)
                                                                                  gTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTACCCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTc                                                             <  1:1680966/150‑1 (MQ=255)
                                                                                   tGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCa                                                            <  1:2133688/150‑1 (MQ=255)
                                                                                       ggggggATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCaa                                                        >  2:659431/1‑150 (MQ=255)
                                                                                        gggggATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTGTCAGCAAc                                                       >  2:613419/1‑150 (MQ=255)
                                                                                         ggggATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACa                                                      >  2:533358/1‑150 (MQ=255)
                                                                                                ttGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAg                                               >  2:1882058/1‑150 (MQ=255)
                                                                                                  ggCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCg                                             >  2:669038/1‑150 (MQ=255)
                                                                                                  ggCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCg                                             <  2:1061439/150‑1 (MQ=255)
                                                                                                   gCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCgg                                            >  2:2336527/1‑150 (MQ=255)
                                                                                                    cGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGa                                           <  2:1621670/150‑1 (MQ=255)
                                                                                                       aTGTTGGCTTAGTGCCCAGGGTTCCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAgg                                         <  1:1362991/150‑1 (MQ=255)
                                                                                                          ttGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCaa                                     <  1:1765206/150‑1 (MQ=255)
                                                                                                            ggCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGa                                   <  2:272006/150‑1 (MQ=255)
                                                                                                              cTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATg                                 >  2:1585070/1‑150 (MQ=255)
                                                                                                              cTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATg                                 >  2:78701/1‑150 (MQ=255)
                                                                                                                 aGTGCCCAGTGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGttt                              <  1:819914/150‑1 (MQ=255)
                                                                                                                   tGCCCAGGGGT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGa                            <  1:313441/150‑1 (MQ=255)
                                                                                                                    gCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                           <  2:415100/150‑1 (MQ=255)
                                                                                                                    gCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                           <  2:1918549/150‑1 (MQ=255)
                                                                                                                    gCCCAGGGGT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                           <  2:1649807/150‑1 (MQ=255)
                                                                                                                     cccAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGCGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTACCGGAGGCCAAGATGTTTGAtt                          >  1:2382832/1‑150 (MQ=255)
                                                                                                                     cccAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAtt                          >  1:2191413/1‑150 (MQ=255)
                                                                                                                     cccAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGAGGTTTGAtt                          >  1:2168683/1‑150 (MQ=255)
                                                                                                                       cAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCa                        >  1:342073/1‑150 (MQ=255)
                                                                                                                       cAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTATGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCa                        >  1:341904/1‑150 (MQ=255)
                                                                                                                        aGGGTG‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCAc                       <  1:1988536/150‑1 (MQ=255)
                                                                                                                         gggTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACt                      >  1:518772/1‑150 (MQ=255)
                                                                                                                            tt‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGc                   >  2:455748/1‑150 (MQ=255)
                                                                                                                               cccTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCaa                 >  1:364163/1‑150 (MQ=255)
                                                                                                                               cccTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCaa                 <  2:296011/150‑1 (MQ=255)
                                                                                                                               cccTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCaa                 >  1:2371961/1‑150 (MQ=255)
                                                                                                                                  tctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAg              <  1:1451170/150‑1 (MQ=255)
                                                                                                                                  tctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAg              <  1:544750/150‑1 (MQ=255)
                                                                                                                                   ctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGc             <  2:1427758/150‑1 (MQ=255)
                                                                                                                                   ctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGc             <  1:2444218/150‑1 (MQ=255)
                                                                                                                                     cACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGGTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCACCATGGTTAGCGGAGGCCAAGATGGTTGATTCACGGGCAAAAGCCg           >  2:2154730/1‑150 (MQ=255)
                                                                                                                                      aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg          <  2:370779/150‑1 (MQ=255)
                                                                                                                                      aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg          <  2:184502/150‑1 (MQ=255)
                                                                                                                                      aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg          >  2:1557500/1‑150 (MQ=255)
                                                                                                                                      aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACAGGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg          >  2:1101567/1‑150 (MQ=255)
                                                                                                                                        ccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaa        <  2:943616/150‑1 (MQ=255)
                                                                                                                                         agtccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaaa       >  2:1080178/5‑150 (MQ=39)
                                                                                                                                           tccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAatat     >  2:1132640/3‑150 (MQ=38)
                                                                                                                                              ccTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATAtt  <  1:2524278/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATT  >  NC_000913/631925‑632215

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