Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,261,8420TG25.8% 114.4 / 10.2 66T12P (ACC→CCC) prsribose‑phosphate diphosphokinase
Reads supporting (aligned to +/- strand):  ref base T (9/40);  new base G (16/1);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 3.37e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCAT  >  NC_000913/1261601‑1262089
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          
aaGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTAcc                                                                                                                                                                                                                                                   <  1:1579/251‑1 (MQ=255)
aaGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGGTGGCAATATGTTTTGCTAGTTCCGGGGGGGGGTTAcc                                                                                                                                                                                                                                                   >  2:374947/1‑251 (MQ=255)
 aGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGGGGCGTTAccc                                                                                                                                                                                                                                                  >  1:418278/1‑250 (MQ=255)
  gTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGTCTTCGCCATCGCTAAAGCGACCTACAGCGGCTTCGCCGAGTGAAGTGTACAGGCGGTTGTCAATACGTTGTGCTAGTTCCGGGGGGGCGTTACCAg                                                                                                                                                                                                                                                 <  1:408721/251‑1 (MQ=255)
  gTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCAGCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAg                                                                                                                                                                                                                                                 <  2:528639/251‑1 (MQ=255)
        gATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCCATACGTTGTGCTAGTTTCGGGGTGGGGGTACCAGCaaaa                                                                                                                                                                                                                                            >  1:942/1‑250 (MQ=255)
             cAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTc                                                                                                                                                                                                                                      <  2:294971/251‑1 (MQ=255)
                         gACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTCACCAGCAAAAAGCTTCATATCAGGCACg                                                                                                                                                                                                                         <  1:2500/251‑1 (MQ=255)
                              gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGATAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATAAGTTGTGGTAGTTCCGGGGGGGGGTTACCAACAAAAAACTTCATATCAGGcccggga                                                                                                                                                                                                                      >  2:333872/1‑243 (MQ=255)
                              gCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATGCCATCAGGTTGTCGGTAGTAGGGGCACAAGTGGACTGGATGATGAAAAGATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGagaag                                                                                                                                                                                                                    <  2:49459/251‑1 (MQ=255)
                                  aaGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGTGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAAcc                                                                                                                                                                                                                <  1:445078/251‑1 (MQ=255)
                                        gacgCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGc                                                                                                                                                                                                          <  1:298356/251‑1 (MQ=255)
                                                gCATCAACCATAACGACTAATTCCATCAGTTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCCCCACCGCGTTCATTTTCATTAATTTGTACGCTGACTTCGCCATCGCGAAAGCGACCTACATCGGCGTCTCCGAGTGAAGTGTACAGGCGGTTGGCAATACTTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCGTCATATCAGGCACGATAAAAACCTCAGGCATGCGTcc                                                                                                                                                                                                  <  2:357127/251‑1 (MQ=255)
                                                     aaCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTACATtgg                                                                                                                                                                                             <  1:324464/251‑1 (MQ=255)
                                                      aCCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATtggt                                                                                                                                                                                            <  1:444110/251‑1 (MQ=255)
                                                           aaCGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCTCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGaaa                                                                                                                                                                                       <  1:293736/251‑1 (MQ=255)
                                                               aCTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGGTGTGCTAGTTCCGGGGTGGGGTTACCAGCAAAAAACTTCATATCCGGGACGGGAAGAAACTCCGGCATGCggccctttgggggaagaa                                                                                                                                                                                    >  1:416291/1‑233 (MQ=255)
                                                                  aaTTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGGCCATTGGTGGAAAGAATCt                                                                                                                                                                                 <  1:431693/250‑1 (MQ=255)
                                                                          cAGGTTGTAGTTAGTAGGGGCCAAAGTGGACTGGATGATGAAAATATCACCACCGCGTCCATTTTCATTAATTTGTACGCTGACTTAGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCGATACGTTGTGCTAGTTCCCGGGGGGGGGTTCCAGCAAAAAGCTTTCTATCAGGAAAGAGAAGAACCACCCGCAGGCGTCCACTGGAGGAAAGaaattaccgcaaac                                                                                                                                                                        >  2:486945/1‑239 (MQ=255)
                                                                            ggTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAACTATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTAGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGCGTGAAGTGTACAGGCGGTTGGGAATAAGTTGTGCTTGTTCCGGGGGGGGGTTACCCGCAAAAAGCTTTATATCCGGCACGAGAAGAAACTCCGGCATGCGTCCCTTGGTGGGAAGAAACTGCCGAAAACtg                                                                                                                                                                     >  1:508705/1‑251 (MQ=255)
                                                                               tGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCGCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACtgtg                                                                                                                                                                   <  1:474018/250‑1 (MQ=255)
                                                                                     tagtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGcc                                                                                                                                                             <  2:418278/250‑1 (MQ=255)
                                                                                     tagtagGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGcc                                                                                                                                                             <  2:437558/250‑1 (MQ=255)
                                                                                              cacaAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATc                                                                                                                                                   <  2:942/251‑1 (MQ=255)
                                                                                               acaAGTGGACTGGATGATGACAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTTCGGGGTGGCGGTACCAGCAAAAAGCTTCATATCCGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAAACTGCCGGAAAACGTGCGGGCCAGGCATGATcc                                                                                                                                                  >  1:164451/1‑251 (MQ=255)
                                                                                               acaAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATcc                                                                                                                                                  <  1:6402/251‑1 (MQ=255)
                                                                                                caAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCt                                                                                                                                                 <  2:143865/251‑1 (MQ=255)
                                                                                                    tGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCCGGGGGGCGTTACCAGCAAAAAGATTCATATACAGGACGGGAGGAACCTCCCGCACGCGTCCATTGTATGTAAGAAACTGCCGAAAACACTGCGGGCCAAGCaagaaccctttc                                                                                                                                             >  2:406944/1‑240 (MQ=255)
                                                                                                     ggACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTGACCAGCAAAAAGCTTCATATCAGTCACGAGAAGAACCTCAGGCATGCGTCCAGTGGTGGAAAGAATCTGCCGAAAACTGGGCGGGCCAGGCATGATCCTTTCCa                                                                                                                                            <  2:508705/251‑1 (MQ=255)
                                                                                                          ggATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                                                                                                                                        <  2:85820/250‑1 (MQ=255)
                                                                                                                 gaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAgg                                                                                                                                                                                                           <  1:348306/177‑1 (MQ=255)
                                                                                                                     atatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGGAATACGTTGTGCTAGTTCCGGGGGGGCGGTACCCGCAAAAAACTTCATATCAGGCACGGGAAGAAACTCAGGCATGCGTCCATTGGGGGAAAGAATATGCCGAAAACTGTGCGGGCCAGGCATGTTCCTTTTCACGCGGTGTAtttaac                                                                                                                             >  1:173085/1‑246 (MQ=255)
                                                                                                                      tatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCCGGGGGGCGGTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACACCAGGCAGGCGTCCCTTGGTGGAAAAAATCTGCCGAAAACTGTGCGGGGCCGGCACGATCCTTTCCAGGCGGGGTGTTaaaa                                                                                                                             >  2:164149/1‑248 (MQ=255)
                                                                                                                         caccacCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTTGGAATACGTTGTGCTAGTTCCGGGGGGGCGTTACCCGCAAAAAGCTTCATATCAGGCACGGGAAGAACCCCAGGCATGCGGCCATTGGTGGGAAGAATCTTCCGGAAAATGTGCGGGCCAGGGAAGATCCTTTCCAGGAGGGGGATTaaaaagggc                                                                                                                        >  1:318989/1‑245 (MQ=255)
                                                                                                                              ccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGc                                                                                                                  <  1:418007/251‑1 (MQ=255)
                                                                                                                                  gTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGATCAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGGCATCCGTTGTGCTAGTTCCCGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGGACGAGAAGAACCTCAGGCATGCGTCCCTTGGGGGAAAGAATCTGCCGAAAACTGTGCGGGCCCGGCATGATTCTTTCCCGGCGGGGTATTAAAGAGCGCGATGCaa                                                                                                                >  1:52617/1‑249 (MQ=255)
                                                                                                                                     cATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATAGCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCt                                                                                                           <  2:52617/251‑1 (MQ=255)
                                                                                                                                                   tGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGCAGTGTACAGGCGGTTGGGAATACATTGTGCTAGATCCCGGGGGGCGGTACCCGCCCAAAGCTTCCAAACAGGCACCAGAAGAACCTCCCGCATGCGTCCATTGGTGGAAAGAATATGCCGGAAACGGTGCGGGCCAGGCCTGCTCCCTTCCCGGCGGGGTATTTAAGAGGCAGATGCAACGTATGGGACaaaggaccg                                                                                             >  1:202487/1‑244 (MQ=255)
                                                                                                                                                    gTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGt                                                                                            <  1:86192/251‑1 (MQ=255)
                                                                                                                                                    gTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGGGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCCTTTGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGGGTATTAAAGAGCGCGATGCAACGTCTGGAAAAAGGTGACgg                                                                                            >  1:164024/1‑250 (MQ=255)
                                                                                                                                                           gACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGGTTCGGGGGGGGGGTACCCGCAAAAAGATTCATATCAGGCACGAGAAGAACCTCAGGCAAGCGTCCATTGGTGGGAAAAAACTGCCGAAAACTGTGCGGGCCAGGCAGGATCCTTTCCAGGCGGTGGATTAAAAAGCGCGCTGGCACGTCTGGAACAAGGTggcggtgttccc                                                                                     >  2:177088/1‑240 (MQ=255)
                                                                                                                                                            aCTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGGGGCGTTACCCGCAAAAAGCTTCCTATCAGGCACGAGAAGAACCTCCGGCATGCGTCCCTTGGGGGGAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCCGGCGGGGTATTAAAAAGGGCGATGCAACGTCTGGAAACAGGGGACGTTGTCAccc                                                                                    >  1:120150/1‑250 (MQ=255)
                                                                                                                                                              ttCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAAAACGTTGTGCTTGTTCCGGGGTGGCGTTATCAGCCAAAAACTTCATATCAGGGACCAGAAAAACCTCCGGCCTGCGTCCCTTGGTGGAAAGAATCTTCCGGAAAATTGTCGGGCCCGGCATGATCCTTTCCTGGCGGGGTATTTAAAAGCGCGATGCAAAGTCTGGAACaaagggaccttgttaccca                                                                                   >  1:114875/1‑233 (MQ=255)
                                                                                                                                                                    aTCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGGGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCCTTGGTGGGAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAAAGCGCGGTGCAACGTCTGGAAAAAGGTGACGTTTTTACCGAAACTc                                                                              >  1:388381/1‑249 (MQ=255)
                                                                                                                                                                    aTCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGGTGTGGTAGTTTCGGGGGGGCGTTACCAGCAAAAAGATTCATATCAGGCAAGGGAAGAACCTCAGGCATGCGTCCATTGGTGGGAAGAATCTGCCGAAAACTGTGGGGGGCAGGGATGGTCCTTTTCAGGCGGGGGATTTAAGAGGGGGGGGCAACGTCTGGGACAAGGGGACGtttgcccccaaactccg                                                                            >  1:45441/1‑236 (MQ=255)
                                                                                                                                                                       gCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCtt                                                                         <  1:224971/251‑1 (MQ=255)
                                                                                                                                                                           aaGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGccc                                                                     <  1:164149/251‑1 (MQ=255)
                                                                                                                                                                                          gcgTCGCCGAGTGAAGTGTACAGGCGGTTTGCAATACGTTGTGCTAGGTCCCGGGGGGCGTTACCAGAAAAAAGCCTCATATCAGGCACGGGAAGAAACCCCGGCATGCGTCCATTGGGGGGAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCCGGCGGGGTATTAAAGAGGGCGGTGCAAAGTCTGGAACAAGGGGGCGTTGTCACCGCAACTCCGCTTGCCCGGCTTAAAGCCTggg                                                      >  1:464640/1‑250 (MQ=255)
                                                                                                                                                                                              cGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGctct                                                   <  1:370511/250‑1 (MQ=255)
                                                                                                                                                                                                ccGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCAGGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCtgtg                                                <  1:486945/251‑1 (MQ=255)
                                                                                                                                                                                                caggttGAAGTGTACAGGCGATTGGCAATACGTTGTTCGGGGTCCGGGGTGGCGTTACAAGCAAATAGCTGCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGTGCCAGGCATGATCCTTTCCAGGCGGGGTATTAAAGAGCGCGATGCAACGTCTGGATCAAGGTGACGTTGACACCGAAACACAACTTGCCCGGATTACAGCATGGCTCtgtg                                                <  2:419126/246‑1 (MQ=255)
                                                                                                                                                                                                 cGAGTGAAGGGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGc                                               <  2:120150/251‑1 (MQ=255)
                                                                                                                                                                                                    gTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAAt                                            <  1:488721/251‑1 (MQ=255)
                                                                                                                                                                                                       aaGTGTACAGGCGGTTGGCAATCCGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATAGCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATggg                                         <  2:18716/251‑1 (MQ=255)
                                                                                                                                                                                                        aGTGTACAGGCGGTTGGCAAGACGTTGTGCTAGGTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAATAGCGCGATGCAACGTCTGTAACAAGGTGACGTTGTCACCGAAACTGAGCTGGCCCGGGTTAAATCATGGCTCTGTACAATgggg                                        <  2:491525/251‑1 (MQ=255)
                                                                                                                                                                                                                    gTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAAc                                                                                <  1:19520/199‑1 (MQ=255)
                                                                                                                                                                                                                    gTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGGTGTCACCGAAAc                                                                                >  2:19520/1‑199 (MQ=255)
                                                                                                                                                                                                                      tGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGcc                           <  1:401051/250‑1 (MQ=255)
                                                                                                                                                                                                                      tGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGcc                           <  1:333872/250‑1 (MQ=255)
                                                                                                                                                                                                                        gCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCCACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCttt                        >  2:378880/1‑251 (MQ=255)
                                                                                                                                                                                                                         cAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCCCCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTc                       <  1:116499/251‑1 (MQ=255)
                                                                                                                                                                                                                          aaTACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCg                      >  2:415153/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaTACGTTGTGCTAGTGCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCg                      <  2:164451/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCcac                  >  1:92445/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ccGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAg      <  1:453679/250‑1 (MQ=255)
                                                                                                                                                                                                                                               gggTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCAt  <  1:175522/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          
AAGTAGGGGATAACAGCGGTGATACGACCTGCGGAAGCACGACGCAGGGCATCAACCATAACGACTAATTCCATCAGGTTGTCGTTAGTAGGGGCACAAGTGGACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTGGCAATACGTTGTGCTAGTTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGGGGAAAGATTAGCGCCTTTCGCCACAAAGCCATTGAGCCAT  >  NC_000913/1261601‑1262089

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