Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,465,7870GA30.8% 55.4 / 10.1 39R132C (CGC→TGC) treBtrehalose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base G (19/8);  new base A (0/12);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 4.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCC  >  NC_000913/4465537‑4465997
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                   
acTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCg                                                                                                                                                                                                                     >  1:102524/1‑251 (MQ=255)
     cGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAAc                                                                                                                                                                                                                <  2:153777/251‑1 (MQ=255)
      ggCAGCTGCTGCCGGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCCTTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGGTCAACCACAGAAAATCGTAGCTCGTTTGCAGGGACGGGTCCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCCCCTTGGGAAAAcc                                                                                                                                                                                                               >  2:62989/1‑251 (MQ=255)
       gCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCg                                                                                                                                                                                                              <  2:295521/251‑1 (MQ=255)
           gtgctgcACGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGGGCCGCCCATGTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAAACACAGAAAATCGTAGAGCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTATGACCGTGGCTCAAGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGa                                                                                                                                                                                                         <  2:188043/250‑1 (MQ=255)
             gctgcCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATc                                                                                                                                                                                                       >  1:152848/1‑251 (MQ=255)
                     aGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAAccgcc                                                                                                                                                                                               >  2:117970/1‑251 (MQ=255)
                      gCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCTTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAAGCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAAccgccg                                                                                                                                                                                              <  1:263093/251‑1 (MQ=255)
                         gATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCt                                                                                                                                                                                            >  2:14465/1‑250 (MQ=255)
                               cGTTCATCAGCTGTGGAGAAACCAGTGTCACGACAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCt                                                                                                                                                                                            >  1:280432/1‑244 (MQ=255)
                                  tcatcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAAcgcg                                                                                                                                                                                  >  1:346023/1‑251 (MQ=255)
                                  tcatcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑GAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAAcgcg                                                                                                                                                                                  >  2:238390/1‑251 (MQ=255)
                                       ggttggggAGAAACCAGTGTCACGCAAAGCACGATACCAAGGATCGGCGTGCCCCCCATTTTTTC‑CACCGCTGACCAAAAAATACCGCCGGGCAGG‑AAGAAGAAGATCGCTTCGCCGATCAACCAAAAAAAAGGGTAGATCGTTTGCAGGGAAGGGTACATTGGCGCCGGCGTCTGGCCGTTGCTCAGGGGAAAA‑TCGCCGATCACATTGGAAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGgcag                                                                                                                                                                             <  1:214557/245‑1 (MQ=255)
                                               agaAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACCCGGGCAGCAACggg                                                                                                                                                                      >  2:250894/1‑250 (MQ=255)
                                                   aCCAGTGTCACGCCAAGCACAATACCAAGGATCGGCGT‑CCCCCCATTTTTTTTCCCCGCTGACAAGAAAATACCGCCCGGAAGG‑AAGAAGAAGACCGCTTCCCCTATCAACCACAAAAAATCGTAGATCGTTTGCAGGGAAGGGTACTTTGGCGCCAGCGTCTGCCCGTTGCTCATGGGAAAA‑GCGCCGATCACATTGAGAAAACCG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaagaa                                                                                                                                                                 <  1:333271/251‑1 (MQ=255)
                                                    ccAGGGTCACGCAAAGCAGGATCCCAAGGAGCGGCGTGCCGCCCTTTTTTTC‑CCCCGCTGACCAGAAAATACCGACCGGAAGG‑AAAAAAAAGACCCCTCCACCAATAAACAACAAAAAAGCGTAGATCGTTTGCAGGGAAGGGTACATTGGCGCCAGCGTCTGACCGTGGCACATGGGAAAA‑CCGCCAATCACATTG‑GAAAACCGGAGGAAAAATCCGCCGCTAATCAACGCGGGCAGCAACGGGaagaag                                                                                                                                                                <  1:161881/251‑1 (MQ=255)
                                                                     aCGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CAACGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATAGCTTCACCGATCAACCACAGCAAATCGTAGATCGTTTGCAGGGCAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATg                                                                                                                                                  >  1:206363/1‑248 (MQ=255)
                                                                         acccagggacctgctggccgccCTTTTTTTT‑CCCCGCACACCACAAAACCCCGCCCGGGAGA‑TAGAAAAAGAGCGCTTCCCCAACCACCCACAAAAAATGGTAGATCGTTTGCAGGAAAGGGAACATTTGCGCCAGCGTCGGCCTGTTGCCCAGGGGAAAA‑CCGCCAATCACATTGCAAAAACCG‑AGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATc                                                                                                                                           <  2:154243/235‑1 (MQ=255)
                                                                           ccAAGGAGCGGCTTGCCGCCCTTTTTTTC‑CACCGTTGCCCAGCAAATACCGACCGGCGGG‑AAGAAGAAGTTCCCTCCCCGGAACAACCAAAAAAAGTCGTAGTTTTTGTCCGGGGAGGGGTACATTTGCCCCAGCTCCTGCCCGTGGCTCAGGGGCAAA‑TCGCCACTCAAATGGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCaa                                                                                                                                         <  1:312297/251‑1 (MQ=255)
                                                                            cAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAAc                                                                                                                                        >  1:21238/1‑251 (MQ=255)
                                                                             agggATCGGCGCGCCGCCCATTTTTTC‑CCCCGCTGCCCACAAAATCCCGCCCGGCAGG‑TAGAAAAAGATCGCTTCACCGACCAACCACAAAAAATCGTAGATCGTTGGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCTTGGCTCAGGGGAAAA‑CCGCCGATCCAATTGCAAAACCCG‑AGGACCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACt                                                                                                                                       <  1:364984/249‑1 (MQ=255)
                                                                                 accGGCGTGCCGCCCATTTTTTC‑CCCCGCTGACCAGCAAATCCCGACCGGCAGG‑TAGAAGAAGATCGCTCCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGATACTTTTGCCCCAGCGTCTGCCCGTTGCTCA‑GGGCAAAATCCCCGATCACATTGCAAAACCCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTc                                                                                                                                   <  1:228450/249‑1 (MQ=255)
                                                                                     gCGTGCCCCCTATTTTTTC‑CACCGCGGACCAACAAAAACCGACGGGCAGG‑TAAAAGAAGATCGTTTAACCGATCACCCAAAGAAAATGGTAGATCGTTGGCAGGGAAGGGTACATTTGCGCCAGCGTCGGACCGTTGCTCAGGGAAAAA‑TCGCCGATCACATTGCAAAAACCG‑AGAATCAACCCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGc                                                                                                                               <  1:335672/251‑1 (MQ=255)
                                                                                        ggccgccCTTTTTTTC‑CCCCGCTGACCAGAAAATACCG‑CCGGGAGGGAAGAAGAAGGTGGCTTCCCCAACCAACCACAAAAAATCGTAGATCGTTTGGAGGAAGGGGTACATTTCCGCCGGCGGCGCCCCGTGGCTCAGGGGCAAA‑TCGCCAACCACATTGAAAAAACGG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCAt                                                                                                                            <  2:98198/250‑1 (MQ=255)
                                                                                            gccCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑GAGAAGAAGATCGCTTCACCGATCCACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATAACCGCGAAATAAGAGATCAACCGCTCATGCCATTTc                                                                                                                         >  2:211781/1‑250 (MQ=255)
                                                                                                        cccccgctaccaaGCAAAAACCGACGGGCAGG‑AAAAAGAAGATCGTTACACCGACCCACCAAAAAAAATCGAAAATCTTTGGCGGGGAAGGTTAAATTTCCGCCGGCGCCTGCCCTTGGCTAAGGGGCAAA‑CCGCCGACCACATTGCAAACCCCG‑AGAAACAAACCGCCGCTAACCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCgg                                                                                                             <  2:222132/239‑1 (MQ=255)
                                                                                                          acggCTGACCAGCAAATACCGACCGGCCGG‑CAGAAGAACATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCCGGGAAGGGTACCTTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCAAATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCCACGGGAAAAAGATCACCGCGAAATGATAGAACAACTGCTCATGCCATTTCATATTATGCCGggc                                                                                                           >  2:61264/4‑251 (MQ=255)
                                                                                                          aCCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGggc                                                                                                           >  2:392360/1‑251 (MQ=255)
                                                                                                                   cAGCAAATACCGACCGGCAGG‑GAGAAGAAGATCGCTTCACCGATCAACCACAGAAACTCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGCGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCtttttt                                                                                                  >  1:466877/1‑251 (MQ=255)
                                                                                                                    aGCAAATACCGACCGGCAGG‑TAGAAGGAGAACGCATCACCGCTCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGGCATCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGCAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCAtattttg                                                                                                                 >  2:212905/1‑232 (MQ=255)
                                                                                                                                   gCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATc                                                                                   >  1:314885/1‑250 (MQ=255)
                                                                                                                                      ataaagaagaagaTCGCTTCACCGACCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCCCCAGCGTCGGCCCTTTGCTCAGGGGAAAA‑TCGCCGATCACATTGAAAAAACCG‑AGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAAc                                                                                <  1:258365/247‑1 (MQ=255)
                                                                                                                                          agaagaagaTCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAAc                                                                                >  2:258365/1‑247 (MQ=255)
                                                                                                                                          aaaagaagaTCGCTACACCGATCACCCAAAAAAAATCGTAGATCGTTTGCAGGGAGGGGTACATTTGCCCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCAAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTgc                                                                            <  1:117970/249‑1 (MQ=255)
                                                                                                                                              gaagaTCGCTTCACCAATCAACCACAGAAAATCGAAGATCGTTTGCAGGGAAGGGTACTTTTGCGCCAGCGTCTGACCGTTGCGCAGGGGAAAA‑TCGCCGATCAAATTGGAAAAACCG‑AGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCt                                                                        <  1:238390/251‑1 (MQ=255)
                                                                                                                                                                           aaTCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGAt                                             >  2:216960/1‑249 (MQ=255)
                                                                                                                                                                                                            gcgcCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTg           >  1:101360/1‑250 (MQ=255)
                                                                                                                                                                                                             cgcCAGCTTCTGACCGTTGCTCATGGACAAA‑TCCCCGATCAAATTGCGAAAACCG‑AGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGaa         <  1:373796/251‑1 (MQ=255)
                                                                                                                                                                                                             ccccggcctctgcccgttgctgtggggCAAA‑CCCCCGATCACATTGCAAAAACCG‑AGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGaa         <  2:21238/227‑1 (MQ=255)
                                                                                                                                                                                                                    gTCTGCCCGTTGCTCATGGGCAAA‑TCCCCAATCACATGGCAAAACCCG‑AGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGccc  <  1:306175/251‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                   
ACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTT‑CACCGCTGACCAGCAAATACCGACCGGCAGG‑TAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAA‑TCGCCGATCACATTGCGAAAACCG‑AGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCC  >  NC_000913/4465537‑4465997

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