Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6340TG35.2% 82.2 / 23.4 91G60G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (14/43);  new base G (30/1);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.84e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCG  >  NC_000913/433400‑433848
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                       
tttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGTCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACACCATGCCGCAGTACCCGCGTTGCGTCTCGCGGGTGaaa                                                                                                                                                                                                                     >  2:1955517/1‑239 (MQ=255)
tttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGCAGTACACGCGTTGCGTATGGCGGGTGgaa                                                                                                                                                                                                                     >  1:1272308/1‑239 (MQ=255)
tttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTCTGGCGGGTGaaa                                                                                                                                                                                                                     >  1:1544732/1‑239 (MQ=255)
tttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGaaa                                                                                                                                                                                                                     >  1:459811/1‑239 (MQ=255)
tttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGaaa                                                                                                                                                                                                                     >  1:1676763/1‑239 (MQ=255)
     aaGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACCCCATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGca                                                                                                                                                                                                                >  1:2290067/1‑238 (MQ=255)
          atCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAGAAAACCcaacg                                                                                                                                                                                                           >  1:63712/1‑237 (MQ=255)
              aggAATTGGGCGAAGAGCTCGCGCGCCTGGAGGACTAAGCCGGGCAGGAGGATTTTTACATGGCGCGGGCGCTAAAGCGGGCGCAACGAGGACTTTTTCCCCCGCACCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTCCCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCg                                                                                                                                                                                                       <  2:2290067/237‑1 (MQ=255)
                  aTTTGCGGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGAGGAGTATTACATGGCGCGGGCGCAAAAGCTGGCGCAACGAGGACGTTTTCCCACGCCTCCCACCCCGAATGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGTTTCCCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGGGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAc                                                                                                                                                                                                     <  2:806861/237‑1 (MQ=255)
                   tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAAGGGCGGCCg                                                                                                                                                                                                   >  1:2020255/1‑238 (MQ=255)
                   tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCAGAAGGACACGCGTTGCGTATGGCGGGCGAAAAAGCCAAAAGCGAGACCg                                                                                                                                                                                                   >  2:1935651/1‑238 (MQ=255)
                    tgggCGAAGAGATCGCGCGCATGAGGGACGAGGCCGGGCGGGACGGGTATACAAGGGCGGGGGCGCTAAAGCGGGCGCAAGGAGGACGATTTCCCACGCATCCCACCCCGAAGGTGGGGTGCGTCATTGTAAAAGAGGGCAAAATGGTCGGTGAAGTTTACCACCACCGTGCGGGTGAACCACTGGCGGAAGTCCCCGCGTTGCGTAGGGGGGGTGAAAAAGCCAAAGGTGCGACCGcc                                                                                                                                                                                                 <  2:1893563/237‑1 (MQ=255)
                     tGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGCGACCGccc                                                                                                                                                                                                >  1:629066/1‑238 (MQ=255)
                           agagATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACAACAGGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCCAAAGTGGGAcccactcgttac                                                                                                                                                                                          >  2:2094346/1‑229 (MQ=255)
                                    gcgcCTGAGAGACTAAGCCGGGCAGGAGGAGTATTACGTGGCGCGGGCGCTAAAGCTGGCGCAAGGAGTCGGTTTTCCCCCGCACCCCAACCCAAATGTGGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAGGTTTCCCACCAACGGGCGGGTAAACCACAGGCCGAAGTACCCGCGTTGCGTATGGGGGGTAAAAAACCAAAAGGTGCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                 <  2:1413410/238‑1 (MQ=255)
                                            gggaCTAAGCCGTGCAGGACGATTATTACATGGCGCGGGCGCAAAGGCGGGCGAAACGAGGACGTTTTCCCCCGCACCCCACCCCGAATGTGGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGACCCACATGCCGAAGTACACGCTTTGCGTATGGGGGGTAAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCg                                                                                                                                                                                   <  2:279387/227‑1 (MQ=255)
                                                cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGGGACCGCCCATTGTACACTCGAACCCTGTAAccct                                                                                                                                                                     >  1:1699916/1‑236 (MQ=255)
                                                  aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                                                                                                                                                                                                                >  1:1509276/1‑194 (MQ=255)
                                                  aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                                                                                                                                                                                                                <  2:1509276/194‑1 (MQ=255)
                                                      cGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAg                                                                                                                                                                         >  1:279387/1‑228 (MQ=255)
                                                       gTGCGGGACGGTTAATACAGGGCGCGGGCTCTAAGGCGGGCGCAACGAGGGCGTTTTCCCACACATCCCAACCCAAATGGCGGGGGCGGCATTGTAAAAGATGGGGAAATTGTCGGTGAGGGTTACCCCCAACGTGGGGGTGACCCACATGCCGAAGTACGCGCGTTGCGTAGGGCGGGTGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                              <  2:1544732/238‑1 (MQ=255)
                                                           aGGACGGGTTTTACAGGGCCGGGGCTGTAAAGCTGGCGCAACGAGGACGTTTAACCCCGCATCCAAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAGGGTTCCCCCCAACGTGGGGGTGACCCACATGCCGAAGTACACGCGTTGCTTATGGGGGGTAAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                          <  2:1836599/222‑1 (MQ=255)
                                                             agggAGTATTACATGGCGCGGGCGCAAAAGCGGGCGCAACGAGGACGTTTTACCACGCACCCCAACCCGAATGGCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAAGGTGCGGGTGAACCACAGGCCGAAGTCCACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTa                                                                                                                                                           <  2:33809/232‑1 (MQ=255)
                                                               cGAGTATTACGGGGCGCGGGCGCTAACGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTGGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGACCCACAGGCAGAAGTCCACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                      <  2:1828790/238‑1 (MQ=255)
                                                               cGAATATTACTTGGCGGGGGCGCTAAAGCTGGCCAAACGGGGCCTTTTTCCCACCCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTAAAGGTTCCCACCAACGGGGGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                      <  2:2020255/238‑1 (MQ=255)
                                                                 aGTATTACAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTCCCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                              <  1:1612482/228‑1 (MQ=255)
                                                                  gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCGCCTCGGCGCCACTCCAACCCCGGAGCACACATGGTCGTAAGCCACCgc                                                                                                                                                  >  1:2217828/1‑238 (MQ=255)
                                                                   tATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTg                                                                                                                                                                                                         <  1:2093901/184‑1 (MQ=255)
                                                                   tATTACATGGCGCGGGCGCAAAAGCGGGCGCAACGAGGGCGTTTTACCACGCACCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCAAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAACCCAAAGGTCCGACCGCCTATGTCACACTCGAACCCTGTAGc                                                                                                                                                                        <  2:348749/216‑1 (MQ=255)
                                                                    agcacgtGGCGCGGCCGCAAAGGCGGGCGCAACGGAGCCTTTTTCCCACGCCACCCACCCCGAAGGTCGGGTGCGTCATTGTCAAAGAGGGCGAATTTGTCGGTGAAGGTTACCACCAGCGGGGGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                 <  2:2258575/232‑1 (MQ=255)
                                                                     ttACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGTGCGACCGCCTATTTCACACTCCAACCCTGTAGCCACAATgg                                                                                                                                                                >  1:1836599/1‑222 (MQ=255)
                                                                       aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGGGGGACCGCCCATGTCACACTCGAACCCCGTAGCCATCATGGACGTACGCCAcccg                                                                                                                                                  >  1:38809/1‑232 (MQ=255)
                                                                           ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCCCCCATGTCCCACCCGAACCCTGTAGCCATCATGGGCGTACCCCACCCTg                                                                                                                                                 >  1:33809/1‑231 (MQ=255)
                                                                           ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCCCCCATGTCACACCCGAACCCCGTAGCCATCAAGGTGGTACGCCCCCCAGCtgtg                                                                                                                                            >  1:586813/1‑236 (MQ=255)
                                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAACCAAAGGTCGACCCCCtg                                                                                                                                                                                               >  2:2093901/1‑183 (MQ=255)
                                                                               cgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGACAAAGGGGCGACCGCCTATGTCACACTCGAAACCTGTAGGCATCATGGTCTGACCCCACCGGGCGGTGGCGCaa                                                                                                                                      >  1:2197743/1‑237 (MQ=255)
                                                                               cgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAACCCAAAGGGGCGCCCGCCTATGTCACACTCGAACCCCGTAGCCATCATGGTCGGACGCCACCGTGCTGCGACGCAc                                                                                                                                      >  1:777467/1‑238 (MQ=255)
                                                                                  gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGTGCGACCGCCTATGTCACACTCGAACCCCGTAGCCATCATGGGCGGACGCCACCGTGCTGGTACGCACTCa                                                                                                                                   >  1:1266071/1‑238 (MQ=255)
                                                                                                           gTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCCCAGGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                            <  1:442197/157‑1 (MQ=255)
                                                                                                                 cccgccATCCCACCCCGAATGTCGGGTGCGTCAGTGTAAAAGAGGGCGAAATGGGCGGTGAAGGTTCCCACCACCGTGCGGGGAAACCACATGCCGAGGTACCCGGTTTGCGTATGGCGGGTGAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCg                                                                                                   <  2:629066/234‑1 (MQ=255)
                                                                                                                     gCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGGCCGCCCCTGTcaccaccgaa                                                                                                                                                                                  >  2:442197/1‑150 (MQ=255)
                                                                                                                      cATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCCCCGGCGTAGCGCGCGGGGTTGCCTCCATGCa                                                                                              >  1:847998/1‑239 (MQ=255)
                                                                                                                       accccaacccaaATTGCGGGGGGGTATATGTAAAAGATGGGGAAATTGTCGGTGAAGGTTACCCCCAACGTGGGGGGAAACCACAGGCAGAAGTACCCGGGTTGCGTATGGGGGGGAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                              <  2:670555/236‑1 (MQ=255)
                                                                                                                       aTCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                                                                                                                                                                                                    <  1:61692/137‑1 (MQ=255)
                                                                                                                            aaCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAGGTGCGACCGCCCATGTCACACCCCAACCCCGTAGCCATCCTGGTCGTACGCCCCCCGTCTGTGCCGCACCCATCGCCCCCGGCGTAGCGCGCGCGGTTGCCCCCGTGCCAGATc                                                                                         >  2:2119153/1‑238 (MQ=255)
                                                                                                                                 gAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGGGACCCCCCATGTca                                                                                                                                                                                          >  2:61692/1‑137 (MQ=255)
                                                                                                                                  aaTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCCCCCATGTCCCACCCCAAACCCTTAACCCACCAGGGCGTCACCCACCGCGCGGTGAGGCCCCCACCCCCGCCTGCGTAGCGGGCGGTGGTGCCTCGATTGACGATccccacc                                                                                   >  2:1268621/1‑233 (MQ=255)
                                                                                                                                   aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAGGTGCGACCGCCCATGTCACACTCGAACCCCGTAGCCACCATGGTCGTACGCCACCGTGCTGTGACGCACTCCTCGCCCCTGGCGGCGCGCGCGGGGTGTCCTCCATGCCAGATCCCAAccc                                                                                  >  2:2081115/1‑238 (MQ=255)
                                                                                                                                    tgtgGGGTGCGTCATTGGAAAAGTGTGCGAAATGTGCGGTGGAGGTTACCACCAGCGTGCGGGTGACCCACAGGCGGAGGTACACGTGTAGCGTATGGGGGGTGAAAAACCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGc                                                                                <  2:1326779/239‑1 (MQ=255)
                                                                                                                                                 atggacAAAGAGGGCAAATGTGGCGGTGAGGGTCACCGCAAAGGGGGGGGAAAACCACAGGCCGAAGTACCCGCGGTGCGTAGGGGGGGTAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                    <  2:38809/233‑1 (MQ=255)
                                                                                                                                                  ttGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCGCATGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                  <  1:1268621/239‑1 (MQ=255)
                                                                                                                                                  ttGACAAAGATGGAGAAATTGTCGGTGAAGGTTCCCACCAAGGGGCGGGTGAACCACATGCCGAAGTACCCGGGTTGCGTATGGCGGGTGAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                  <  2:847998/239‑1 (MQ=255)
                                                                                                                                                   tttCAAGGATGCCGAAATTGTCGGTGAGGGTTCCCGCCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGTCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                     <  1:1529267/201‑1 (MQ=255)
                                                                                                                                                     tCAAAGATGGCGAAATTGTCGGTGAGGGTTACCCCCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGt                                                                                                             <  2:427380/193‑1 (MQ=255)
                                                                                                                                                             ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                            <  2:1376321/107‑1 (MQ=255)
                                                                                                                                                             ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                            >  1:1376321/1‑107 (MQ=255)
                                                                                                                                                             ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGTGCGGCCGCCCCTGTCACCCCCCAACCCCGTAGACATCACGGGCGTACCCCCCCGTGCTGTGTCGCCCCCACCGCCCCCCGCGCACCGCGCGGGGTTGCCCTCATGCACGAGCCCACACCCCACGTCGCCGGGCgggggggttact                                                        >  2:1471686/1‑230 (MQ=255)
                                                                                                                                                             ggCGAAATTGTCGGGGAGGGTTACCCCCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                     <  1:37367/193‑1 (MQ=255)
                                                                                                                                                               cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGGGTAGCGCGCGGGGTTGCCTCg                                                                                                   >  1:427380/1‑193 (MQ=255)
                                                                                                                                                                    ttGTCGGTGAAGGTGACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCaa                                                <  1:2081115/239‑1 (MQ=255)
                                                                                                                                                                       tCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCCATGTCACCCCCGAAACCCGTAGCCACCATGGGCGTACGCCACCGTGCTGTGACGCACCCATCGCCGCTGGCGTAGCGCGCGGGGTTGCCCCGCTggcaga                                                                                           >  2:37367/1‑188 (MQ=255)
                                                                                                                                                                          gTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGAACCCGCTTTGCGTATGGCGGGTAAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                          <  1:1740907/111‑1 (MQ=255)
                                                                                                                                                                           agAAGGTTACCACAAGCGTGCGGGTGACCCACATGCCGAAGTACCCGCTTGGCGTAGGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGt                                                                                                                                                                           <  2:881786/108‑1 (MQ=255)
                                                                                                                                                                                    ccaccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCCCCCATGTCACACTCGAACCCCGTAGGCAACATgg                                                                                                                                                                >  2:1740907/1‑111 (MQ=255)
                                                                                                                                                                                     caccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGTGGGACCGCCCATGTCCCACTCGGACCCTGTAGCCAACATg                                                                                                                                                                 >  1:881786/1‑109 (MQ=255)
                                                                                                                                                                                          aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGGGCGAACGCCCATg                                                                                                                                                                                             >  1:1204498/1‑77 (MQ=255)
                                                                                                                                                                                           ggggggggTGACCCACAGGCCAAAGTCACCGCGTTGCGTAGGGCGGGTGAAAAAGCAAAAGGGGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                         <  1:350098/234‑1 (MQ=255)
                                                                                                                                                                                           cGTGCGGGTGAACACCAGGCCGAAGTCCCCGCGTTGCGAATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                         <  1:1261089/239‑1 (MQ=255)
                                                                                                                                                                                              gCGGGAGACCCACATGCCAAAGGACAGCGGTGGGTTAGGGAGGGTAAAAAACCCAATGGGCCAACCGCCTATTTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                       <  2:178805/238‑1 (MQ=255)
                                                                                                                                                                                                ggggGAACCACAGGCCGAAGTACACGCTTTGCGTATGGCGGGTAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                     <  1:1280923/238‑1 (MQ=255)
                                                                                                                                                                                                gggTAACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                          <  2:2247044/89‑1 (MQ=255)
                                                                                                                                                                                                  gtaaaCCCCAGGCAGAAGTACCCGCGTTGGGGATGGCGGTAAAAAAAGCCAAAGGGGCGACCGCCTATGTCGCACTCGAACCCTGTAGCCATCATGGTCGTACg                                                                                                                                                         <  1:2058003/101‑1 (MQ=255)
                                                                                                                                                                                                  gTGAACCACATGCCGAAGTACACGCGTTGCGTTTGGGGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                 <  2:1928937/112‑1 (MQ=255)
                                                                                                                                                                                                     aaCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                            <  2:1237890/67‑1 (MQ=255)
                                                                                                                                                                                                     aaCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                            >  1:1237890/1‑67 (MQ=255)
                                                                                                                                                                                                          caTGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCCCCTATGTCACAATCGAACCCTGTAACCATCATgg                                                                                                                                                                >  1:2247044/1‑89 (MQ=255)
                                                                                                                                                                                                           aTGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatc                                                                                                                                                                    <  2:320081/84‑1 (MQ=255)
                                                                                                                                                                                                            tGCCGAAGTCCGCGCGTTGCGTAGGGGGGGTGAAAAAGCCAAGGGGGCGGCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTg        <  1:46666/239‑1 (MQ=255)
                                                                                                                                                                                                            tGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGACAAAAGGGCGGACCCCTATTTCACACACGAACCCTGGAGCCCTCATTGTCGTACACCACCGTGCt                                                                                                                                               >  2:2058003/1‑104 (MQ=255)
                                                                                                                                                                                                            tGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGGAGCCATTATGGTCGTACGCCACCGTGCTGTGACGCa                                                                                                                                       >  1:1928937/1‑112 (MQ=255)
                                                                                                                                                                                                            tGCCGAAGTACACGCGTTGCGTAGGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                      <  1:1108339/97‑1 (MQ=255)
                                                                                                                                                                                                            ggCCGAAGTACACGCGTTGGGTATGGGGGGTGAAAAAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                         <  2:221903/173‑1 (MQ=255)
                                                                                                                                                                                                             gCCGAAGTCCCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                     <  2:144873/145‑1 (MQ=255)
                                                                                                                                                                                                                    tACACGGGTTGCGTAGGGCGGGTAAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGt                                                                                                          <  1:695502/133‑1 (MQ=255)
                                                                                                                                                                                                                    tACACGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                      <  1:509826/89‑1 (MQ=255)
                                                                                                                                                                                                                     acacGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGGCGTAc                                                                                                                                                          >  1:320081/1‑84 (MQ=255)
                                                                                                                                                                                                                      cacGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCCCCCATGTCACACTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGTGCtgtg                                                                                                                                            >  2:1108339/1‑97 (MQ=255)
                                                                                                                                                                                                                      cacGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGGCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGa                                                               >  1:221903/1‑174 (MQ=255)
                                                                                                                                                                                                                       acGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATTTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGa                                                                                           >  1:144873/1‑145 (MQ=255)
                                                                                                                                                                                                                          gtgtggctaggggCGGGAGAAACAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTGGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg  <  2:777467/221‑1 (MQ=255)
                                                                                                                                                                                                                              gCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCCCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACCCATCGCCCCTGGCGTAGCGCGCGGGGTTGCCTCGATg                                                                                                >  2:695502/1‑133 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                       
TTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCG  >  NC_000913/433400‑433848

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