Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6490TG46.4% 20.0 / 19.8 60G65G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (8/22);  new base G (17/9);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTC  >  NC_000913/433425‑433871
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              
gaggagaCCGGGCGCCGGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGGGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTCCCCCGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGTAAAAAACCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                  <  2:1094889/236‑1 (MQ=255)
 aggagaGCGCGCGCCTGGGGAACGAAGCCGTGCAGGAGGAGTATTACATGGGGCGGGCGCTAAAGCTGGCGCAACGGGGACGTTTTCCCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                  <  2:288539/236‑1 (MQ=255)
    ggacggcgcgcCGGGAGGACCAACCCGGGCAGGACGAGTATTACAGGGCGCGGGCGCAAAAGCGGGCGCAACGAGGAGTTTTTACCACGCATCCCAACCCGAATGGCGGGGGGGTAATTGTAAAAGATGGCAAAATTGTGGGTGAAGGTTACCCCCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGTAAAAAAGCAAAAGGTGCGACCGc                                                                                                                                                                                                                        <  2:1042659/224‑1 (MQ=255)
             gcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAAcc                                                                                                                                                                                                      <  1:927703/238‑1 (MQ=255)
              ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCCAAAGGGCGACCGCCTATGTcaca                                                                                                                                                                                                              >  1:1042659/1‑229 (MQ=255)
                               tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGACAAAGAGGTGAACCCCCATTGCACACCCGCACCCCGTAGCCCTCATgttggt                                                                                                                                                                                   >  1:1091920/1‑234 (MQ=255)
                               tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGt                                                                                                                                                                                   >  1:902372/1‑239 (MQ=255)
                                         gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCCATGTCACACTCGAAACCCGTAGCCATCATGGTCGTACGCCAc                                                                                                                                                                            >  1:1071804/1‑236 (MQ=255)
                                         gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGCGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCGCCTATGTCACACTCGAACCCCGTAGCCAACATGGGCGTACGCcccccg                                                                                                                                                                         >  1:18482/1‑237 (MQ=255)
                                          ttttACAGGGCGGGGGCGCAAAAGCGGGCGCACGGAGGTCGTTTTCCCACGCACCCCACCCCGAATGTGGGGTGCGTAATTGCAAAAGATGGCAAAATTGTCGGTGAAGGTTACCACCAGCGGGCGGGTAAACCACATGCCGAAGTCCCCGCGTTGCGTAGGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcat                                                                                                                                                                                         <  2:1000915/220‑1 (MQ=255)
                                           aTTACAGGGCGGGGGCGCTAAAGCGGGGGAAAGGAGGACTTTTTACCACGCACCCCACCCCGAATGTCGGGGGCGTCATTGTCAAAGATGGCAAAATTGTCGGGGAAGGTTACCACCAACGGGGGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                        <  2:221713/238‑1 (MQ=255)
                                                agggCGCGGGCACAAAAGCTGGCCCAACGAGGACTTTTTCCCACGCACCCCACCCCGAAGGTCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTAAACCACATGCCAAATTACACGCGTTGCGTAGGGGGGGTAAAAAAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcat                                                                                                                                                                                         <  2:70071/214‑1 (MQ=255)
                                                  ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAACCAAAGGGCGACCCCCCACATCTCACCCCCAACCCGCAGCCCTCAAGGTCGTACGCCACCggcccgttac                                                                                                                                                                 >  2:688885/1‑229 (MQ=255)
                                                  ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGACAAAGGGGGGACCGCCCAAGTCACCACCGAACCCCGTAg                                                                                                                                                                                                >  1:1477806/1‑207 (MQ=255)
                                                  ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGTGCGACCGCCCATGTCACAACCGAACCCTGTAGCCATCATGGGCGTACGCCACCGGGCTGTGAc                                                                                                                                                                 >  1:258871/1‑238 (MQ=255)
                                                    cgcgGGCGCTAAGGCGGGCGCAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  2:223653/214‑1 (MQ=255)
                                                    cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGGGCGACCGCCCATTTTACCCTCGAACCCCGGAGCCATCATGGGCGGAccc                                                                                                                                                                               >  1:1000915/1‑220 (MQ=255)
                                                    cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCCATCCTGGTCGTACGCCACCGTGCTGTGACgg                                                                                                                                                               >  1:1035009/1‑237 (MQ=255)
                                                    cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGCGGAAAAGACCAAAGAGCGACCGCCCATTTCACACTCGAACCCTGTAGCCACCATGGGCGCACTCCACCGCGCtctgtggc                                                                                                                                                               >  2:251123/1‑234 (MQ=255)
                                                          cgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCAAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                         <  1:907199/190‑1 (MQ=255)
                                                          cgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCTATTTCACACCCGAACCCCGTAGCCATCATGGGCGTACGc                                                                                                                                                                               >  1:70071/1‑216 (MQ=255)
                                                              aaaGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAGGTGGGACCGCCTATGTCACACTCGAAACCTGTAGCCACCATGGTCGGACGCCACCCTGCGGTGACGCCCTCaccacc                                                                                                                                                     >  1:1495544/1‑233 (MQ=255)
                                                              aaaGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCGCCCATGTCCCACTCGAACCCCGTAGGCCTCATGGTCGTAAGCCa                                                                                                                                                                             >  1:223653/1‑214 (MQ=255)
                                                                                    tttccccgccATCCCACCCGAAGTGTCGGGTGCTTCATTTTAAAGGTGGCCAAAATTGTCGGTGATGTTTCCCACCACCGGGGGGGTGACCCACATGCCGAAGTCCACGCGTTGCGTATGGCGGGTAAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGc                                                                                                                              <  2:1438577/230‑1 (MQ=255)
                                                                                     tttccaccgcccccAACCCGAAAGGTGGGGTGGTTCATGGTAAAAGGGGGCGAAATTGTCGGGGAGGGTTACCCCCAACGGGGGGGTAACCCAAAGGCCGAATGCCCCGCGTTGCGTAGGGGGGGTAAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGc                                                                                                                              <  2:1492409/228‑1 (MQ=255)
                                                                                             cATCCCAACCCGAAGGTGGGGTGCGTAATTGTAAAAGATGGCAAAATTGCCGGTGAAGGTTCCCACCACGGTGCGGGTGACCCACAGGCCGAAGTACCCGCGTTGCGTATGGCGGGTAAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  2:258871/239‑1 (MQ=255)
                                                                                               tCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  1:97617/171‑1 (MQ=255)
                                                                                                cccAACCCGAATGTCGGGTGCGTAATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgc                                                                                                                                          <  1:389440/215‑1 (MQ=255)
                                                                                                    aCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCt                                                                                                              <  1:1442660/239‑1 (MQ=255)
                                                                                                          aaTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGGGGCCGCCCTAGTCCCACTCGAAACCCGTAACCCTCATGGGCGTGAGCCa                                                                                                                                                                             >  2:97617/1‑171 (MQ=255)
                                                                                                           aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGCCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGGCAAAAGTGCGAACCCCCATGTCACACTCGAACCCTGTAGCCATCAATGGCGTACCCCCCCGTGCCGTGACGCACTCATCCCCCCCCGCGCGGCGCGCGCGGTTGCCCCGCTGCAAAATCCTAAccc                                                                                                         >  2:911937/1‑238 (MQ=255)
                                                                                                                            tgccAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCAc                                                                                                                                                             <  1:1140425/167‑1 (MQ=255)
                                                                                                                                aaaGGTGGCGAAATGGTGGGTGAAGGTTCCCACCGACGTGCGGGCGAACCACAGGCCGAAGTACCCGCGTGGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACt                                                                                    <  2:18482/239‑1 (MQ=255)
                                                                                                                                      ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCAAAATCGAACCCTGGAGCCATCATGGGCGGACGCCCCCGCGCTGTGGCGCCCTCCTcgcccc                                                                                                                                                   >  2:1140425/1‑168 (MQ=255)
                                                                                                                                         cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAAGGGCGACCGCCTCTGTCCCCCTCTAACCCCGTAGCCCTCAAGGGCGGTCGCCACCGTGCCTGGACGCACTCACTCCCGCCGGCGCAGCGCGCGTGGGTGCCCCCCTGCAAGAACCTCACCCCCCGGTCGCCggggcggggctct                                                                                  >  2:784248/1‑223 (MQ=255)
                                                                                                                                            aaTTGTCGGTGAAGGTTACCCCCAACGTGGGGGTAAACCACATGCCGAAGTCCACGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                 <  2:128558/119‑1 (MQ=255)
                                                                                                                                                  gggaagagtttACCCCCAGCGTGCGGGTGACCCACAGGCAAAAGTACACGCTTTGCTTAGGGGGGATAAAAAAGAAAAAGGCCGCACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAgg                                                                   <  2:1321206/230‑1 (MQ=255)
                                                                                                                                                      gAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCCAAGGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCCCTCATgg                                                                                                                                                                                       >  1:128558/1‑119 (MQ=255)
                                                                                                                                                          gTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTAGGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                 <  1:1305236/105‑1 (MQ=255)
                                                                                                                                                             cccaccaACGTGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  1:526287/111‑1 (MQ=255)
                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  2:1496431/105‑1 (MQ=255)
                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGCGGGACCGCCCATGTCCCACTCGGACCCCCGAGCCATCAAgg                                                                                                                                                                                       >  2:1305236/1‑105 (MQ=255)
                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       >  1:1496431/1‑105 (MQ=255)
                                                                                                                                                                     ggTGGGGGTCAACCATACGCGAAATTCCCCGGGTGGCGTAGGGGGGGAAAAAAAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                <  2:1275275/238‑1 (MQ=255)
                                                                                                                                                                     cGTGCGGGTGAACCACAGGCCAAAGTCCACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  2:1102462/170‑1 (MQ=255)
                                                                                                                                                                       tGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCCATCATGGGCGTACGccc                                                                                                                                                                             >  2:526287/1‑111 (MQ=255)
                                                                                                                                                                            gtaaaCCACATGCCAAAGTACACGCGTTGCGTAGGGCGGGTAAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  2:1392494/94‑1 (MQ=255)
                                                                                                                                                                               aaCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAc                                                                                                           >  1:1102462/1‑170 (MQ=255)
                                                                                                                                                                                  cacaTGCCGAAGTACGCGCTTTGCGTAGGGCGGGGAAAAAAGCAAAAGGGGCGACCGCCTATGTCACACGCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg                                      <  1:713147/236‑1 (MQ=255)
                                                                                                                                                                                  cacaTGCAGAAGTACGCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  1:1202867/157‑1 (MQ=255)
                                                                                                                                                                                      tGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCCAAGGTGCGACCGCCCATTTCACACCCGAACCCTGTAGCCATTATGGGCGTACGCCa                                                                                                                                                                             >  1:1392494/1‑97 (MQ=255)
                                                                                                                                                                                                 acggtttGCGTAGGGGGGGTAAAAAGGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                    <  1:672496/234‑1 (MQ=255)
                                                                                                                                                                                                 acGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGa                                                                                                                                                                                                                            <  1:274492/39‑1 (MQ=255)
                                                                                                                                                                                                     gTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                        <  1:1082469/87‑1 (MQ=255)
                                                                                                                                                                                                     gTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                        >  2:1082469/1‑87 (MQ=255)
                                                                                                                                                                                                       tGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                            <  1:1228678/129‑1 (MQ=255)
                                                                                                                                                                                                        gCGTATGGCGGGTGAAAAAGCCCAAGGTGCGACCGCCTCTGTCACACCCCGACCCTGTAGCCCTCATGGTCGTACGCCCCCGTGCTGTGACCCCCCCCTCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGCTGCAAGATCCCAACCCGCAGGTCGCTGGGCGTGGGCTTTACCGTCCGCAACAGGCTGGCCTTGACGTCAGCCCCGGCCCGGTGTTTGGCGAAGCCGCGGAATTTAATa             >  2:979070/1‑239 (MQ=255)
                                                                                                                                                                                                        gCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCCCTGTCACACCCGAAACCTGTAGCCACCCTGGTCGTACGCCACCGGGCTGGGACGCACTCATCCCCCCTGGCGTAGCGCGCGTGGTTGCCTCGATGCCAAATCCCAACCCGCAGGTCGCTGGGCGGGGGCTTTACCGTCGCCAACAGGCTGGCAATGACGCCGCCCACGGCCTGCTGTTGAGTGAAACCGAGCGATTgcaaa             >  2:843549/1‑235 (MQ=255)
                                                                                                                                                                                                                gggggTAAAAACGCCAAGGGGGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGCACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCt       <  2:1495544/236‑1 (MQ=255)
                                                                                                                                                                                                                    gggAAAAAACCAAAGGGGCGACCCCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGGCGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGGTGCAAGATCCTAACCCGCCGGTCGCCGGGCGGGGGCTTTACCGTCCGCCACACGCTGGCATTGACGTCAGCCACCGCCCGCTGATTGGGGAAGCCGGGCAATTGAATAAAGGCTTtttc  >  2:1516994/3‑239 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              
GAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTC  >  NC_000913/433425‑433871

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