Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3600A.21.4% 103.7 / 15.8 42intergenic (‑487/+922)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base A (27/6);  new base . (1/8);  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 = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACG  >  NC_000913/2173114‑2173601
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
aaCGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGGGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATaccaccgc                                                                                                                                                                                                                                                <  2:360870/251‑3 (MQ=255)
               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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
AACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCT‑TGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTT‑AATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACG  >  NC_000913/2173114‑2173601

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