Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6490TG32.5% 72.1 / 19.3 83G65G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (14/42);  new base G (23/4);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 2.48e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAG  >  NC_000913/433419‑433863
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                      
tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCCCATGCCGAAGGACACGCGTTGCGTCTGGCGGGGGTGAAAAACAAAAGGTggagcc                                                                                                                                                                                                                  >  2:1494977/1‑233 (MQ=255)
        agagATCGGGCGCCTGAGGGCATAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTCCCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTAAACCACAGGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTca                                                                                                                                                                                                        <  2:254881/239‑1 (MQ=255)
               gcgcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                 <  1:1534709/239‑1 (MQ=255)
                   gcCTGGAGGACTAAGCCGGGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAAGGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGGCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAAcc                                                                                                                                                                                              <  1:285171/238‑1 (MQ=255)
                     cTGGAGGACTAAGCCGTGCAGGACGAATATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTCCCACGCATCCCAACCCAAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTg                                                                                                                                                                                           <  1:448339/239‑1 (MQ=255)
                         aggaCTAAGCCGTGCAGGACAGGTATTACAGGGCGAGGGCGCTAAAGCTGGGGCAACGAGGACGTTTTACCACGCACCCCAACCCGAATGGCGGGTGCGTCATTGTCAAGGATGGAGAAATTGTCGGTGAGGGTTACCACCATCGGGCGGGTGAACCACATGCCAAGGTACCCGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAg                                                                                                                                                                                        <  2:335197/238‑1 (MQ=255)
                               aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGTGCGACCCCCCATGTCACACTCGAACCCCGTAGCcatcat                                                                                                                                                                                 >  1:796039/1‑239 (MQ=255)
                                     tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGa                                                                                                                                                                                                                    <  1:235040/198‑1 (MQ=255)
                                           aCGATTATTACATGGCGCGGGCGTTAAAGGTGGCGCAACGGGGACGTTTTACCCCGGCACCCAGCCCGAATTTGGGTTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGACCCACAGGCGAAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTa                                                                                                                                                                          <  2:64568/234‑1 (MQ=255)
                                           aCGAGTATTCAATGGCCCGGGCGCTAAAGCTGGCGCAAGGAGGACGTTTTCCCACGCATCCCACCCCGAATGTCGGGGGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCCCCAACGGGGGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAGGGGGCGACCGcc                                                                                                                                                                                                               <  2:1031517/197‑1 (MQ=255)
                                           aCGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACg                                                                                                                                                                        <  1:1766531/236‑1 (MQ=255)
                                            cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAAAGAGGACGTTTCACCACGCACCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATGGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                     <  1:379627/238‑1 (MQ=255)
                                            cGAGTATTAAAGGGCGCGGGCTCTAAGGCGGGCGCAACGAGGACTTTTTCCCACACATCCCACCCCAAATGTGGGGGGCGTCATTGGCAAAGATGGCGAAATTGTGGGGGAAGGTTCCCACCACCGGGCGGGTGAACCACATCCCGAAGTACACGCGTGGCGTATGGCGGGTAAAAAACCAAAAGGTGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                     <  2:17844/238‑1 (MQ=255)
                                            cGAATATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                     <  1:1303127/238‑1 (MQ=255)
                                               gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGGGCCCGCCtcttg                                                                                                                                                                                                          >  2:235040/1‑194 (MQ=255)
                                               gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGGGCGACCGCCTATGTCACACTCGCAACCTGTAGCCATCCTGGt                                                                                                                                                                              >  1:1192323/1‑226 (MQ=255)
                                               gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACCCCCGAACCCTTTAGCCATCATGGGCGTACGCCCCCgc                                                                                                                                                                 >  1:445074/1‑238 (MQ=255)
                                               gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGGCACACTCGAACCCCGTAGCCCTCAACGTCGTACGCCAc                                                                                                                                                                    >  1:64568/1‑236 (MQ=255)
                                                tATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCAGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCt                                                                                                                                                                                            <  1:1825900/211‑1 (MQ=255)
                                                tATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCAGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                         <  1:673501/214‑1 (MQ=255)
                                                 aTTACAGGGCGCGGACGCTAAAGCGGGCGCAACGAGGACGTTTTACCAGGCATCCCAGCCCGAATGTCGGGTGCGTCATTGTCAAAGAGGGAGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCAATTGCCGAAGTACGCGCGTGCGGAATGGCGGATGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                <  2:445074/238‑1 (MQ=255)
                                                    aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCgg                                                                                                                                                                 >  1:278756/1‑233 (MQ=255)
                                                     cagggCGCGGGCGTTAAAGCTGGCGAAACGAGGACTTTTTCCCACGCATCCCACCCCGAATGTCGGTTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGACCCACATGCCGAAGTCCCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                         <  2:82107/206‑1 (MQ=255)
                                                     cATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCccac                                                                                                                                                                                                     >  1:1031517/1‑197 (MQ=255)
                                                          cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGTCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGc                                                                                                                                                       >  2:199854/1‑238 (MQ=255)
                                                          cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGGGACCGCCCATTTCACACTCCAACCCCGTAGCCATCATGGTCGTACGCCACCGTGCTGTGAAGc                                                                                                                                                       >  1:462287/1‑238 (MQ=255)
                                                          cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCCCACTCGAAACCTGTAGCCATCATGGCCGTACGCCACCCTGCTGTGCCGc                                                                                                                                                       >  1:1185923/1‑238 (MQ=255)
                                                          cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTATCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCCGTAGCCATCATGGTCGTACGCCCCCGTGCGGTGACGc                                                                                                                                                       >  1:1693702/1‑238 (MQ=255)
                                                               gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTTCACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAACCATCATgg                                                                                                                                                                               >  1:82107/1‑209 (MQ=255)
                                                                         tGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGGCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgccgc                                                                                                                                           <  1:308222/235‑1 (MQ=255)
                                                                         gagcgcAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTGGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCTTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTgg                                                                                                                                        <  1:1112930/236‑1 (MQ=255)
                                                                                 cGAGGACTTTTTACCACCAATCCCAACCGGAAGGTCGGGGGGTTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTAACCCACATGCCGAGGAACACGCGGGGCGTATGGCGGGTGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTg                                                                                                                                                                                           <  2:1795685/179‑1 (MQ=255)
                                                                                     gACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCCGTAGCCATCATGGTCGTACGCCACCGGGCTGTGAAGCACTCATCGCCGCTGGCGTAGCGCGCGTg                                                                                                                           >  1:1004509/1‑239 (MQ=255)
                                                                                         ttttACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAAGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGGAGCGCGCCTGGtt                                                                                                                        >  1:590221/1‑238 (MQ=255)
                                                                                           ttACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTAACACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTCTGTCACACTCGAACCCTGTAGcccccat                                                                                                                                                                                 >  1:1795685/1‑174 (MQ=255)
                                                                                                  gCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCCCTCATCGCCGCTGGCGTAGCGCGCGGGGTTGCCTCGATGc                                                                                                              >  1:939938/1‑239 (MQ=255)
                                                                                                   cATCCCACCCGAAAGGTCGGGTGCGTCATTGGAAAAGAGGGCGAAATTGTCGGTGAAGGTTACCCCCAGCGTGCGGATGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                             <  2:590221/239‑1 (MQ=255)
                                                                                                   cATCCCACCCCGAATGGCGGGGCGGTCATTGTCAAAGATGGCGAATTTGTCGGGGAAGGTTCCCACCACCGTGGGGGTGAACCACATGCCGAAGTCCCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                             <  2:36835/239‑1 (MQ=255)
                                                                                                   cATCCCAACCCGAATGTCGGGGCCGTCATTGTCAAAGATGGCAAAATTGTCGGGGAAGGTTACCCCAAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                             <  2:939938/239‑1 (MQ=255)
                                                                                                   actCCCACCCGGAAGCTCGGGTACTCCATTGACAAAGGTGGCGAAATGGTGGGTAAGGGTTCCCCCCACGGTGGGGGTGACCCACAGCCCGAAGTACCCTCTTTGCGTATGGCGGGTGAAAAGGCCAAGGGTGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                             <  2:1211666/237‑1 (MQ=255)
                                                                                                    aTCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                             <  1:741649/238‑1 (MQ=255)
                                                                                                      cccAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGa                                                                                                          >  1:36835/1‑239 (MQ=255)
                                                                                                           cccGAATGTCGGGTGGGTCATTCGCAAATGTGTCGAAATTGCCGGTGAAGGTCACCCCCAGCGTGCGGTTGAACCGCTTGCCGAAATCAACGAGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTAGGCCCGCCTATGTCACACTCGAACCCTGt                                                                                                                                                                                          <  1:1078728/154‑1 (MQ=255)
                                                                                                                      ggggCGTAATTGTCAAAGAGGGGGAAATTGGGGGGGAAGTTTACCACCAACGGGCGGGTAAACCAAAGGCCGAAGTACACGCGTTGCGTAGGGGGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCg                                                                                                                                      <  2:1401522/192‑1 (MQ=255)
                                                                                                                         gCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTAAAAAGGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                        <  1:930949/238‑1 (MQ=255)
                                                                                                                                tGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                     <  1:905064/154‑1 (MQ=255)
                                                                                                                                tGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGt                                                                                                                            >  1:1401522/1‑195 (MQ=255)
                                                                                                                                ggaaaaaGATGGAGAAATGGTGGGTGATGGTGACCACCAGCGTGGGGGGGACCCACATGCGAAAGTACACGCGTTGCGTGTGGGGGGTGAAAAACCCAAGGTTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                 <  2:441018/234‑1 (MQ=255)
                                                                                                                                    aaaGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCCATGTCACACCCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGCGGACt                                                                            >  1:1178784/1‑239 (MQ=255)
                                                                                                                                          ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGACAAAGGGGCGACCGCCCATGTCCCACTCGCACCCCGTAGCCATCATGGTCGTAAGCCCCCGGGCtgtg                                                                                                                                                           >  2:905064/1‑154 (MQ=255)
                                                                                                                                            cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGGGGTTGCCCCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGGGGACttt                                                                          >  1:1828427/1‑233 (MQ=255)
                                                                                                                                             gAAATTGTCGGGGAAGGTTACCCCCAACGGGCGGGGAAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTc                                                                    <  1:125162/238‑1 (MQ=255)
                                                                                                                                                  tGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGGGCGACCCCCCATGTCACACTCCAACCCTGTAGCCATCATTGTCGTACGCCCCCGCGCTGTTTCCCACCCCTCGCCCCAGGCGTAGCGCGCGGGGGTTCCCCGCTGCAAGATCCCAACCCGCCCGTCGCGGGGCGCGGCCTTTACCGCCTGCac                                                               >  2:740489/1‑237 (MQ=255)
                                                                                                                                                             gTTACCCCCAAGGTGGGGGTAAACCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCa                                                     <  2:278756/237‑1 (MQ=255)
                                                                                                                                                                     aaaCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCAtt                                                   <  1:1354646/230‑1 (MQ=255)
                                                                                                                                                                       aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                               >  1:1415023/1‑105 (MQ=255)
                                                                                                                                                                       aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                               <  2:1415023/105‑1 (MQ=255)
                                                                                                                                                                        cGTGCGGGTGACCCCCATGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGGCCGCCTATGt                                                                                                                                                                                                          <  2:1072443/77‑1 (MQ=255)
                                                                                                                                                                         gTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                     <  1:370687/113‑1 (MQ=255)
                                                                                                                                                                            cGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTca                                                                                                                                                                                                        <  2:137230/75‑1 (MQ=255)
                                                                                                                                                                            cGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTca                                                                                                                                                                                                        >  1:137230/1‑75 (MQ=255)
                                                                                                                                                                             ggttGACCCACATGCCGAAGTACGCGCGTTGCGTATGGCGGTTGAAAAACCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACg                                                                                                                                                                        <  2:1834923/103‑1 (MQ=255)
                                                                                                                                                                               gTGAACCCCATGCCGAAGTCCACGGGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                   <  1:1158577/189‑1 (MQ=255)
                                                                                                                                                                               gTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCCCCCATGTCCCACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTTACGCACTCATCGCCGCTGGCGGAGCGCGCGCGGGTGCCTCCATGCCAGATCCTAACCCGCAGGGCGCCGGGCGTGGGCTTTTCCCGCCTCCACAGGCTGGCATTGTCgctcgccccct                                     >  2:1354646/1‑225 (MQ=255)
                                                                                                                                                                                  aacaacATGCGAAAGTACGCGGGTGGCGTAGGGCGGGGAAAAAACCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGCCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg                              <  2:1185923/239‑1 (MQ=255)
                                                                                                                                                                                  aaCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                                >  1:1072443/1‑77 (MQ=255)
                                                                                                                                                                                  aaCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg                              >  1:1638805/1‑239 (MQ=255)
                                                                                                                                                                                   aCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGCGACCGCCTATGTATCACTCGAACCCTGTAGCCATCAAGGTCGTACGCCACCCGGCtgtg                                                                                                                                                           >  2:370687/1‑113 (MQ=255)
                                                                                                                                                                                       caTGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGGCAAAGGGGCGACCGCCCATGGCACACTCGACCCCTGTAGCCATCATGGGCGTACGCCACCGGGCt                                                                                                                                                              >  1:1834923/1‑106 (MQ=255)
                                                                                                                                                                                         tGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGTGCGACCGCCCATGTCACACTCGAAACCTGTAGCCATCATGGTCGTACGCCACCCTGCTGTGACGCCCTCATCGCCCCCGGCGGAGCGCCCGTGGGTTCCTCCATGCAAGAACCTAACCCGCCGGGCGCCGGGGGGGGAGtttt                                                                         >  2:1158577/1‑188 (MQ=255)
                                                                                                                                                                                            cGAAGTCCCCGTGTTGCGTAGGGGGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAg                    <  2:1828427/239‑1 (MQ=255)
                                                                                                                                                                                               aGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCCAAGGGGCGACCCCCCATTTCACCCACCAAACCCGTAGGCATCATGGGCGTAAGGCACCGCCCTGTGACGGACCCCTCGCCCCCGGCGTAGCGCGCGCGGGTGCCCCTGTTCCAGATCCTTACCCGCCCGCCGCCGGGCCGCGGCTTTCCCCGCGGCAACACGCCGGCCCTGCGGTCCGCCACGGCGCGGTGATTAGTGGAGccc                 >  2:574692/1‑238 (MQ=255)
                                                                                                                                                                                                gTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                      <  1:287094/217‑1 (MQ=255)
                                                                                                                                                                                                 tACCCGCGTTGCGTAGGGCGGGGAAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgccg                                                                                                                                            <  2:1328189/114‑1 (MQ=255)
                                                                                                                                                                                                 tACACGCGTTGCGTATGGCGGGGGAAAAAGCAAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGc                                                          <  1:442910/196‑1 (MQ=255)
                                                                                                                                                                                                    ccgggTTGCGTATGGGGGGTGAAAAACCAAAAGGTGGCACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa            <  1:844053/235‑1 (MQ=255)
                                                                                                                                                                                                           gCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGGCGCACCCATCGCCGCTGGCGTAgc                                                                                                                                  >  1:1328189/1‑114 (MQ=255)
                                                                                                                                                                                                           gCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCCCTCTCACACTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGTGCTGTGACGCACTCATCGCCGCCGGCGTCGCGCGCGTGGTTGCCTCGATGCAAGAGCCTAAACCCCAGGGCGCGGGGCGGGGGCTTTa                                                                         >  2:104391/1‑171 (MQ=255)
                                                                                                                                                                                                             gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACttt                                                                          <  1:604269/168‑1 (MQ=255)
                                                                                                                                                                                                                      gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAg  >  1:1561913/1‑231 (MQ=255)
                                                                                                                                                                                                                       tGAAAAAACCAAAGGTGCGACCCCCCATGTCACACTCGAACCCCGTAGCCATCAAGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGGAGCGCGCGTGGGTGCCTCGGTGCAAGAACCCAACCCGCAGGTCGCCGGGCGCGGACTTTACCGTCTGCa                                                                >  2:604269/1‑168 (MQ=255)
                                                                                                                                                                                                                                  aaaGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                        <  1:63805/134‑1 (MQ=255)
                                                                                                                                                                                                                                  aaaGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                        >  2:63805/1‑134 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                      
TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAG  >  NC_000913/433419‑433863

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