Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6490TG35.1% 36.1 / 10.5 57G65G (GGT→GGG) ribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (10/27);  new base G (12/8);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAG  >  NC_000913/433425‑433874
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
aaagaacgcgcgcgCTGGGAGGACTAAGCCGTGAAGAGCGAGTTTT‑CCATGCCGCGGGCGAA‑AAAGCTGTCGCAACGAGGACTTTTTACCCCGCATCCCACCCCAAATGCCGGGGGCG‑TAATGGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTCCCACCACCGGGCGGGTGACCCACAGGCGGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAGGGT‑GCGa                                                                                                                                                                                                                                <  2:1757200/221‑1 (MQ=255)
              ccTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGT‑GAGACCGCCt                                                                                                                                                                                                                          >  1:1931443/1‑221 (MQ=255)
                  gaggcttacgCCGTGCAGGACG‑GTATTGACATGGCCGGGGCGCT‑AAAGCGGGCCAAACGAGGACGTTTTACCACGCATCCAAACCGGATGGGCGGGTGGG‑TCATTGTAAAAGATGGCAAAATGGTC‑GTGGAAGTTTCCCCCCAACGTGCGGGTGAACCAAAGGCCGAAGTACCCGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCg                                                                                                                                                                                                              <  2:1835612/220‑1 (MQ=255)
                   gggccTAAGCCGTGCAGGAGGAGTTTT‑ACAGGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGATGTTTTACCACGCATCCAACCCCGAATGTCGGGTGCG‑TCTTTGTCAAAGAGGGCAAAATTGTC‑GGTGAGGGTTCCCACCAACGGGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTAGGGGGGGTAAAAAAGCCAAGGGT‑GCGACCGCCTATGTCACACTCg                                                                                                                                                                                                              <  2:1890195/224‑1 (MQ=255)
                    ggaCTACGCGGGGCAGGAAGATTATT‑CAATGGCGCGGGCGAT‑AAAGCTGGCGCAACGAGGGATTTTTCCCCCGCACCCCAGCCGGAATGTCGGAGGCG‑TCATTGTCAAAGGTGGC‑AAAATGTCGGGTGAAGGTTCCCACCAGCGTGGGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGGGGGTAAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                                            <  2:894341/229‑1 (MQ=255)
                         aaGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGG‑GCGGCCGCCCCTGTCCCACTCCAACCCTg                                                                                                                                                                                                       >  1:1022155/1‑229 (MQ=255)
                            ccgagcggggagAGTTTT‑ACA‑GGGGCGGGCCCTAAAAGCGGGCGAAACGGGCACGTTTTACCACGCACCCAAACCCAAATGCGGGGTGCT‑TCATTGTCAAAGATGGCGAAATGGTG‑GGTAAAGTTTACCACCACGGTGCGGGTGAACCAAAGGCCGAAGTACCCGCGTTGCGTAGGGCGGGTAAAAAACCAAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                     <  2:634126/217‑1 (MQ=255)
                               tGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAACAAAGAT‑GCGa                                                                                                                                                                                                                                >  1:654441/1‑198 (MQ=255)
                               ggCAGGACGAGTTTG‑ACAGGGCGGGGGCGGT‑AAAGCGGGCGAAACGAGGAGGTTTTCCCCCGCATCCCAACCCTAATGTGGGGTGCG‑TAATTGTAAAAGATGGAGAAATTGCC‑GGTGAGGGTTACCCCCAAGGTGGGGTTAAACCAAATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGT‑GCgg                                                                                                                                                                                                                                <  2:653289/197‑2 (MQ=255)
                                  agaaCGTGTATT‑CCATGGCGCGGTCGCT‑AAAGTTGGCGCAAGAAGGACGTTTTCCCACGCATCCCAACCCGATTGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGTTTACCACCACCGTGGGGGTGAACCACAGGCCGAAGTCCACGCGTGGCTTATGGCGGGTAAAAAACCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                                              <  2:1196728/226‑1 (MQ=255)
                                       ttgTATC‑ACAGGGCGGGGGGGCA‑AAAGCGGGCGCACGGAGGCCGTTTTCCCACGCCTCCAACCCAGAATGCCGGGTGCG‑TCATGTGCAAAGAGGGCGAAATTGTC‑GGGGAGGGTTACCACCACCGGGGGGGTGAACCACAGGCGGAGGTACCCGCGTTGCGTAGGGGGGGTAAAAAAGCCAAGGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                                                         <  2:726838/227‑1 (MQ=255)
                                         gTATT‑CCAGGGCGCGGGCGCA‑AAGGGTGGCCGAACGAGGCTGTTTTCCCACGCACCCCACCCCGAATGTCGGGGGCG‑TCATTGTAAAAGAGGGCGAAATTGTC‑GGTGAAGGTTCCCACCACCGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTAAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAAcc                                                                                                                                                                                                          <  2:1059573/210‑1 (MQ=255)
                                         gTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGG‑GCGACCCCCTAGTGCCCACTCCAAcc                                                                                                                                                                                                          >  1:1061441/1‑210 (MQ=255)
                                           tgtcaCAGGGCGCGGACGAT‑AAGGTGGGCGCAACGAGGCATTTTTCCCCCCCCCCCCAGCCCGAATG‑CGGGGGCGTTCAGTGGAAAAGAGGGGGAAATTGTG‑GGTAAAGGTTCCCCCCAGGGGGCGGGAGAACCCAATGCCGAAGTACAGGCGTTGCGTAGGGCGGGTAAAAAACCCAAGGGT‑GCAACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                           <  2:5082/220‑1 (MQ=255)
                                                    gcgcGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGCGGAAAAGCCCAAGGT‑GCGGCCCGCCATGTCACTCTCTGACCCCGTCACCACCATAGTccgtc                                                                                                                                                                                     >  2:549868/1‑217 (MQ=255)
                                                     cgcgGGGGCT‑AAGGCGGGGGCAAGGGGGACGTTTTACCACGCATCCCACCCCGAATGTGGGGTGCG‑TAATTGTAAAAGATGGCAAAATTGTC‑GGTGAAGGTTCCCACCAACGGGCGGGTGAACCACATGCGGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGT‑GCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                             <  2:372333/196‑1 (MQ=255)
                                                     cgcgGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGG‑GCGACCCCCCATGTCacacccca                                                                                                                                                                                                             >  1:373006/1‑192 (MQ=255)
                                                                                                  ccaaccccGAAGGGCGGGTGCG‑TCATTGTCAAAAGGGGCGAAATTGGC‑GGTGAGGGTTACCACCAACGTGCGGGTAACCCACATGCCGAAGAACACGCGTTGCGTAGGGCGGGTAAAAAAGCAAAGGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCt                                                                                                                               <  2:690712/224‑1 (MQ=255)
                                                                                                           aaggTCGGGTGCG‑TCATTGTAAAAGAGGGCGAAATTGTC‑GGTGAAGTTTACCACCACCGTGCGGGTAAACCACATGCCGAAGTCCCCGCTTTGCGTAGGGGGGGGGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCg                                                                                                                                                                                                              <  2:1797133/139‑1 (MQ=255)
                                                                                                           aaTGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGG‑GCGACCCCCCATTTCACAATCg                                                                                                                                                                                                              >  1:1800190/1‑142 (MQ=255)
                                                                                                                 gggTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACtt                                                                                                                                                                                                               >  1:689956/1‑134 (MQ=255)
                                                                                                                 gggTGCG‑TCATTGTAAAAGATGGGGAAATTGTC‑GGGGAAGGTTACCACCAAGGTGGGGGTGACCCACATGCCGAAGTACACGGGTTGCGTATGGCGGGTGAAAAAGCCAAGGGT‑GCGACCGCCTATGTCACACTc                                                                                                                                                                                                               <  2:688745/135‑1 (MQ=255)
                                                                                                                  ggTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCCAAGGG‑GCGCCCCCCTATGTCCCCATCGAACCCTTGAGCCATCATTGCCGTCCGCCCCCGTGCTGTGACGCACCCATCGCCGCTGGCGGAGCGCGCGCGGGTG‑CCTCCATGCAAGACCCTaa                                                                                                               >  2:1254233/1‑229 (MQ=255)
                                                                                                                      cG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGG‑GGGAACCCCCATGTCACACTCGGACCCCGTa                                                                                                                                                                                                     >  1:40997/1‑140 (MQ=255)
                                                                                                                      cG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCCAAGGT‑GCGAACCCCTATGACACACTCGGACCCTGGAgaccccat                                                                                                                                                                                             >  1:1236294/1‑141 (MQ=255)
                                                                                                                      cG‑TCATTGTAAAAGATGGCAAATTGGTC‑GGGGAAGGTTACCACCAACGTGGGGGTGAACCACATGCCGAAGTCCCCGTGTTGCGTATGGCGGGGAAAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcat                                                                                                                                                                                             <  2:1234145/148‑1 (MQ=255)
                                                                                                                         tCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCa                                                                                                         <  1:134101/229‑1 (MQ=255)
                                                                                                                                    ggtGGCGAAATTGGC‑GGTAAAGGTTCCCCCCAACGTGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                            <  1:778339/150‑1 (MQ=255)
                                                                                                                                       ggCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                           >  1:393296/1‑134 (MQ=255)
                                                                                                                                       ggCAAAATTGTC‑GGTGAGGTTTCCCACCAGCGGGGGGGTAAACCACAGGCCAAATTCCACGCGTTGCGTTGGGGGGGTGAAAAAGCCAAAGGT‑GCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                           <  2:392585/134‑1 (MQ=255)
                                                                                                                                         cGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGACCAAGGG‑GGGAACGCCCATGTCACACCCGAACCCCGTAGCCACCCTGGCCGGACCCCACCGGGCTGTGACGCACCCATCGCCCCTGGCGGAGCGCCCGCGGGTG‑CCCCGAGTGAAGATCCCAACCCCCCGGTCgacgtgggtgga                                                                                         >  2:1365718/1‑218 (MQ=255)
                                                                                                                                         cGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCAAGGTCGGACGCCACCCGGCTGTGACg                                                                                                                                                                    >  1:405879/1‑155 (MQ=255)
                                                                                                                                         agAAATTGTC‑GGTGAAGTTTCCCACCACCGTGCGGGTAAACCACAGGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACg                                                                                                                                                                    <  2:405147/154‑1 (MQ=255)
                                                                                                                                                                     aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGT‑GCGACCGCCCATGGCACACTCGAACCCTGTAGCCATCATGGGCGTACGccc                                                                                                                                                                                 >  1:1274973/1‑114 (MQ=255)
                                                                                                                                                                     aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGG‑GCGACCGCCCATGTCACACCCGAACCCCGTAGCCCTCATGGTCGGACGCCCCCCTGCCGTTACCCACCCACCGCCCCTGGCGTAGCCCGCCGGGTTGCCCGCGTTCAAAGTCCCCAACCC‑CCGGGGGCTGGGCGCGGGCTTTACCGCCCGCAAcctgcgtgcta                                                               >  2:1616900/1‑218 (MQ=255)
                                                                                                                                                                     aCGGGCGGGTAACCCACATCCCGAAGTACACGCGTTGCTTAGGGCGGGTAAAAAACCAAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                 <  2:1272757/115‑1 (MQ=255)
                                                                                                                                                                                acccACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGT‑GCACCCGCCTATGTCACACTCGAATCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                   <  1:1727341/227‑1 (MQ=255)
                                                                                                                                                                                aaCCACGTGCGGAAGTACACGCGTGGCGTATGGCGGGTGAAAAGGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                   <  1:1475723/229‑1 (MQ=255)
                                                                                                                                                                                    ccaTGCCGAAGTCCCCGCGTGGCGTAGGGGGGGGTAAAAACCCAAGGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                   <  1:1256415/224‑1 (MQ=255)
                                                                                                                                                                                    acaTGCCGAGGTACCCGCGTTGCTTAGGGCGGGTGAAAAAGCAAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                   <  1:40124/225‑1 (MQ=255)
                                                                                                                                                                                               tACCCGCTTTGCGTATGGGGGTTAAAAAACCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCt                                                                                                                               <  1:1030688/138‑1 (MQ=255)
                                                                                                                                                                                                       ttGCGTTGGGCGGGTAAAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                                 <  1:630278/208‑1 (MQ=255)
                                                                                                                                                                                                       ttGCGTATGGCGGGTGAAAAAGCCAAAGGT‑GCGACCGcc                                                                                                                                                                                                                           >  1:615564/1‑39 (MQ=255)
                                                                                                                                                                                                       ttGCGTATGGCGGGTGAAAAAGCCAAAGGT‑GCGACCGcc                                                                                                                                                                                                                           <  2:614486/39‑1 (MQ=255)
                                                                                                                                                                                                         gCGTATGGCGGGTGAAAAAGCCAAAGGG‑GCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                                            >  1:379373/1‑52 (MQ=255)
                                                                                                                                                                                                         gCGTATGGCGGGGGAAAAAACCCAAAGG‑GCGACCGCCTATGTCACACCCGAACCCCGTAACCATCAAGGTCGTACGCCACCGTGCTGTGACACACACATCCCCCCTGGCGTAGCGCGAGTGGGTG‑CCCCGGTGCAAGAACCTAACCCACACGGCGCCGGGCGTGGGGTTTTGCGTCTGGAACAGGCTGGGCTTGACCTCAGCCACCGCCTGGTTGTGAGTGAAGCCgcg                          >  2:698928/1‑227 (MQ=255)
                                                                                                                                                                                                         gCGTATGGCGGGGGAAAAAACCAAAAGT‑GCGACCCCCCATGTCACCCCCGAACCCCGTAGCCATCATGGTCGTACGCCACCCTGCTGTGGCGCACCCCTCGCCGCCGGCGTAGCGCGCCTGGTTG‑CCTCGATGCAAGATCCTAACCCCCAGGGCGCTGGGCGTGGGCTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGGGGAGCCgcg                          >  2:438827/1‑227 (MQ=255)
                                                                                                                                                                                                          cGTATGGGGGGGGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                          <  1:439610/228‑1 (MQ=255)
                                                                                                                                                                                                           cttagggcggTTAAAAAGCCCAAGGTGGCAGCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                          <  1:1725523/220‑1 (MQ=255)
                                                                                                                                                                                                             aTGGCGGGAGAAAGGGCCAAGGGT‑CCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAt                      <  1:1101205/229‑1 (MQ=255)
                                                                                                                                                                                                                ggggggGAAAAAGCAAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCa                                                                                                                        <  1:335093/126‑1 (MQ=255)
                                                                                                                                                                                                                    ttaaagaggCCAAGGGT‑GCGGCCGCCTATGTCACGCTCGACCCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATa                <  1:1256466/220‑1 (MQ=255)
                                                                                                                                                                                                                    gTGAAAAAGCCAAGGGG‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                            <  1:514265/152‑1 (MQ=255)
                                                                                                                                                                                                                    gTGAAAAAGCCAAAGGG‑GCGACCCCCTATTTCCCACCCCAACCCCGTAGCCATCATGGTCGTACGCCACCCTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCCCGATGCAAAATCCTAACCCCCAGGTCGCTGGGCGGGGGCCTTTCCCTCTGCAACAGGCTGGCATTGACGTCAGCCCACGCCCGAAGATGAATGAAACCGAGCAATTGaaaa                >  2:905191/1‑226 (MQ=255)
                                                                                                                                                                                                                    gTGAAAAAGCAAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCaa                                                                          <  1:311635/170‑1 (MQ=255)
                                                                                                                                                                                                                           aGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCCTGGGCGTACGCCACCGGGCTGTGACGCACTCATCCCCGCCGGCGTAGCGCGCGTGGGTG‑CCTCGATGGAAGGTCCTAACCCGCCGGGCCGCGGGCGGGGGCCTTTCCGTCCGCAAAAGGGCGGCATTGACGGCCGCCACGGCCTGGTGGTGGGGGGAGCCGCGGCATTGGATTAAgggctt        >  2:187619/1‑225 (MQ=255)
                                                                                                                                                                                                                                 gggT‑GGGACCGCCTATGTCCCACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg  <  1:906770/228‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
GAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATT‑ACATGGCGCGGGCGCT‑AAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCG‑TCATTGTCAAAGATGGCGAAATTGTC‑GGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGT‑GCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTG‑CCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAG  >  NC_000913/433425‑433874

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