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→GGG) ribDfused 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.

GAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTC  >  NC_000913/433425‑433871
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              
gaggagaCCGGGCGCCGGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGGGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTCCCCCGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGAGGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTAAAAAACCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                  <  2:1094889/236‑1 (MQ=255)
 aggagaGCGCGCGCCTGGGGAACGAAGCCGTGCAGGAGGAGTATTACATGGGGCGGGCGCTAAAGCTGGCGCAACGGGGACGTTTTCCCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                  <  2:288539/236‑1 (MQ=255)
    ggacggcgcgcCGGGAGGACCAACCCGGGCAGGACGAGTATTACAGGGCGCGGGCGCAAAAGCGGGCGCAACGAGGAGTTTTTACCACGCATCCCAACCCG‑AA‑TGGCGGGGGGGTAATTGTAAAAGATGG‑CAAAATTGTGGGTGAAGGTTACCCCCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGC‑GGGTAAAAAAGCAAAAGGTGCGACCGc                                                                                                                                                                                                                        <  2:1042659/224‑1 (MQ=255)
             gcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAAcc                                                                                                                                                                                                      <  1:927703/238‑1 (MQ=255)
              ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAACCCAAAGGGCGACCGCCTATGTcaca                                                                                                                                                                                                              >  1:1042659/1‑229 (MQ=255)
                               tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGACAAAGAGGTGAACCCCCATTGCACACCCGCACCCCGTAGCCCTCATgttggt                                                                                                                                                                                   >  1:1091920/1‑234 (MQ=255)
                               tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAACCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGt                                                                                                                                                                                   >  1:902372/1‑239 (MQ=255)
                                         gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGGGCGACCGCCCATGTCACACTCGAAACCCGTAGCCATCATGGTCGTACGCCAc                                                                                                                                                                            >  1:1071804/1‑236 (MQ=255)
                                         gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGCGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAAGGGCGACCGCCTATGTCACACTCGAACCCCGTAGCCAACATGGGCGTACGCcccccg                                                                                                                                                                         >  1:18482/1‑237 (MQ=255)
                                          ttttACAGGGCGGGGGCGCAAAAGCGGGCGCACGGAGGTCGTTTTCCCACGCACCCCACCCCG‑AA‑TGTGGGGTGCGTAATTGCAAAAGATGG‑CAAAATTGTCGGTGAAGGTTACCACCAGCGGGCGGGTAAACCACATGCCGAAGTCCCCGCGTTGCGTAGGGG‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcat                                                                                                                                                                                         <  2:1000915/220‑1 (MQ=255)
                                           aTTACAGGGCGGGGGCGCTAAAGCGGGGGAAAGGAGGACTTTTTACCACGCACCCCACCCCG‑AA‑TGTCGGGGGCGTCATTGTCAAAGATGG‑CAAAATTGTCGGGGAAGGTTACCACCAACGGGGGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                        <  2:221713/238‑1 (MQ=255)
                                                agggCGCGGGCACAAAAGCTGGCCCAACGAGGACTTTTTCCCACGCACCCCACCCCG‑AA‑GGTCGGGTGCGTCATTGTAAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTAAACCACATGCCAAATTACACGCGTTGCGTAGGGG‑GGGTAAAAAAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcat                                                                                                                                                                                         <  2:70071/214‑1 (MQ=255)
                                                  ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAACCAAAGG‑GCGACCCCCCACATCTCACCCCCAACCCGCAGCCCTCAAGGTCGTACGCCACCggcccgttac                                                                                                                                                                 >  2:688885/1‑229 (MQ=255)
                                                  ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGACAAAGGGGGGACCGCCCAAGTCACCACCGAACCCCGTAg                                                                                                                                                                                                >  1:1477806/1‑207 (MQ=255)
                                                  ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAACCAAAGGTGCGACCGCCCATGTCACAACCGAACCCTGTAGCCATCATGGGCGTACGCCACCGGGCTGTGAc                                                                                                                                                                 >  1:258871/1‑238 (MQ=255)
                                                    cgcgGGCGCTAAGGCGGGCGCAACGAGGACGTTTTACCACGCATCCCACCCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  2:223653/214‑1 (MQ=255)
                                                    cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAACCAAAAGGGCGACCGCCCATTTTACCCTCGAACCCCGGAGCCATCATGGGCGGAccc                                                                                                                                                                               >  1:1000915/1‑220 (MQ=255)
                                                    cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAAGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCCATCCTGGTCGTACGCCACCGTGCTGTGACgg                                                                                                                                                               >  1:1035009/1‑237 (MQ=255)
                                                    cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGCGGAAAAGACCAAAGAGCGACCGCCCATTTCACACTCGAACCCTGTAGCCACCATGGGCGCACTCCACCGCGCtctgtggc                                                                                                                                                               >  2:251123/1‑234 (MQ=255)
                                                          cgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCAAAGTACACGCGTTGCGTATGGC‑GGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                         <  1:907199/190‑1 (MQ=255)
                                                          cgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAACCAAAGGGGCGACCGCCTATTTCACACCCGAACCCCGTAGCCATCATGGGCGTACGc                                                                                                                                                                               >  1:70071/1‑216 (MQ=255)
                                                              aaaGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAACCAAAGGTGGGACCGCCTATGTCACACTCGAAACCTGTAGCCACCATGGTCGGACGCCACCCTGCGGTGACGCCCTCaccacc                                                                                                                                                     >  1:1495544/1‑233 (MQ=255)
                                                              aaaGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAACCAAAAGGGCGACCGCCCATGTCCCACTCGAACCCCGTAGGCCTCATGGTCGTAAGCCa                                                                                                                                                                             >  1:223653/1‑214 (MQ=255)
                                                                                    tttccccgccATCCC‑ACCCG‑AAGTGTCGGGTGCTTCATT‑TTAAAGGTGGCCAAAATTGTCGGTGATGTTTCCCACCACCGGGGGGGTGACCCACATGCCGAAGTCCACGCGTTGCGTATGGC‑GGGTAAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGc                                                                                                                              <  2:1438577/230‑1 (MQ=255)
                                                                                     tttccaccgcccccAACCCGAAA‑GGTGGGGTGGTTCATGGTAAAAGGGGG‑CGAAATTGTCGGGGAGGGTTACCCCCAACGGGGGGGTAACCCAAAGGCCGAATGCCCCGCGTTGCGTAGGGG‑GGGTAAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGc                                                                                                                              <  2:1492409/228‑1 (MQ=255)
                                                                                             cATCCCAACCCG‑AA‑GGTGGGGTGCGTAATTGTAAAAGATGG‑CAAAATTGCCGGTGAAGGTTCCCACCACGGTGCGGGTGACCCACAGGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTAAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  2:258871/239‑1 (MQ=255)
                                                                                               tCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  1:97617/171‑1 (MQ=255)
                                                                                                cccAACCCG‑AA‑TGTCGGGTGCGTAATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgc                                                                                                                                          <  1:389440/215‑1 (MQ=255)
                                                                                                    aCCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCt                                                                                                              <  1:1442660/239‑1 (MQ=255)
                                                                                                          aa‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAACCAAAAGGGGGGCCGCCCTAGTCCCACTCGAAACCCGTAACCCTCATGGGCGTGAGCCa                                                                                                                                                                             >  2:97617/1‑171 (MQ=255)
                                                                                                           a‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGCCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGGCAAAAGTGCGAACCCCCATGTCACACTCGAACCCTGTAGCCATCAATGGCGTACCCCCCCGTGCCGTGACGCACTCATCCCCCCCCGCGCGGCGCGCGCGGTTGCCCCGCTGCAAAATCCTAAccc                                                                                                         >  2:911937/1‑238 (MQ=255)
                                                                                                                            tgccAAAGATGG‑CGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCAc                                                                                                                                                             <  1:1140425/167‑1 (MQ=255)
                                                                                                                                aaaGGTGG‑CGAAATGGTGGGTGAAGGTTCCCACCGACGTGCGGGCGAACCACAGGCCGAAGTACCCGCGTGGCGTATGGG‑GGGTGAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACt                                                                                    <  2:18482/239‑1 (MQ=255)
                                                                                                                                      gg‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCAAAATCGAACCCTGGAGCCATCATGGGCGGACGCCCCCGCGCTGTGGCGCCCTCCTcgcccc                                                                                                                                                   >  2:1140425/1‑168 (MQ=255)
                                                                                                                                         cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAAGGGCGACCGCCTCTGTCCCCCTCTAACCCCGTAGCCCTCAAGGGCGGTCGCCACCGTGCCTGGACGCACTCACTCCCGCCGGCGCAGCGCGCGTGGGTGCCCCCCTGCAAGAACCTCACCCCCCGGTCGCCggggcggggctct                                                                                  >  2:784248/1‑223 (MQ=255)
                                                                                                                                            aaTTGTCGGTGAAGGTTACCCCCAACGTGGGGGTAAACCACATGCCGAAGTCCACGCGTTGCGTATGGC‑GGGTGAAAAAGCAAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                 <  2:128558/119‑1 (MQ=255)
                                                                                                                                                  gggaagagtttACCCCCAGCGTGCGGGTGACCCACAGGCAAAAGTACACGCTTTGCTTAGGGG‑GGATAAAAAAGAAAAAGGCCGCACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAgg                                                                   <  2:1321206/230‑1 (MQ=255)
                                                                                                                                                      gAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCCAAGGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCCCTCATgg                                                                                                                                                                                       >  1:128558/1‑119 (MQ=255)
                                                                                                                                                          gTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTAGGGC‑GGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                 <  1:1305236/105‑1 (MQ=255)
                                                                                                                                                             cccaccaACGTGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  1:526287/111‑1 (MQ=255)
                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  2:1496431/105‑1 (MQ=255)
                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAACCAAAAGCGGGACCGCCCATGTCCCACTCGGACCCCCGAGCCATCAAgg                                                                                                                                                                                       >  2:1305236/1‑105 (MQ=255)
                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       >  1:1496431/1‑105 (MQ=255)
                                                                                                                                                                     ggTGGGGGTCAACCATACGCGAAATTCCCCGGGTGGCGTAGGGG‑GGGAAAAAAAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                <  2:1275275/238‑1 (MQ=255)
                                                                                                                                                                     cGTGCGGGTGAACCACAGGCCAAAGTCCACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  2:1102462/170‑1 (MQ=255)
                                                                                                                                                                       tGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCAAAAGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCCATCATGGGCGTACGccc                                                                                                                                                                             >  2:526287/1‑111 (MQ=255)
                                                                                                                                                                            gtaaaCCACATGCCAAAGTACACGCGTTGCGTAGGGC‑GGGTAAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                       <  2:1392494/94‑1 (MQ=255)
                                                                                                                                                                               aaCCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAc                                                                                                           >  1:1102462/1‑170 (MQ=255)
                                                                                                                                                                                  cacaTGCCGAAGTACGCGCTTTGCGTAGGGC‑GGGGAAAAAAGCAAAAGGGGCGACCGCCTATGTCACACGCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg                                      <  1:713147/236‑1 (MQ=255)
                                                                                                                                                                                  cacaTGCAGAAGTACGCGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  1:1202867/157‑1 (MQ=255)
                                                                                                                                                                                      tGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAGCCCAAGGTGCGACCGCCCATTTCACACCCGAACCCTGTAGCCATTATGGGCGTACGCCa                                                                                                                                                                             >  1:1392494/1‑97 (MQ=255)
                                                                                                                                                                                                 acggtttGCGTAGGGG‑GGGTAAAAAGGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                    <  1:672496/234‑1 (MQ=255)
                                                                                                                                                                                                 acGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGa                                                                                                                                                                                                                            <  1:274492/39‑1 (MQ=255)
                                                                                                                                                                                                     gTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                        <  1:1082469/87‑1 (MQ=255)
                                                                                                                                                                                                     gTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                        >  2:1082469/1‑87 (MQ=255)
                                                                                                                                                                                                       tGCGTATGGC‑GGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                            <  1:1228678/129‑1 (MQ=255)
                                                                                                                                                                                                        gCGTATGGC‑GGGTGAAAAAGCCCAAGGTGCGACCGCCTCTGTCACACCCCGACCCTGTAGCCCTCATGGTCGTACGCCCCCGTGCTGTGACCCCCCCCTCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGCTGCAAGATCCCAACCCGCAGGTCGCTGGGCGTGGGCTTTACCGTCCGCAACAGGCTGGCCTTGACGTCAGCCCCGGCCCGGTGTTTGGCGAAGCCGCGGAATTTAATa             >  2:979070/1‑239 (MQ=255)
                                                                                                                                                                                                        gCGTATGGC‑GGGGGAAAAAACCAAAGGGGCGACCGCCCCTGTCACACCCGAAACCTGTAGCCACCCTGGTCGTACGCCACCGGGCTGGGACGCACTCATCCCCCCTGGCGTAGCGCGCGTGGTTGCCTCGATGCCAAATCCCAACCCGCAGGTCGCTGGGCGGGGGCTTTACCGTCGCCAACAGGCTGGCAATGACGCCGCCCACGGCCTGCTGTTGAGTGAAACCGAGCGATTgcaaa             >  2:843549/1‑235 (MQ=255)
                                                                                                                                                                                                                gggggTAAAAACGCCAAGGGGGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGCACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCt       <  2:1495544/236‑1 (MQ=255)
                                                                                                                                                                                                                    gggAAAAAACCAAAGGGGCGACCCCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGGCGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGGTGCAAGATCCTAACCCGCCGGTCGCCGGGCGGGGGCTTTACCGTCCGCCACACGCTGGCATTGACGTCAGCCACCGCCCGCTGATTGGGGAAGCCGGGCAATTGAATAAAGGCTTtttc  >  2:1516994/3‑239 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              
GAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCG‑AA‑TGTCGGGTGCGTCATTGTCAAAGATGG‑CGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTC  >  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: ‑