Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6490TG34.9% 72.1 / 19.6 83G65G (GGT→GGG) ribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (16/38);  new base G (24/5);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.50e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCA  >  NC_000913/433433‑433884
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
cgctcgCCTGGAGTCCTAGGCCGTGAAGGCCGAGTATTACAGGGCGGGGGCGCTAAAGCGGGCGCAACGAGGCCGTTTTCC‑CACGCATCCCACCCCGAATGTGGGGTGCTTCATGTTAAAAGAGGGCGAAATTGTCGGTGAAGGTTCCCCCCAACGTGCGGGTGACCCACATGCCGAAGTCCACGCGTTGCGTATGGGGGGTGAAAAAGCAAAAGGGGCGACCGCCTATGTCACACTCg                                                                                                                                                                                                                       <  2:1880541/239‑1 (MQ=255)
               cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGGGAACCCCCATGTCACACCCGAACCCTGTAGccca                                                                                                                                                                                                         >  1:1347674/1‑236 (MQ=255)
               cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAACAAAAAGGCGACCGCCTCTGTCACACACCAACCCCTTAGCcatc                                                                                                                                                                                                        >  2:2261204/1‑239 (MQ=255)
                tacgCCGTGCGGGACGATTATTACATGGCGCGGGCGCTAAAGCGGCCGCAACGAGGACGTTTTCC‑CACCCATCCCAACCCAAATGTCGGGTGCGTCATTGTCAAAGCGGGCGAAATTGTCGGTGACGTTTCCCCCCACCGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCTTATGGCGGTTGAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                                                       <  2:2101975/236‑1 (MQ=255)
                 aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGTGCGACCCCCTATGTCACACTCCAAACCCGTCGAcatcat                                                                                                                                                                                                      >  1:1950683/1‑239 (MQ=255)
                 aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAGCAAAAGGGGGACCCCCCATTTCACAACCCCAACCTGTAGCcaacaa                                                                                                                                                                                                      >  1:2344177/1‑238 (MQ=255)
                     cGTGCAGGCGTATTATTACAGGGCGGGGGATCTAAGGCTGGGCCAAC‑AGGAACTTTTACACCCGCTTCCCATCCGGAAAGTGGGGTGCGTTATGTGAAAAGAGGGCGAATTTGGCGGTAAAGGTGACCCCCAGCGGGCGGGTGACCCACATCCCAAAGTACCCGCTTGGCGTATGGCGGGTGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                                                                  <  2:1271702/239‑1 (MQ=255)
                      gTGCGGGACGAGTATTACATGGCGCGGACGCTAAAGCTGGAGCGACGAGGACGTTTTAC‑CACGCATCCCACCCCGAATGTGGGGTGGGTCATAGACAAAAAGGGCGAATTTGGGGGTAAAGGTTACCACCAACGTGGGGGTAAACCACATGCCGAAATGCACGCTTAGCGTATGGCGGGTGAAAAACCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCAACATGGTCg                                                                                                                                                                                                 <  2:1427387/239‑1 (MQ=255)
                       tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGGGACCCCCCATTTCACACTCGAACCCCGTAGCCATCATGGGCGt                                                                                                                                                                                                >  1:1905590/1‑239 (MQ=255)
                            gACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAACGGGCGCGCGCCTATGTCCCCCACGCACCCCGTTGCCATCATCGTctgaacc                                                                                                                                                                                            >  2:1776854/1‑232 (MQ=255)
                                aGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAc                                                                                                                                                                                                                                        <  1:1277603/190‑1 (MQ=255)
                                            cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAc                                                                                                                                                                                                                                        <  2:177301/178‑1 (MQ=255)
                                            cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGAc                                                                                                                                                                                                                                        >  1:177301/1‑178 (MQ=255)
                                                           tggggaaaCGAGTGCTTTTTAC‑CGCGCATCCCATCCCAAATGTCGGGGGCGTCATTAGCAAGGATGGCGAAGGTGTCGGTGAAGGTTACCACCAAGGTGCGGGTGACCAGCATGACGAAGTACACAGGTTGCGTGGGGCGGGTGAAAACGCAAAGGGGGCGACCGCCTATGTCTCACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgccg                                                                                                                                                                 <  2:430506/228‑1 (MQ=255)
                                                                    gAGGACGTTTTAC‑CCCGCACCCCACCCGGAATGTGGGGGGCGCCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGACCCCCATGCCGAAGTACACGGTTTGCGTATGGGGGGTAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                                                     <  2:2172880/173‑1 (MQ=255)
                                                                        aCGTTTTAC‑CCCGCACCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAGGCCAAAGGTGCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                                                     <  1:591077/169‑1 (MQ=255)
                                                                            tttAC‑CACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGt                                                                                                                                                                                                               <  1:295099/171‑1 (MQ=255)
                                                                              tAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGTGCGACCGCCTATGTCACACCCGAAACCTGTAGcc                                                                                                                                                                                                           >  1:2172880/1‑173 (MQ=255)
                                                                                   aCGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGGGGGACCGCCCATTTTACAACCGAACCCTGTAGcc                                                                                                                                                                                                           >  2:591077/1‑169 (MQ=255)
                                                                                     gCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGGGGGACCGCCCATGTCccacc                                                                                                                                                                                                                         >  1:1146955/1‑152 (MQ=255)
                                                                                       aTCCCACCCCGAATGTCGGTGGCGGCATTGTAAAAGAGGGCGAAATTGTCGGTGAAGGTTACCCCCACCGGGCGGGTGAACCACATGCCGAAGTACACGCTTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                  <  2:2344177/238‑1 (MQ=255)
                                                                                       aTCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGGGACCGCCCATGTTACACCCGAACCCCGTAGCCAACATg                                                                                                                                                                                                     >  2:295099/1‑171 (MQ=255)
                                                                                            aaCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGTTTCCCCCCAACGTGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                               <  1:2241415/140‑1 (MQ=255)
                                                                                                  aaTGTCGGGTGCGTCATTGTCAAAGAGGGCAAAATTGTCGGTGAGGTTTACCACCAACGTGGGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                              <  1:1130896/151‑1 (MQ=255)
                                                                                                   aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCGCCCATGTCACCAACGAAACCCGTAGCCATCAAGGGCGGACGGCACCCGGCCGGGAAGCACACTACGCCGCCGGCGGTACGCGCGGCGGTGGCTCGGCGCAAGAAccgcaacc                                                                                                                      >  2:2362064/1‑232 (MQ=255)
                                                                                                   aTGGCGGGTGCGTCATTGTAAAAGATGGCAAATTTGTCGGTGAAGGTTCCCCCCAACGTGGGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGCGGGTAAAAAAGCAAAGGGTGCGACCGCCt                                                                                                                                                                                                                                   <  2:386050/129‑1 (MQ=255)
                                                                                                        gggTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTAAGGGTTCCCCCCAACGGGCGGGTGACCCACATGCCGAAGTACACGCTTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGcc                                                                                                                                                                                                                                    <  2:150324/123‑1 (MQ=255)
                                                                                                         ggTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAAGTGCGAACGCCCATg                                                                                                                                                                                                                                >  1:1020598/1‑126 (MQ=255)
                                                                                                             cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTTTGTCACACt                                                                                                                                                                                                                         >  1:386050/1‑129 (MQ=255)
                                                                                                                  ttGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCacac                                                                                                                                                                                                                          >  1:150324/1‑123 (MQ=255)
                                                                                                                   tGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGGGGGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTAAAAAAGCCAAGGGGGCg                                                                                                                                                                                                                                          <  2:476368/106‑1 (MQ=255)
                                                                                                                             ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATg                                                                                                                                                                                                                                >  1:476368/1‑106 (MQ=255)
                                                                                                                                           tGAAGTTTCCCCCCAACGTGCGGGGGAACCAAATGCCGAAGTACCCGCGTTGCGTATGGGGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                                      <  1:1031567/102‑1 (MQ=255)
                                                                                                                                           tGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGGGTGGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                                      <  1:793876/102‑1 (MQ=255)
                                                                                                                                               ggTTCCCACCAACGTGGGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                                      <  1:2307094/98‑1 (MQ=255)
                                                                                                                                                gTTCCCCCCAACGGGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGc                                                                                                                                                                                                                                     <  1:1172885/82‑1 (MQ=255)
                                                                                                                                                 ttACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                                                                                                                                                                                       >  2:1914225/1‑95 (MQ=255)
                                                                                                                                                 ttACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                                                                                                                                                                                       <  1:1914225/95‑1 (MQ=255)
                                                                                                                                                   accaccAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTTGAGCCCTCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAACCGCAGGTCGCTGGGCGGGGACTTTACCGTCTGCAACAGGCTGGCATTGAc                                                                     >  1:959705/1‑239 (MQ=255)
                                                                                                                                                     ccccaACGGGCGGGTGAACCACTTGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGt                                                                   <  1:2362064/237‑1 (MQ=255)
                                                                                                                                                     caccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAAACAAAAGTGCGACCGCCCATGTCACACCCGAACCCTGTAGc                                                                                                                                                                                                            >  2:1031567/1‑102 (MQ=255)
                                                                                                                                                        caACGTGCGGGTGAACCACATGCCGAAGTCCCCGCGTTGCGTATGGCGGGTAAAAAACCCAAAGGTGCGACCGCCTATGTca                                                                                                                                                                                                                             <  2:308685/82‑1 (MQ=255)
                                                                                                                                                        caACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                               <  1:1058303/80‑1 (MQ=255)
                                                                                                                                                         aaCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGGGCGACCGCCCATGTCACACCCGCACCCTGTAGc                                                                                                                                                                                                            >  2:2307094/1‑98 (MQ=255)
                                                                                                                                                          aCGTGCGGGTGACCCACTTGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCaa                                                                                    <  1:1026852/217‑1 (MQ=255)
                                                                                                                                                          aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGGGCGACCGCCTATGTcaca                                                                                                                                                                                                                           >  2:1172885/1‑82 (MQ=255)
                                                                                                                                                          aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACCCGAACCCTGTAGcc                                                                                                                                                                                                           >  1:1291420/1‑98 (MQ=255)
                                                                                                                                                           cGTGGGGGTAAACCACATGCCGAATTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAGGGTGCAACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                             <  1:1859793/239‑1 (MQ=255)
                                                                                                                                                               cGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCCCCGCTGGCGGAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCCGGGCGGGGACTTTACCGTCTGCAAAAGGCTGGGCTTTACGCTGGCCACgg                                                          >  2:1659567/1‑238 (MQ=255)
                                                                                                                                                                  gTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCCCCTATGTCACACTCCAAcc                                                                                                                                                                                                                   >  1:308685/1‑82 (MQ=255)
                                                                                                                                                                  gTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCCATGTCACACTCGaa                                                                                                                                                                                                                     >  2:1058303/1‑80 (MQ=255)
                                                                                                                                                                   agcaaaacaTGCCGAAGTCCCCGCGTTGCGTAGGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                                         <  1:1567120/229‑1 (MQ=255)
                                                                                                                                                                    gAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGTGCGACCGCCCATGTCACACTCGAACCCTCTAGCCCTCATGGTTGTACGCCCCCGAGCTGTGACGCACCCATCGCCGCCGGGGGCGCGCGCGGGGGTGCCCCGGTGCAAACACCCAACCCCCCGGGCGCCGGGCGGGGACGTTTCCGTCGTCAAccagcttgcg                                                                          >  2:1026852/1‑208 (MQ=255)
                                                                                                                                                                          caTGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGGGACCCCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGGGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGGGGGTGCCTCGATGCAAGATCCTCACCCGCAGGTCGCCGGGGGTGGACTTTTCCGTCTGCAACAGGCTGGCATTGACGTCAGCCCCGGCCTGTTGATGAg                                              >  2:481685/1‑239 (MQ=255)
                                                                                                                                                                            tGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAACGGGCGACCGCCCTTTTCACACCCGAAACCCGTCGCCCCCATTGGCGCACGCCACCCTGCTGTGACGCACCCATCCCCCCTGGCGTAACGCGCGTGGTTTCCTCGATTCCCGATCCTAAACCCCAGGGCGCCGGGGGGGGGCTTTATCGTCTCCACCAAGCGGGCATTTGCGCACGCGCCGGCCtggtgatgcggg                                            >  2:261099/1‑235 (MQ=255)
                                                                                                                                                                                    tccccgcgTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                        <  1:1025152/167‑1 (MQ=255)
                                                                                                                                                                                    tacgggcgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCTCCGTGCTGTGACGCCCTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                         <  1:672381/133‑1 (MQ=255)
                                                                                                                                                                                    tACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                             <  1:2164049/166‑1 (MQ=255)
                                                                                                                                                                                    tACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTg                                                                                                                                                                                                                <  1:93245/67‑1 (MQ=255)
                                                                                                                                                                                    tACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGGGACCCCCCATGTCACACTCGAACCCCGTAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACTCATCGCCCCCGGCCGAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCCCAGGGCGCCGGGCGGGGGCTTTTCCGTCTGCAACAGGCCGGCCTTGACGCCAGTCCCGGCCcgattatgggctgttc                                        >  2:1567120/1‑227 (MQ=255)
                                                                                                                                                                                     acccgggtggCGTATGGGGGGTAAAAAAGCAAAGGGGCCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                                    <  1:1659567/229‑1 (MQ=255)
                                                                                                                                                                                             tGCGTATGGCGGTGAAAAAACCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                         <  2:2080484/129‑1 (MQ=255)
                                                                                                                                                                                              gCGTATGGCGGGTGGAAAAGACAAAGGTGAGACCGCCTATGTCCCACTCGGACCCCGGAGCCCCCATGGTCGTACGCCACCGTGCTGTGTCGCACTCTTCGCCGCTGGCGCAGCGCGCGGGGGTGCCTCGATGCCAGa                                                                                                                               >  2:672381/1‑138 (MQ=255)
                                                                                                                                                                                              gCGTATGGCGGGTGAAAAAGCCCAAGGTGCGACCGCCCATGTCACACTCGAACCCTGTAGGCAACAAGGGCGTACGCCACCGCGCTGTGACGCACTCATCCCCGCCGGCGTGGCGCGCGGGGTTGCCCCGCTGCAAGAACCTAACCCCCAGGTCCCAGGGGGGGGa                                                                                                   >  2:2164049/1‑166 (MQ=255)
                                                                                                                                                                                              gCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                         <  2:1654629/128‑1 (MQ=255)
                                                                                                                                                                                              gCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                         >  1:1654629/1‑128 (MQ=255)
                                                                                                                                                                                              gCGTATGGCGGGTGAAAAAGCCAAAAGTGCGACCCCCTATGTCACACTCCAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCCGGCGGAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACGTTACCGTCTGCAACAGGCTGGCATTGACGTCCGCCACGGACCGATGATGAGGGAAGCCGAGCAAATGAATa                          >  2:1470825/1‑239 (MQ=255)
                                                                                                                                                                                              gCGTATGGCGGGGGAAAAAACCAAAAGTGCGAACGCCCATGTCACACCCGAACCCCCTAGCCAACATGGTCGGACGCCCCCGGGCTGTTACGCACTCATCGCCGCCGGCGTAGCGCGCGTGGGTTCCTCGACGCAAGATCCTAACCCGCAGGTCGCTCGGGGGGGACTTTACCGCCTGCAACAGGCTGGGCCTGACGGGCGCCCCGGCCCGCGGGTGAGTGGGGACGAGGAAGTGaaaa                          >  2:1839344/1‑237 (MQ=255)
                                                                                                                                                                                              gCGTATGGCGGGGGAAAAAACCAAAAGGGCGACCCCCTATGTCCCACTCGAAACCTGTAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCCGGCGTAGCGCGCGGGGGTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCCGGGCGTGGGCTTTa                                                                                              >  2:1025152/1‑171 (MQ=255)
                                                                                                                                                                                               ggTAGGGGGGGTGAAAAAGCCAAAGGTGGGACCGCCTATGTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                                                       <  1:289998/64‑1 (MQ=255)
                                                                                                                                                                                               cgatgggCGGGTAAAAACGCAAAGGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATa                          <  2:191991/233‑1 (MQ=255)
                                                                                                                                                                                               cGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                                                                                                                                                                                       <  1:1020965/49‑1 (MQ=255)
                                                                                                                                                                                               cGTATGGCGGGTGAAAAAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACg                                                                                                                                                                                             <  2:1135937/75‑1 (MQ=255)
                                                                                                                                                                                                       gggTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCCTCGCCGCTGGCGTAGCGCGCGGGGTTGCCCCGATGCAAGa                                                                                                                               >  1:2080484/1‑129 (MQ=255)
                                                                                                                                                                                                        ggTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTtc                 <  1:1410252/238‑1 (MQ=255)
                                                                                                                                                                                                         gggAAAAAGCCAAAGGGGCGACCGCCCATGTCACACTCGAACCCTGTAg                                                                                                                                                                                                             >  2:1020965/3‑49 (MQ=255)
                                                                                                                                                                                                         gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                         <  1:2138917/117‑1 (MQ=255)
                                                                                                                                                                                                         gTGAAAAAGCCAAAGGTGCGACCCCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCt                                                                                                                                                                                   >  1:1135937/1‑75 (MQ=255)
                                                                                                                                                                                                         gTGAAAAAACCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                         >  2:2138917/1‑117 (MQ=255)
                                                                                                                                                                                                          taaaaaacaaaaaGGTCCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTtctc               <  2:1347674/228‑1 (MQ=255)
                                                                                                                                                                                                                   cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATgc     >  2:1963521/1‑239 (MQ=255)
                                                                                                                                                                                                                      tggTGCTACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCa  <  1:1325340/238‑1 (MQ=255)
                                                                                                                                                                                                                        gTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCa  <  1:1470825/237‑1 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
CGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTAC‑CACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCA  >  NC_000913/433433‑433884

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