Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.93.2% 311.0 / 9.0 73pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (51/17);  total (51/22)
Fisher's exact test for biased strand distribution p-value = 1.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.97.1% 308.8 / 0.5 73pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (51/17);  total (51/22)
Fisher's exact test for biased strand distribution p-value = 7.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCACCCCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCA  >  NC_000913/2173125‑2173599
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                              
aTTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATaccac                                                                                                                                                                                                                                                 <  1:344063/237‑1 (MQ=255)
    gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATaccaccgccgt                                                                                                                                                                                                                                           >  1:845666/1‑234 (MQ=255)
         cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGt                                                                                                                                                                                                                                    >  1:1500257/1‑238 (MQ=255)
             ggggTGGGCAACCCAATAGTTGGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCTACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTccc                                                                                                                                                                                                                                 <  1:64731/236‑1 (MQ=255)
                                  tGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                          >  1:1319003/1‑239 (MQ=255)
                                  tGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGAATCAAGGcc                                                                                                                                                                                                          >  1:784165/1‑239 (MQ=255)
                                    cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                          >  1:1434388/1‑237 (MQ=255)
                                       cgcgATCCACGGGGTGATGCCAATAATAGTGACCCCCGAGATAAAGGCGCGGAAAAGTTTCCCACGATGCAGGGCGCCCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                      <  2:1319003/238‑1 (MQ=255)
                                                     cgATCCTCATAAGGATGCCCCCCGAGATGAAGGTGTGAAGCAGATTCCCACGATGCACGGCGACCGCCATCTCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGTTGAGTGGTATAAAAATCA‑GCCTGGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                                 <  2:1239734/212‑1 (MQ=255)
                                                      gATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGCAACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                                                                       >  1:1562807/1‑238 (MQ=255)
                                                      gATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                                                                       >  1:20357/1‑238 (MQ=255)
                                                        tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                    >  1:422891/1‑239 (MQ=255)
                                                        tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                    >  1:234318/1‑239 (MQ=255)
                                                        tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAAAGCCGGATCAAGGCCAATCCGGAACTCCTGAccgccg                                                                                                                                                                                    >  1:1109659/1‑239 (MQ=255)
                                                               aaTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                                       >  1:1239734/1‑213 (MQ=255)
                                                                  gatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTACGTCCCAGCAGAAACGCCGGATCAAGACCAATCAGGAACGCCACACCGCCGAACATCg                                                                                                                                                                             >  2:914969/1‑236 (MQ=255)
                                                                     gccacCCGAGATTAAGGTGCGGAACAGATTACCACGATGCACGGCGCCCGCCATCGCCAGGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                                       <  1:1735740/237‑1 (MQ=255)
                                                                                    ggTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                                                                                                           >  1:1516711/1‑236 (MQ=255)
                                                                                        cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                    >  1:1300677/1‑239 (MQ=255)
                                                                                              agtttaccACAAGGCACGCCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCCCCGGATTACCCGGCACACAAACAGCAATAAAAATGGTGAGGGGAAAAAAAATCA‑GGCTGGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                              <  2:422891/233‑1 (MQ=255)
                                                                                                     cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                                         <  1:1435463/173‑1 (MQ=255)
                                                                                                     cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                                         >  2:1435463/1‑173 (MQ=255)
                                                                                                             cgcgCGACCGCCCTCGCAAGGAGGACCCCGAGGGGGGCAAGACCTCAAAGGGGCAGCGCCTGATACCCCGTCACACAACCAG‑AATGAAAAGTGGGAAGGGGGTAAAAAGCAGGGCTGGCCGATAC‑‑CGCCGCC‑GTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                                                                                                                               <  2:1562807/236‑1 (MQ=255)
                                                                                                                   gccaccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                         >  1:1849847/5‑239 (MQ=255)
                                                                                                                    ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                                                    <  1:1602034/163‑1 (MQ=255)
                                                                                                                    ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                                                    >  2:1602034/1‑163 (MQ=255)
                                                                                                                           cGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                                                                                                                                   >  1:953601/1‑221 (MQ=255)
                                                                                                                                   agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                                    >  1:609648/1‑212 (MQ=255)
                                                                                                                                       gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCa                                                                                                     >  1:40840/1‑239 (MQ=255)
                                                                                                                                            tgggggAAAGATCGCCAAACGGCAGCACCTAATTCCCCGGCCCACAAACAGAAATTAAAAGGGGGAGTGGGAAAAAAATCA‑GGCTGGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag                                                                                                 <  1:556250/238‑1 (MQ=255)
                                                                                                                                                 gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                >  1:1053872/1‑186 (MQ=255)
                                                                                                                                                 gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                >  1:815769/1‑186 (MQ=255)
                                                                                                                                                         gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                                          >  2:312635/1‑216 (MQ=255)
                                                                                                                                                         gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                                          <  1:312635/216‑1 (MQ=255)
                                                                                                                                                           cAAACGGCGGACCCTGATTACCCGGAACACAAACGACAATAAAAATGGTGAGTGGCATAAAAACCA‑GGCTTGCCGACAC‑‑CGCCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                                         <  2:1745496/119‑1 (MQ=255)
                                                                                                                                                             aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGc                                                                                                      >  1:723762/1‑216 (MQ=255)
                                                                                                                                                             aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGt                                                                               >  2:338216/1‑239 (MQ=255)
                                                                                                                                                                     caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                                               >  1:1745496/1‑119 (MQ=255)
                                                                                                                                                                       ccTGATTCCCCGGCCAAAAAACAGAAATTAAAAGGGTGAGTGGGAAAAAAACCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                     <  1:1588191/239‑1 (MQ=255)
                                                                                                                                                                       ccTGATTCCCCGGCACAAAAACAGAAATAAAAAGGGTGAGGGGAAAAAAAATCA‑GGCTTGCCGA‑ACC‑CCCCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                     <  1:59507/239‑1 (MQ=255)
                                                                                                                                                                        cttaTTCCCGGGCAAAAAAAAAACAATTAAAATGTTGAGGGGGAAAAAAACCA‑GGCTTGCCGATCCCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                                          <  2:734458/182‑1 (MQ=255)
                                                                                                                                                                           atgaccggccacaaaAACGGCAAGAAAAAGGGGGAGTGGAATAAAAAACG‑GGCTTGCCAATCCCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                         <  2:1107513/174‑1 (MQ=255)
                                                                                                                                                                            ttccccGGCACACAAACAGCAATAAAAATGGTGGGTGGGAAAAAAATCA‑GGCTGCCCG‑‑ACCCCCCCGCC‑GTATGTCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                          <  2:290956/146‑1 (MQ=255)
                                                                                                                                                                               cccGGCCCCCAAACAGCAATTAAAAGGGTGGGTGGGATAAAAATCA‑GGCGTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGc                                                                                                                                                                                                                      <  2:1725157/86‑1 (MQ=255)
                                                                                                                                                                                 cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                            >  1:1023148/1‑238 (MQ=255)
                                                                                                                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                                >  1:734458/1‑185 (MQ=255)
                                                                                                                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTa                                                                          >  1:134259/1‑223 (MQ=255)
                                                                                                                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATa                                                           >  1:1399517/1‑238 (MQ=255)
                                                                                                                                                                                  ggCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCAGCCGTTTTTACCGCCAGCTGCAATaa                                                          >  1:1967867/1‑238 (MQ=255)
                                                                                                                                                                                    cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATaca                                                         >  1:1301620/1‑238 (MQ=255)
                                                                                                                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCCTGCTTAGCg                                                                                                                                                >  1:2021739/1‑150 (MQ=255)
                                                                                                                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                               >  1:1107513/1‑183 (MQ=255)
                                                                                                                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGCAAACGACTACGAACCTGCTTAGCGATGGGCGTTACACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATAAATGCCGCCGTTTTTACCGCCAGCCGCacacacca                                                       >  2:1029398/1‑232 (MQ=255)
                                                                                                                                                                                     acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGAAACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCCTAGCGATGGGCGTTAAACCATCCATGATGGGTCTAATCACACGTGGCATTATCAGCATCACTACCGCCGTTTTTACCGCCAGCTGCAATaat                                                         >  2:237212/1‑235 (MQ=255)
                                                                                                                                                                                      cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                >  1:290956/1‑149 (MQ=255)
                                                                                                                                                                                         aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCcagcag                                                                                                                                                                                                                            >  2:1865728/1‑70 (MQ=255)
                                                                                                                                                                                         aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCcagcag                                                                                                                                                                                                                            <  1:1865728/70‑1 (MQ=255)
                                                                                                                                                                                         aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GGATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                            >  1:1725157/1‑86 (MQ=255)
                                                                                                                                                                                               cAATAAAAATGGGGAGGGGGATAAAAACCA‑GGCTTGCCGATCCCA‑‑CCGCC‑GTATGCCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                              <  2:112036/142‑1 (MQ=255)
                                                                                                                                                                                                aaTTAAAAGTGGGAGGGGGAAAAAAAACC‑GGGCTTCCGA‑ATCCCACCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                           <  2:1165167/177‑1 (MQ=255)
                                                                                                                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATaccac                                                                                                                                                                                                                                                 <  1:914421/44‑1 (MQ=255)
                                                                                                                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATaccac                                                                                                                                                                                                                                                 >  2:914421/1‑44 (MQ=255)
                                                                                                                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                                                                                                            <  1:108285/110‑1 (MQ=255)
                                                                                                                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                                                                                                            >  2:108285/1‑110 (MQ=255)
                                                                                                                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGc                                                                                                      >  1:1239804/1‑180 (MQ=255)
                                                                                                                                                                                                 aTTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                           >  1:1886907/1‑239 (MQ=255)
                                                                                                                                                                                                         ggTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                                    >  1:112036/1‑142 (MQ=255)
                                                                                                                                                                                                         ggTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag                                                                                                 >  1:1165167/1‑177 (MQ=255)
                                                                                                                                                                                                         ggTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCACTAAACCCTTGACATCGTAACCcgcca                                   >  2:2114865/1‑237 (MQ=255)
                                                                                                                                                                                                          gTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                                                      >  2:691764/1‑155 (MQ=255)
                                                                                                                                                                                                          gTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                                                      <  1:691764/155‑1 (MQ=255)
                                                                                                                                                                                                             aGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                               >  1:501298/1‑239 (MQ=255)
                                                                                                                                                                                                             aGTGGGATAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                               >  1:133014/1‑239 (MQ=255)
                                                                                                                                                                                                                    tAAAAATCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATa                        >  1:1055517/1‑239 (MQ=255)
                                                                                                                                                                                                                       aaaTCA‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATAGTAACCAGCGAGGCTGCCGATAAta                     >  1:1699733/1‑238 (MQ=255)
                                                                                                                                                                                                                            a‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcc                 >  1:284868/1‑238 (MQ=255)
                                                                                                                                                                                                                            a‑GGCTTGCCGATACCA‑‑CCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGAAACCCGCGAGGATGCCGATAATCAGcc                 >  1:1853683/1‑238 (MQ=255)
                                                                                                                                                                                                                                cgcttagccaacaccCCGCCGGTATGCCCCGGCACAACCCCCGGATCAAGGCCATCCAGAACCTCCGGCCCGCCGACCTTCCCCTGTAAACGACTACGAGCCTGCTTGGCGGTGGGCGTTAACCCAGCCATGAGGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGccc                <  2:1516711/228‑1 (MQ=255)
                                                                                                                                                                                                                                             ccccGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCa  <  1:2114865/239‑1 (MQ=255)
                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                              
ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGCTTGCCGATACCACCCCGCC‑GTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCA  >  NC_000913/2173125‑2173599

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