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

ACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGG  >  NC_000913/4465537‑4466028
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
acTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCg                                                                                                                                                                                                                                                   <  1:158076/251‑1 (MQ=255)
    ccGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGaaaa                                                                                                                                                                                                                                               <  2:632300/251‑1 (MQ=255)
      ggCAGCTGCTGCCCGGGAAGATATGCGTTCATAAGCTGTGGAGATACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTATAAGAAGATCGCTTCACCGATAAACACC‑AGAAAATCGGAGATCGTTTGCAGGGGAAGGT‑ACATTTGAGCCCGCGTTTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAAcc                                                                                                                                                                                                                                             >  2:220841/1‑251 (MQ=255)
         aGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTGTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAg                                                                                                                                                                                                                                          <  2:129212/251‑1 (MQ=255)
             gctgcCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAATAAGATCGCTTCACCGATCAACCCA‑AGAAAATCGTAGATCGTTTGCCGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGGTTCCCCTGGGCA‑AATCGCCGATCACATTGCGAAAACCgggcagc                                                                                                                                                                                                                                      >  2:460956/1‑245 (MQ=255)
                gcCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGGAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGAACaa                                                                                                                                                                                                                                    >  2:233100/1‑250 (MQ=255)
                 cccGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAAc                                                                                                                                                                                                                                  <  1:233100/251‑1 (MQ=255)
                                gTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGGTCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGAGTCAAACCGCCGCTAATCCAcg                                                                                                                                                                                                                   >  2:63680/1‑251 (MQ=255)
                                                   cccAGTGTCCCGCCAACCACAATCCCAAGGACCGGCGTCCCCCCTATTTTTTCCCCCGCTAACAAGCAAAACCCGCCCGGCAGGTAGAAGAAGACCGCTTCACCGATCAACAAA‑AGAAAACCTTAGATCGTTTGCGGGGAGGGTT‑CCATTTGCCCCACCGTCTGCCGGTGGCTCAGGGGCA‑AATCGCCGAACCAATGGCAAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaagaa                                                                                                                                                                                                <  1:17096/250‑1 (MQ=255)
                                                    ccAGTGCCACGCAAAGCACTATACAAGGGATGGGCGCCCCCCCTTTTTTTTTCACCGCTGACCAGAAAATCCCGCCCGGCAGGTAAAAGAAGATCGCTCCCCCGATCAACCCC‑AGAAAATCGTAGATCGTTTGCAGGGAGGGTT‑ACATTTGCGCCAGCGTCTGCCGGTTGCTCAGGGGAA‑AATCGCCAATCAAATGGCAAAAACCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGaagaag                                                                                                                                                                                               <  1:15871/251‑1 (MQ=255)
                                                          tCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCAcc                                                                                                                                                                                         >  1:460568/1‑251 (MQ=255)
                                                                  accactaTACAAAGGATCGGCGTGCCGCCCATTTTTTCCACCGCTGCCCAGCAAATACCGCCCGGCAGGTAGAAGAAGACCGCTTCCCCGATCAACCAC‑AAAAAATCGTAGATCGTGTGC‑GGGAAGGGTTACATTTGCGCCGGCGTCGGCCCGTTGCTCAGGGGCA‑AAGCGCCGATCACATTGCGAAACCCGAGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATg                                                                                                                                                                                 <  1:229997/245‑1 (MQ=255)
                                                                        atacaaaggagaggcggccccaccATTTTTTTCACCGCTGACCACAAAATACCGCCGGGCGGGTAAAAGAAGATCGCTTCCCCGAAAAACAAA‑AGAAAACCGTAGATCTTTTCGAGGGAAGGTT‑ACATTTGCACCGGCGTCTGACGGTTGGTCATGGGCA‑ACTCGCCAATCACATTGCAAAACCCGGGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGAt                                                                                                                                                                           <  1:347778/229‑1 (MQ=255)
                                                                          aCCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACCTTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGGGAGATc                                                                                                                                                                          >  1:217564/1‑250 (MQ=255)
                                                                             aaGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATAGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCAGCTAATCAACGCGGgcagca                                                                                                                                                                                                          >  1:155270/1‑215 (MQ=255)
                                                                                 acgggCGTGCCGCCCATTTTTTCCCCCGCTGACCAGAAAATCCCGACCGGCAGGAAAAAGAAGACCGCTTCCCCGATCAACCAA‑AGAAATTCGTAGTTCGTTTGCAGGAAAGGGT‑CCATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGACCACATGGCGAAAACGGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTc                                                                                                                                                                  <  1:207298/248‑1 (MQ=255)
                                                                                       gtcccgccCATTTTTTCCCCCGCGGCCCAGCAAACCCCGCCCGGCAGGTAAAAGAAGATCGCTCCCCCGACCACCCAA‑AGAAAATCGTAGATCGTTTGAAGGAAGGGGT‑ACATTGGCCCCAGCGTCTGACCGTGGCTCAGGGGAA‑AACCGCCGATCACATTGCAAAACCCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCa                                                                                                                                                            <  1:26123/248‑1 (MQ=255)
                                                                                       gTGCCGCCCTTTTTTTCCCCCGCGGACCAGCAAACCCCGACCGGCAGGAAGAAGAAGACCGCTCCCCGGATCACCCAC‑AAAAAATGGTGGTTCGTTTGCGGGGAGGGGT‑ACTTTGGCGCCAGCGTCTGACCGTGCCCCGGGGGAA‑AACCGCCGACCACTTGGAAAAAACGGAGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCa                                                                                                                                                            <  2:231065/251‑1 (MQ=255)
                                                                                          ccgccCATTTTTTTCACCGCTGACCAGCAAATCCCAACCGGCAGAAAAAAGAAGATCGCTTCACCGACCAACCAA‑AGAAAATCGAAGATCTTTTGCAGGGAAGGGT‑ACATTGGCGCCAGCGTCTGACGGTTGCTCAAGGGAA‑AACCCCCGATCACATGGCAAAAACCGAGTATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCAttt                                                                                                                                                         <  1:358321/251‑1 (MQ=255)
                                                                                                aTTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCAtatt                                                                                                                                                   >  1:311717/1‑251 (MQ=255)
                                                                                                  tttttcccccgccgcccaacaattACCGACCGGAAGGTAGAATAAGACCGCTTCACCGCCAAACCAC‑AAAAAAGTGTAGTTCTTTTGCAGGGAAGGTT‑CCATTGGCGCCCGCGACGGACCTTTGCTCAGGGGCA‑AATCGCCATTCACATTGCAAAACCCGAGAATAACACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCAtatta                                                                                                                                                  <  1:239865/227‑1 (MQ=255)
                                                                                                         cccGCTGACCACAAAATCCCGACCGGGAGGTAAAAAAAGACCGCTCCACCGA‑CAACCCCCAGAAAACCGTGGATCGTTTGCAGG‑AAGGGTTCCATTTGCGCCAGCGTCGGCCGGGTGCGCGTGGGAA‑AGCCGCCGCCCACATTGCAAAAACCGAGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGggc                                                                                                                                          <  1:542836/250‑1 (MQ=255)
                                                                                                           cGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGG‑ACATTTGCGCCAGCGGCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGACGATAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCCCAATCCATTTCaacttaagccgcgcgg                                                                                                                                        >  2:361267/1‑236 (MQ=255)
                                                                                                                  cagaaaaTCCCGACCGGCAGGTAAAAGAAGATCGCTTCACCGATAACCCAA‑AAAAATTCGTAGATCGTTTGCGGGGAGGGGT‑ACATTTGCGCGAGCGTCTGACGGTTGCTCAGGGGCA‑AACCGCCGATCACATTGCGAAAACCGAGAATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCtttttt                                                                                                                                 <  1:603863/247‑1 (MQ=255)
                                                                                                                       aaTACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCt                                                                                                                            >  1:125737/1‑251 (MQ=255)
                                                                                                                        aTACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTg                                                                                                                           >  1:286406/1‑251 (MQ=255)
                                                                                                                          accgaccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAAAGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGAGCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGct                                                                                                                         >  1:113869/1‑251 (MQ=255)
                                                                                                                                    gggTAGAAGAAGTCCGCTCCACCGACCAACCAA‑AAAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACTTTTGCGCCGGCTTCTGCCCGTTGCTCATGGACA‑AACCGCCGATCACATTGCAAAACCCGAGAATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAAc                                                                                                               <  1:273486/250‑1 (MQ=255)
                                                                                                                                              agaTCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGt                                                                                                     >  1:270053/1‑251 (MQ=255)
                                                                                                                                                                       gAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTg                                                                              >  1:582835/1‑250 (MQ=255)
                                                                                                                                                                             cGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTa                                                                        >  1:330140/1‑250 (MQ=255)
                                                                                                                                                                                          cAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTgg                                                          >  2:111199/1‑251 (MQ=255)
                                                                                                                                                                                             ggAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGc                                                       >  1:462/1‑251 (MQ=255)
                                                                                                                                                                                                                   gctataaCCGTCGCTCATGG‑CACAATCGCCAATCACATGGCGAAAACCGAGGATCAAACCGCCGCTACTCAACGCGGGCAGCAACGAGACGAAGATCCCCGCGAACTGAGAGATCCACTGCCCATGCCACTTCATATTATGCCGGGCGGCTCCTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGACAGAAATCACCCACGTTGGTGCCAATCACCACCTgccattcac                                  >  2:603007/7‑243 (MQ=255)
                                                                                                                                                                                                                      cTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAACCCGAGGATCAAACCGCCGCTACTCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGCAATTGCCCgg                               >  2:493504/1‑251 (MQ=255)
                                                                                                                                                                                                                             ttGCTCAGGGACA‑ACTCGCCGATCACATGGCAAAACCCGAGGACAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTgg                         <  1:99405/250‑1 (MQ=255)
                                                                                                                                                                                                                                cctcatggaaaaatccccaaacaaaTTGAGAAAACCGGGGAGAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTg                       <  1:220841/226‑1 (MQ=255)
                                                                                                                                                                                                                                         a‑aaTCGCCGATCACATTGCGAAAACAGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCttt            >  1:186399/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    ccacaTGGCGAAAAACGAGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAggg  <  1:169436/250‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
ACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCAC‑AGAAAATCGTAGATCGTTTGCAGGGAAGGGT‑ACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCA‑AATCGCCGATCACATTGCGAAAACCGAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGG  >  NC_000913/4465537‑4466028

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: ‑