Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑490/+918) gatD ← / ← gatB galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3630C.84.8% 143.8 / 18.1 46intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/6);  new base . (33/6);  total (34/12)
Fisher's exact test for biased strand distribution p-value = 6.02e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3640C.84.8% 143.8 / 18.2 46intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/6);  new base . (33/6);  total (34/12)
Fisher's exact test for biased strand distribution p-value = 6.02e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173120‑2173613
                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                          
ccAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATg                                                                                                                                                                                                                                                    >  1:247055/1‑251 (MQ=255)
         cccgccAGTTGGGTGTGCAAACCAATAGTTTGCTTCGCGATCCAAAGGGTTATGCTCATAATGATGCCCCCCGAGATTAAGGTGCGGAACAGATTCCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATCCCACC‑‑GCCGTATGTCCcagcag                                                                                                                                                                                                                                           <  1:285833/250‑1 (MQ=255)
                       gtgCAAACCACTAGTTTGCGTCGCGATACACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGCTGCACTGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACATGGTTACCCGGCACACAAACAGCAA‑TTAACATGGTGAGTGGGATAAAAAACAGACTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCac                                                                                                                                                                                                                             >  2:505524/1‑250 (MQ=255)
                           aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GCTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                         >  1:223142/1‑251 (MQ=255)
                                       tGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCTGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATAGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGAACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑AATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGCCCAATCAGAAACTc                                                                                                                                                                                                             >  1:486017/1‑251 (MQ=255)
                                         cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGACGAAGCCGATGGAGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGCTAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTCTGTCCCAGCAACAACGCCGGATCACGGCCAATCAGGAACTccc                                                                                                                                                                                                           >  2:205107/1‑250 (MQ=255)
                                              cgATCCACAGGGTGATGCTCATAAGGATGACACCCGAGATTAAGGTGCGGAACATATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                                                                                      <  1:459453/251‑1 (MQ=255)
                                                agccACAGGGTGATGCTCTGAATGAGGACCCCCGAGATAAAGGTGCTGAACAGTTTTCCACGATGCACGGCGCCCGCCACCGCCACGAAGAAGCCGATGGTGCAAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACTCAAACAGCAT‑TTAAAATGGTGATTGGGATAAAACTCAGGCTGGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGATCTCCTGAccgcc                                                                                                                                                                                                    <  1:25437/249‑1 (MQ=255)
                                                 tCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GCTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                                   >  1:316482/1‑251 (MQ=255)
                                                   cacaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                                                                                       >  1:389710/1‑229 (MQ=255)
                                                        ggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                                            >  1:604670/1‑251 (MQ=255)
                                                              gCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                                                       >  2:608177/1‑250 (MQ=255)
                                                                                  agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                                   >  1:236939/1‑250 (MQ=255)
                                                                                     ttcaGGTGCGGAACAGATTTCCACGCTGCACGGCGACCGCAATCGCCACGGAGAAGCCGATGGCGGCAAGATCGCCAAACGACAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGTCTTGCC‑GGTACCAAC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAATCCTGACCGCCGAAATTACACTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                               >  2:197875/4‑251 (MQ=255)
                                                                                         ggTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGCTgg                                                                                                                                                           >  1:384974/1‑251 (MQ=255)
                                                                                         ggTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                                                                                                                                           >  1:419219/1‑251 (MQ=255)
                                                                                             cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCGGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGt                                                                                                                                                       >  1:171620/1‑251 (MQ=255)
                                                                                              ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                                                      >  1:275894/1‑251 (MQ=255)
                                                                                                          cccGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                          <  1:490274/249‑1 (MQ=255)
                                                                                                          cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCCAGATCGCCACACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGAATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCCGGCTTAGCGACGGGCGTTAAACACTCCatg                                                                                                                                          >  2:35886/1‑251 (MQ=255)
                                                                                                              agccccGGCGACCGCCTTCGCCCCGAAAAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCTCTTGATTCCCCGGCACACAAACAGCAA‑TAAAAATGGTGATGGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                                                                                                                                      <  1:112996/246‑1 (MQ=255)
                                                                                                                  aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                  >  1:384595/1‑251 (MQ=255)
                                                                                                                              accgccaaaaagaagCCGGTGGTGGCAAGACCGCAAAGCGGCACCACCTGATCACCCGGCCCACAAACAGCAA‑TAAAAAGGGGGGGGGGAATAAAAATCAGGCTTCCC‑GAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCATCCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACGACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTgg                                                                                                                     <  2:1555/242‑1 (MQ=255)
                                                                                                                                cGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCTTGGc                                                                                                                    >  2:98189/1‑250 (MQ=255)
                                                                                                                                       aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGCGAGTGGGATACAAATCAGGCTTGCC‑GATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATAATGGGTTTAATTCACC‑GTGGCATTagca                                                                                                             >  2:398496/1‑250 (MQ=255)
                                                                                                                                                ggggggacaaagctgcaaaaCGGGACCACCGGATTCCCCGGCCCAAAAACA‑CAACTTAAAATGGGGAG‑GGGTTTAAAAGCAGGGTTGCCCGTTAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTTAGCGTTGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTgc                                                                                                   <  2:223142/234‑1 (MQ=255)
                                                                                                                                                 tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTgcc                                                                                                  >  1:371221/1‑251 (MQ=255)
                                                                                                                                                           aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                         >  1:320368/1‑250 (MQ=255)
                                                                                                                                                               caaaagggCAGCACCGTATTCCCCGGCACAAAAACAGCAA‑TAAAAAGGGGGAGGGGGATAAAAATCAGGCTTCCC‑GATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGc                                                                                    <  1:455745/245‑1 (MQ=255)
                                                                                                                                                                               cgatccccgccccaaaaaCAGAAA‑TAAAAAGGGTGGGTGGAATAAAAACCGGGCTGCCC‑GAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCt                                                                    <  1:585119/236‑1 (MQ=255)
                                                                                                                                                                               ccttccccgatcacaaaaCAGAAA‑TAAAAAGGGTGGGTGGGATAAAATCCAGGCTTGCC‑GAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCt                                                                    <  1:682378/236‑1 (MQ=255)
                                                                                                                                                                                  tACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                                                 >  1:343569/1‑251 (MQ=255)
                                                                                                                                                                                    cccGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                               >  1:171613/1‑251 (MQ=255)
                                                                                                                                                                                    cccGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                               >  2:97578/1‑251 (MQ=255)
                                                                                                                                                                                     ccGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCCTGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCt                                                                    >  2:585119/1‑245 (MQ=255)
                                                                                                                                                                                          acacaAACAGCAA‑TTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                          >  1:85789/1‑250 (MQ=255)
                                                                                                                                                                                           cacaAAAAGAAA‑TTAAAAGGGTGAGGGGGAAAAAAACCAGGCTGGCC‑GAT‑CC‑CCCCGCCGTATGTCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAGCGTaa                                                        <  1:505524/251‑1 (MQ=255)
                                                                                                                                                                                            acaAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCt                                                                    >  2:682378/1‑238 (MQ=255)
                                                                                                                                                                                                aacgccatttAAAATGGTGAGGGGGTAAAAAATCAGTCTTGCC‑G‑‑ACCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcg                                                    <  1:214294/243‑1 (MQ=255)
                                                                                                                                                                                                    cAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAgg                                               >  1:548233/1‑251 (MQ=255)
                                                                                                                                                                                                            aaTGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                                                                                                                                                >  1:254110/1‑132 (MQ=255)
                                                                                                                                                                                                             aTGGTGAGTGGGCTAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCAAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                       >  2:239331/1‑251 (MQ=255)
                                                                                                                                                                                                                 tGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                   >  1:550759/1‑251 (MQ=255)
                                                                                                                                                                                                                  gAGTGGGATAAAAATCAGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGACCGCCGAACTTCGCCTGTAACCGCCTACGAGCCAGCTTAGAGATGGGCGTTAAACCATCCATGATGGGATTACT‑CACCCGTGGAATTAGCAGCATCACTCCCGCAGTTATTACCGCCAGCTGCACTAAACCTTTGACATCGTAACCCGCGAGGATGCCgcaaatc                                   >  2:571768/1‑244 (MQ=255)
                                                                                                                                                                                                                                  aGGCTTGCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAcc                  >  2:627988/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCC‑GATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGt            >  1:344170/1‑251 (MQ=255)
                                                                                                                                                                                                                                           gccaataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGt            <  1:97578/241‑1 (MQ=255)
                                                                                                                                                                                                                                                   caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCt    >  1:19198/3‑251 (MQ=255)
                                                                                                                                                                                                                                                     ccgccgTTTGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACTACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg  <  1:534837/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                          
CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAA‑TTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCC‑GATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAT‑CACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173120‑2173613

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