Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3640C.78.6% 99.5 / 32.5 42intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/8);  new base . (27/6);  total (28/14)
Fisher's exact test for biased strand distribution p-value = 1.93e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173129‑2173606
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGCCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                      >  1:427833/1‑249 (MQ=255)
          gggTGTGCAAACCAATAGTTTGCGTCGCGATCAACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                                                           >  1:296095/1‑250 (MQ=255)
                  aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGCTGCTCCTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGCGTGGGATAAAAATCAGGCT‑TGCCGACACCACC‑‑GCCGTATGTCCCAGCAGCAACCCCGGATCAAGGcc                                                                                                                                                                                                                  >  2:272284/1‑251 (MQ=255)
                      cAATAGTTTGCGTCGCGATCCCCAGGGTGAGGCTCAGAATGATGTCCCCGGTGATTAAGGTGCGGAACGGATTTCCTGGATGCACGGCGACCGCCAGCGCCTTGAAGAAGCCCAGGGTTGCAAGATCGTGAACGGGCGGCACGTGATTACCAGGCCCACAAACAGGAATTAAAATGGTGAGTGGGATAAAAAGCAGTCT‑TGCCGATGCGACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAATGCCAATc                                                                                                                                                                                                              <  1:330164/251‑1 (MQ=255)
                      cAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGACGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                              >  2:362309/1‑251 (MQ=255)
                                 gTCGCGATACACAGGGTGCTGCTCATAATGATGACACCAGAGATTAAGGAGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGACGATGGTGGCAAGATCGCCAAACGGCAGCAACTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAAGCT‑TTCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                    >  2:288989/1‑250 (MQ=255)
                                     cgATCCACAGGGTGATGCCCATAATGATGACACCCGGGATTAAGGTGCGGAACAGTTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAGCGGCAGCACCTTATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAGCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                                                                                                                                 <  1:267529/249‑1 (MQ=255)
                                                   aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGACCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                                                                                                                                 >  1:440306/1‑251 (MQ=255)
                                                          tAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg                                                                                                                                                                           >  1:269468/1‑250 (MQ=255)
                                                                         agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                            >  2:404728/1‑250 (MQ=255)
                                                                           ataaaGGGGCGAAACAGATTCCCCCGAGGCACGCCGACCGCCACCGCCACAAAGAAGCCGAGGGTGGAAAGATCGCCAAACGGCGGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGAATAAAAATCAGGCT‑GGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                                                                                                                                         <  2:341095/248‑1 (MQ=255)
                                                                               aGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                                                                                                                                     >  2:84187/1‑251 (MQ=255)
                                                                                       aaCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                                                                                                                              >  1:69842/1‑250 (MQ=255)
                                                                                           gttttCCACGATGCAAGAAGACCGCCATCGCCACGAATAAGCCGATGGCGGCAAGCTCGCCAAACGGCAACAACTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GACGAATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCATGACCGCCGAACTTCGCATGTAAAAGCCTACGAGCGTGCTTAGCGATGGGCGTTACAcc                                                                                                                                          >  2:270319/3‑250 (MQ=255)
                                                                                            attcccccgaggCACGGCGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                        <  1:12615/240‑1 (MQ=255)
                                                                                            aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAAGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCACCGCCGGATCAAGGCCAAACAGGAACTCCTGACCGCCGCACTTCGCCTGTACACGACAACGAGCCTGCTTAGCGAGGGCCGTTAAAccat                                                                                                                                        >  2:82499/1‑251 (MQ=255)
                                                                                                   cGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                                 >  1:402657/1‑251 (MQ=255)
                                                                                                              gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGCAAACGACTACGAGCCCGCTTAGCGAATGGCGTTAAACCATCCATGATGGGTTT‑AAtac                                                                                                                     >  2:163115/1‑249 (MQ=255)
                                                                                                              gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATc                                                                                                                      >  1:93825/1‑250 (MQ=255)
                                                                                                              gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGAAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATc                                                                                                                      >  2:113709/1‑250 (MQ=255)
                                                                                                               aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCAc                                                                                                                    >  1:276968/1‑251 (MQ=255)
                                                                                                                  gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCg                                                                                                                 >  2:159209/1‑251 (MQ=255)
                                                                                                                              aaaaacccgttggtggCAAGATCGCCAAAGGGCACCCCCTGATTACCCGCCCCACAAACAGCAATTAAAATGGTTTGTGGGATAAAAATCAGGTT‑GGCCGATAC‑‑CCCCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTagcag                                                                                                     <  2:409048/241‑1 (MQ=255)
                                                                                                                                 aagCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCAt                                                                                                  <  1:343688/251‑1 (MQ=255)
                                                                                                                                    ccGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCAc                                                                                               <  1:45380/251‑1 (MQ=255)
                                                                                                                                              aaagagccccAACCGGCAGCACCTGTTCACCCGTCACACAACCAGCATTTAAATGGGTGAGTGTGAAAAAAACCGGGCT‑TGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACTCCGGATCAAGGCCAACTAGGAACTCCTGTCCGCCGAGTTTCTCGTGTAAACGGCTACTAGCCTGCTTGGCGATGGGCGTTAAACCATGCATTA‑GGGTTTGAATCGCCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                     <  1:273818/243‑1 (MQ=255)
                                                                                                                                                 gaccGCCAAACGGCATCCCCTGATTACCCGGCACAAAAACAGCAATAAAAGTGGGGAGTGGGATAAAAATCAGTCT‑TGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                  <  1:402254/248‑1 (MQ=255)
                                                                                                                                                  cccgcaaagcggcctcacctatgtACCCGGCCCAAAAACAGAAATAAAAAGGGGGAGGGGGATAAAAA‑CCTGCTGTGCGGATAC‑‑CCCCGCTGTCTGCCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGTTTATCGATGTTCGTTAAACCATCCATGATGGGTTG‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTa                                                                                 <  2:159987/228‑1 (MQ=255)
                                                                                                                                                          aCGGCAGCACCTGATTACCCGGCACACCAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                                         >  2:408286/1‑251 (MQ=255)
                                                                                                                                                                  aCCTGATTACCCGGCACACAAACCGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCACTac                                                                 >  2:138448/1‑251 (MQ=255)
                                                                                                                                                                   ccTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                                 >  1:261002/1‑250 (MQ=255)
                                                                                                                                                                           cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                        >  2:201108/1‑251 (MQ=255)
                                                                                                                                                                               gCACAAAAACAGCACTTAAAATGGTGAGAGGGATAAAAATCAGGAT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCACCGCCGGATCAAAGCCAATCCGGAACTCCTGACCGCCGAACTTCACCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTACACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGACGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                    >  2:214656/1‑251 (MQ=255)
                                                                                                                                                                                   acaAACAGCAATTAACATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGAATTAGAAGCATCACTGCCGCAGTTTTTACCGCAAGCTGCAATACACCTTTGACATCGTAAc                                                >  1:337014/1‑251 (MQ=255)
                                                                                                                                                                                   aaaaaCAGAAATAAAAAGGGTGAGTGGGATAAAAATCGGGCT‑TGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                                <  1:404728/249‑1 (MQ=255)
                                                                                                                                                                                       aaaGCAATAAAAAGGGGGAGTGGTATAAAAGTCGGCCT‑TGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgc                                            <  1:341676/249‑1 (MQ=255)
                                                                                                                                                                                           cAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAgg                                        >  2:392749/1‑251 (MQ=255)
                                                                                                                                                                                                   aTGGTGAGGGGAATAAAAACCAGGCT‑GGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                <  1:408286/251‑1 (MQ=255)
                                                                                                                                                                                                               aaaaaaaTCAGGTT‑GGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCAt                    <  1:392749/249‑1 (MQ=255)
                                                                                                                                                                                                                tAAAAATCAGGCT‑TGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc                   >  1:85389/1‑251 (MQ=255)
                                                                                                                                                                                                                  aaaaTCAGGCT‑TGCCGATACCACC‑‑GCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac                 >  1:332121/1‑251 (MQ=255)
                                                                                                                                                                                                                            t‑tGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCTTTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCCTCACAAAACCCACg       >  1:354145/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gccggtccccccGCCGTATGTCCCGGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCGGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTATCCCTCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGt     <  1:397003/243‑1 (MQ=255)
                                                                                                                                                                                                                                    gataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACACAACCCACGGTGAc  >  1:230539/8‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
GCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173129‑2173606

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