Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 632,077 T→C 7.7% intergenic (+78/‑105) 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,0770TC 7.7% 546.4 / 26.5 170intergenic (+78/‑105)cstA/ybdDcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein
Reads supporting (aligned to +/- strand):  ref base T (74/81);  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 = 9.37e-01

Only 100 of 170 total aligned reads displayed.

AAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAG  >  NC_000913/631929‑632223
                                                                                                                                                     |                                                                                                                                                  
aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTg                                                                                                                                                   <  2:774865/150‑1 (MQ=255)
 agaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTgg                                                                                                                                                  >  1:2185697/1‑150 (MQ=255)
   aaaCGCCGTATGAACCCATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAACGGCGCACACTACAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTgggg                                                                                                                                                <  2:2474267/150‑1 (MQ=255)
   aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTgggg                                                                                                                                                >  1:2053766/1‑150 (MQ=255)
     aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                                                                                                                              >  1:452813/1‑150 (MQ=255)
     aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGAT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                                                                                                                                              >  1:452858/1‑150 (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                                                                                                                                           <  2:2139229/150‑1 (MQ=255)
         cGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg                                                                                                                                          <  2:1604413/150‑1 (MQ=255)
           tctGAACCAATGCCGGACAATGTCGAGGAGATCGTGGCGCAGGCAAACGGCGCCCACTAACGTCAGAGTGA‑GGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCGCCCCCCCGGGGAGAgggtta                                                                                                                                       <  2:2370249/148‑4 (MQ=255)
           tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                                                                                                                                        >  1:102960/1‑150 (MQ=255)
               aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                                                                                                    >  1:1796956/1‑150 (MQ=255)
               aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                                                                                                    >  1:1796957/1‑150 (MQ=255)
               aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                                                                                                    <  2:631858/150‑1 (MQ=255)
               aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGa                                                                                                                                    <  2:1839602/150‑1 (MQ=255)
                   aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa                                                                                                                                >  1:1686585/1‑150 (MQ=255)
                        cGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGcgcg                                                                                                                           <  2:72462/150‑1 (MQ=255)
                              aTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGACTGTTGGCTTAGTGCCCAGGGGT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATa                                                                                                                     <  2:1448828/150‑1 (MQ=255)
                               tGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg                                                                                                                    <  1:2573110/150‑1 (MQ=255)
                                 tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCa                                                                                                                  <  1:1925282/150‑1 (MQ=255)
                                 tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCCAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCa                                                                                                                  <  1:950042/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)
                                       agaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTa                                                                                                            <  2:494695/150‑1 (MQ=255)
                                          tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGA‑GGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGcc                                                                                                        <  1:477718/150‑1 (MQ=255)
                                            gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCg                                                                                                       >  2:293629/1‑150 (MQ=255)
                                            gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCg                                                                                                       <  1:1278788/150‑1 (MQ=255)
                                            gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGGGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGCCTGACCGAGCGCGGTGAGAGCATTTGTAGGCCg                                                                                                       >  2:1321545/1‑150 (MQ=255)
                                               gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGAt                                                                                                    >  1:109564/1‑150 (MQ=255)
                                               gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGAt                                                                                                    >  2:1377887/1‑150 (MQ=255)
                                                 gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATaa                                                                                                  <  2:741125/150‑1 (MQ=255)
                                                 gcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATaa                                                                                                  <  2:741105/150‑1 (MQ=255)
                                                  cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAg                                                                                                 <  1:346165/150‑1 (MQ=255)
                                                  cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAg                                                                                                 <  1:1079330/150‑1 (MQ=255)
                                                    ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                               >  2:871346/1‑150 (MQ=255)
                                                    ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                               >  1:1286313/1‑150 (MQ=255)
                                                          aGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAc                                                                                         <  1:493527/150‑1 (MQ=255)
                                                              gcACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                     <  2:2205814/150‑1 (MQ=255)
                                                                  acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                 <  1:2400883/150‑1 (MQ=255)
                                                                  acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                 >  2:618293/1‑150 (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)
                                                                      aaGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                             <  2:1835079/150‑1 (MQ=39)
                                                                         tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                          >  1:2157385/1‑150 (MQ=255)
                                                                         tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                          >  1:800562/1‑150 (MQ=255)
                                                                         tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                          >  1:2189333/1‑150 (MQ=255)
                                                                         tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAct                                                                          >  1:2252582/1‑150 (MQ=255)
                                                                             agTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCttt                                                                      >  1:923753/1‑150 (MQ=255)
                                                                             agTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCttt                                                                      <  2:2109339/150‑1 (MQ=255)
                                                                              gTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTACCCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTc                                                                     <  1:1680966/150‑1 (MQ=255)
                                                                                gAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAg                                                                   <  2:1891775/150‑1 (MQ=255)
                                                                                   ggggggATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCaa                                                                >  2:659431/1‑150 (MQ=255)
                                                                                     ggggATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACTCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACa                                                              >  2:1030629/1‑150 (MQ=255)
                                                                                     ggggATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACa                                                              >  2:533358/1‑150 (MQ=255)
                                                                                        gATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATgg                                                           >  2:319691/1‑150 (MQ=255)
                                                                                            ttGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAg                                                       >  2:1882058/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:1803503/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)
                                                                                                             aGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGttt                                      <  1:819927/150‑1 (MQ=255)
                                                                                                               tGCCCAGGGGT‑CCCTCTCACCCTAACCCTCTCCCCGGGGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGa                                    <  2:1976125/150‑1 (MQ=255)
                                                                                                                gCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                                   <  2:415100/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)
                                                                                                                 cccAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGAGGTTTGAtt                                  >  2:1529892/1‑150 (MQ=255)
                                                                                                                   cAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCa                                >  1:342073/1‑150 (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:2371961/1‑150 (MQ=255)
                                                                                                                              tctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAg                      <  1:544750/150‑1 (MQ=255)
                                                                                                                              tctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAg                      <  1:1451170/150‑1 (MQ=255)
                                                                                                                               ctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGc                     <  2:1427758/150‑1 (MQ=255)
                                                                                                                               ctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGc                     <  1:551960/150‑1 (MQ=255)
                                                                                                                                  aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg                  >  2:409108/1‑150 (MQ=255)
                                                                                                                                  aCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCgg                  >  2:1557500/1‑150 (MQ=255)
                                                                                                                                    ccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaa                <  2:943616/150‑1 (MQ=255)
                                                                                                                                     agtccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaaa               >  2:1080191/5‑150 (MQ=39)
                                                                                                                                     agtccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaaa               >  2:1080178/5‑150 (MQ=39)
                                                                                                                                       tccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAatat             >  2:1132640/3‑150 (MQ=38)
                                                                                                                                       tccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAat             >  2:2442395/3‑150 (MQ=39)
                                                                                                                                        aCCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAata            <  1:947644/150‑1 (MQ=255)
                                                                                                                                          ccTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATAtt          <  1:2524278/150‑1 (MQ=255)
                                                                                                                                                  cGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAg  <  1:1246299/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTT‑CCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAG  >  NC_000913/631929‑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: ‑