Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 632,076 G→C 7.6% intergenic (+77/‑106) 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,0760GC 7.6% 615.1 / 35.8 170intergenic (+77/‑106)cstA/ybdDcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein
Reads supporting (aligned to +/- strand):  ref base G (75/82);  new base C (5/8);  total (80/90)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 172 total aligned reads displayed.

GAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAG  >  NC_000913/631927‑632223
                                                                                                                                                      |                                                                                                                                                   
gaaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCgg                                                                                                                                                     <  1:1701103/150‑1 (MQ=255)
  aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTg                                                                                                                                                   <  2:774865/150‑1 (MQ=255)
     aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTgggg                                                                                                                                                <  2:1581657/150‑1 (MQ=255)
       aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                                                                                                                              >  1:490038/1‑150 (MQ=255)
       aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                                                                                                                              >  1:452813/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:2139229/150‑1 (MQ=255)
          ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCGCCCCACCGGGGAGAgg                                                                                                                                           <  1:1786925/150‑1 (MQ=255)
             tctGAACCAATGCCGGACAATGTCGAGGAGATCGTGGCGCAGGCAAACGGCGCCCACTAACGTCAGAGTGA‑GGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCGCCCCCCCGGGGAGAgggtta                                                                                                                                       <  2:2370249/148‑4 (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                                                                                                                                    <  2:631858/150‑1 (MQ=255)
                 aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                                                                                                    >  1:1796957/1‑150 (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)
                          cGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGcgcg                                                                                                                           <  2:1383788/150‑1 (MQ=255)
                               aaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt                                                                                                                      <  2:1545882/150‑1 (MQ=255)
                                aTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGACTGTTGGCTTAGTGCCCAGGGGT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATa                                                                                                                     <  2:1448828/150‑1 (MQ=255)
                                     gaggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCAtt                                                                                                                >  1:370779/1‑150 (MQ=255)
                                      aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCAttt                                                                                                               >  1:1186939/1‑150 (MQ=255)
                                       ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTg                                                                                                              >  2:544750/1‑150 (MQ=255)
                                        gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGt                                                                                                             >  2:2225514/1‑150 (MQ=255)
                                        gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGt                                                                                                             >  2:1453238/1‑150 (MQ=255)
                                              gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCg                                                                                                       >  2:293629/1‑150 (MQ=255)
                                              gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCg                                                                                                       <  1:1278788/150‑1 (MQ=255)
                                                 gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGAt                                                                                                    >  2:1377882/1‑150 (MQ=255)
                                                 gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGAt                                                                                                    >  1:109564/1‑150 (MQ=255)
                                                   gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATaa                                                                                                  <  2:741125/150‑1 (MQ=255)
                                                   gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATaa                                                                                                  <  2:741105/150‑1 (MQ=255)
                                                    cAGGCAAAAGGCGCACACTAACGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAg                                                                                                 <  2:1836030/150‑1 (MQ=255)
                                                    cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAg                                                                                                 <  1:1079330/150‑1 (MQ=255)
                                                      ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                               >  1:1286313/1‑150 (MQ=255)
                                                      ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                               >  2:871346/1‑150 (MQ=255)
                                                      ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                               >  1:992820/1‑150 (MQ=255)
                                                           aaGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCa                                                                                          >  2:1988536/1‑150 (MQ=255)
                                                            aagcgcACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAc                                                                                         >  2:106811/3‑150 (MQ=255)
                                                            aGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAc                                                                                         <  1:493527/150‑1 (MQ=255)
                                                              gcgcACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGCCTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATCAGGCGTTCAcgc                                                                                       >  1:600805/1‑150 (MQ=255)
                                                                gcACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                     <  2:2205814/150‑1 (MQ=255)
                                                                    acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                 <  1:2400879/150‑1 (MQ=255)
                                                                    acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                 <  1:2400883/150‑1 (MQ=255)
                                                                      tAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                               >  2:169437/1‑150 (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:1835273/150‑1 (MQ=39)
                                                                           tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                          >  1:2189305/1‑150 (MQ=255)
                                                                           tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                          >  1:800562/1‑150 (MQ=255)
                                                                           tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                          >  1:2252582/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)
                                                                                  gAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAg                                                                   <  2:1891775/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)
                                                                                          gATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATgg                                                           >  2:319691/1‑150 (MQ=255)
                                                                                                ggCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCg                                                     >  2:669038/1‑150 (MQ=255)
                                                                                                   gAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAg                                                  <  1:2352132/150‑1 (MQ=255)
                                                                                                     aTGTTGGCTTAGTGCCCAGGGTTCCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAgg                                                 <  1:1362991/150‑1 (MQ=255)
                                                                                                     aTGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGc                                                >  1:184502/1‑150 (MQ=255)
                                                                                                          ggCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGa                                           <  2:271869/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)
                                                                                                                 tGCCCAGGGGT‑CCCTCTCACCCTAACCCTCTCCCCGGGGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGa                                    <  2:1976125/150‑1 (MQ=255)
                                                                                                                  gCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                                   <  2:1649632/150‑1 (MQ=255)
                                                                                                                  gCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                                   <  2:415100/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)
                                                                                                                      aGGGTG‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCAc                               <  1:1988536/150‑1 (MQ=255)
                                                                                                                       gggTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACt                              <  1:1398591/150‑1 (MQ=255)
                                                                                                                        ggTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTg                             <  1:323609/150‑1 (MQ=255)
                                                                                                                          tt‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGc                           >  2:455748/1‑150 (MQ=255)
                                                                                                                             cccTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCaa                         >  1:364163/1‑150 (MQ=255)
                                                                                                                                tctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAg                      <  1:544750/150‑1 (MQ=255)
                                                                                                                                tctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAg                      <  1:1451170/150‑1 (MQ=255)
                                                                                                                                 ctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGc                     <  1:2169496/150‑1 (MQ=255)
                                                                                                                                   cACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCg                   >  2:1180205/1‑150 (MQ=255)
                                                                                                                                    aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg                  <  2:370779/150‑1 (MQ=255)
                                                                                                                                    aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg                  <  2:184502/150‑1 (MQ=255)
                                                                                                                                    aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACAGGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg                  >  2:1101567/1‑150 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaa                <  2:943616/150‑1 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaa                >  1:1731512/1‑150 (MQ=255)
                                                                                                                                       agtccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaaa               >  2:1080178/5‑150 (MQ=39)
                                                                                                                                         tccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAat             >  2:2442395/3‑150 (MQ=39)
                                                                                                                                          aCCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAata            <  1:947644/150‑1 (MQ=255)
                                                                                                                                            ccTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATAtt          <  1:2524278/150‑1 (MQ=255)
                                                                                                                                                 ccccGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGt     <  2:1057303/150‑1 (MQ=255)
                                                                                                                                                    cGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAg  <  1:1246321/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAG  >  NC_000913/631927‑632223

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: ‑