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.80.3% 192.5 / 47.5 66intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/12);  new base . (36/17);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 9.72e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01
*NC_0009132,173,3640C.80.3% 192.9 / 47.1 66intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/12);  new base . (36/17);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 9.72e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173127‑2173612
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                        
gaGCCGCCAGTTGGGTTTGCAAACCAATAGTTTGCGTGGCTATCCACAGGGGGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGACCGACAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCcag                                                                                                                                                                                                                                            <  2:643340/249‑1 (MQ=255)
  gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                           >  2:493241/1‑249 (MQ=255)
    cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGAAACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                        >  2:12369/1‑250 (MQ=255)
             ggTGTGCAAACCAATAGTTTGCGTGGCTATCCACAGGGTAATGCTAATAAGTATACCACCCGAGATAAAGGTGCGGAACAGATTCCCACGA‑GGCCCGGCGACCGCCACCGCCACGAAGAAGCCAATGGGGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGAAAAAAAT‑CAGTCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                                              <  2:137932/251‑1 (MQ=255)
                         aaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                   >  1:191305/1‑250 (MQ=255)
                          aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCACGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                   >  1:571034/1‑249 (MQ=255)
                             gTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                                                                                                                                                                                                               >  2:550205/1‑250 (MQ=255)
                                   gTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                                        >  2:636778/1‑251 (MQ=255)
                                      gggATCCACAGGGTGATGCTTATAATGATGACACCCGAGAGTAAGGTGCGGAACAGATTTCCCCGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                     <  1:11078/249‑1 (MQ=255)
                                       cgATCCACAGGGGGATTCCCATAATTAGGCCACCCAAAATAAAGGTGCGGAACAGATTCCCACGA‑GGCCCGGCGACCGCCATCGCCAGGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAC‑CAGGCTGGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                                                                                    <  2:702399/251‑1 (MQ=255)
                                             acaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                                                                                              >  1:9852/1‑251 (MQ=255)
                                                gggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                           >  2:721768/1‑251 (MQ=255)
                                                gggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                           >  1:201851/1‑251 (MQ=255)
                                                    gATGCCCATAATACTGACACCCGAGATTAAGGTACGGAACACATTTCCACGA‑TGAACGTCGACCGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAACCGGCAGCACATGATTACCCGGCACACAAACAGCAATTAAAA‑TGGAGAGTGGGATAAAAAT‑CAGACTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGACGGATCAAGGCCAATCAAGAACTCCTGACAGCCGCACTTCGCCt                                                                                                                                                                                       >  2:695548/1‑251 (MQ=255)
                                                          caaaatgataccacCCGAGATTAAGGTGCGAAACAGATTCCCACGA‑GGCACGCCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGACCGCAAAACGGCAGCACCTGATTACCCGGCAAACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                                                 <  2:8596/240‑1 (MQ=255)
                                                             aaTGCTGACCCCCGATAGTAAGGTGCGGAACAGATT‑CCCCGAGTCCCCGGCGCCCGCCATCCCCACGAAAAAGCCGATGGGGGCAAAAGCGCCAAACGGCGGCACCTGGTTACCCGCCACACAAGCAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGCGTTGCGGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAAGCAGAAACGCCTGTTCGCCGAACTTCGCCTGTAAACGAc                                                                                                                                                                              <  2:351548/251‑1 (MQ=255)
                                                              atgatgACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAC‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                             >  2:536901/1‑251 (MQ=255)
                                                               tgatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                             >  2:375996/1‑250 (MQ=255)
                                                                 aggaCACCCGAGATTAAGGTGCGGAACAGTTTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                                                           <  1:376105/248‑1 (MQ=255)
                                                                  ggaCACCCGAGATTAAGGGGCGAACCAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑GGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa                                                                                                                                                                         <  1:12369/250‑1 (MQ=255)
                                                                       cccGATATTAGGGTGCGGACAAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGAAGCACCTGATTCCCCGGCACTCAAACAGCAATTAAAA‑TGGTGAGGGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGcttg                                                                                                                                                                    <  1:498619/251‑4 (MQ=255)
                                                                             aTTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGc                                                                                                                                                              >  2:632498/1‑251 (MQ=255)
                                                                               tAAGGTGCGGAACAGCTTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCCACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCTTGCTTAGCg                                                                                                                                                             >  2:606873/1‑250 (MQ=255)
                                                                               aaaGGTGCGGAACAGATTCCCACGA‑TGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAACGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCg                                                                                                                                                             <  2:9852/249‑1 (MQ=255)
                                                                                aaGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGAt                                                                                                                                                           >  1:22257/1‑251 (MQ=255)
                                                                                 aGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATg                                                                                                                                                          >  2:614237/1‑251 (MQ=255)
                                                                                    tGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGc                                                                                                                                                       >  1:234555/1‑251 (MQ=255)
                                                                                      cGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGt                                                                                                                                                     >  1:417986/1‑251 (MQ=255)
                                                                                                ttCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                           >  2:479340/1‑251 (MQ=255)
                                                                                                 accACGA‑GGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGTGGGCAAGACCGCCAAACGGCAGCACCTAATTACCCGGCACACAACCGGCAATTAAAA‑TGGGGAGTGGGATAAAAAT‑CAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                           <  2:507983/249‑1 (MQ=255)
                                                                                                  ccACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                         >  2:718122/1‑251 (MQ=255)
                                                                                                      gA‑TGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTCAACCATCCatgatg                                                                                                                                     >  2:696865/1‑251 (MQ=255)
                                                                                                         tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                                                                   >  1:508362/1‑251 (MQ=255)
                                                                                                          gCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                                                                                                                                  >  2:231481/1‑251 (MQ=255)
                                                                                                            aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                >  2:196613/1‑251 (MQ=255)
                                                                                                                 gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGCTGGGTTTAATc                                                                                                                            >  1:465846/1‑250 (MQ=255)
                                                                                                                     cccATCGCCACGAAGAAGCCGATGGGGGAAAGAGCGCAAAACGGCAGCACCTAATTACCCGGCACACAAACAGAAATTAAAA‑TGGTGAGTGGGAAAAAAAT‑CAGGTTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                                                       <  2:191305/250‑1 (MQ=255)
                                                                                                                           gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                                 >  2:39267/1‑251 (MQ=255)
                                                                                                                            ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                >  2:242923/1‑251 (MQ=255)
                                                                                                                                  agaagCCGATGGTGGCAAGATCGCCAACCGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGTCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcagc                                                                                                          <  1:684704/251‑1 (MQ=255)
                                                                                                                                         gagggtggAAAGATCGCCAAACGGCAGCACCTAATTACCCGGCACACAAACAGAAATTAAAA‑TGGTGGGGGGGAAAAAAAT‑CAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTATCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTg                                                                                                   <  2:730370/248‑1 (MQ=255)
                                                                                                                                              gggCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                              <  1:93937/250‑1 (MQ=255)
                                                                                                                                                      ccgcaaaagggCAGCACCTGATTCCCCGGCACACAAACAGAAATAAAAA‑GGGGGAGGGGGA‑AAAAATGCAGGCTTGCGGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAc                                                                                      <  2:593795/242‑1 (MQ=255)
                                                                                                                                                        gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCg                                                                                    >  1:89658/1‑251 (MQ=255)
                                                                                                                                                            aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑GGGTGAGTGGGATAAAAAT‑CAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAGCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg                                                                                <  1:339213/251‑1 (MQ=255)
                                                                                                                                                                  agcccaggATTCCCCGGCACACAACCAGCAATTAAAA‑TGGTGATGGGGATAAAAAC‑CAGGCTTGCCGAAAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATTGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCaa                                                                          <  1:606873/244‑1 (MQ=255)
                                                                                                                                                                  agcaccgaattccccGGCACACAAACGGAAATAAAAA‑GGGTGAGGGGGAAAAAAAT‑CAGGCTGGCCAATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCaa                                                                          <  2:359875/239‑1 (MQ=255)
                                                                                                                                                                      cctgataccccGGCACCCAAACA‑TAATTAAAAAGGGTGAGGGTGATAAAAAC‑CAGGCTGGCC‑A‑ACCCCACCTCCTTATGTCCCAGCAGCACCGCCGTATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATaca                                                                      <  2:417986/243‑1 (MQ=255)
                                                                                                                                                                         gATTCCCCGGCACAAAAACAGCAATTAAAA‑GGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACAcc                                                                    <  1:493241/250‑1 (MQ=255)
                                                                                                                                                                         aaTTACCCGGCACACAAACAGCAATAAAAA‑TGGGGAGTGGTATAAAAAC‑CAGGTTTCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCt                                                                   <  2:201851/250‑1 (MQ=255)
                                                                                                                                                                            tACCCGGCACACAAACAGCAATAAAAA‑GGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                                 <  1:721768/250‑1 (MQ=255)
                                                                                                                                                                            accccGGCACAAAAAAAGCAATAAAAA‑TGGTGAGTGGGAAAAAAAT‑CAGCCTGTCCGAT‑CC‑CCCCGCCGTATGCCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                                                <  2:22257/249‑1 (MQ=255)
                                                                                                                                                                              cccGGCACACAAACAGCAATTAAAA‑GGGTGAGTGGGATAAAAAC‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                              <  1:39267/251‑1 (MQ=255)
                                                                                                                                                                                ggccccaaaaacacaaataaaaa‑GGGTGAGGGGGAAAAAAAT‑CAGGCTTCCCGA‑ACC‑CCCCGCCGTATGTCCCGGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                            <  2:235068/232‑1 (MQ=255)
                                                                                                                                                                                     cacaAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATA‑‑ACACCGCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCTTCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                                        >  1:317873/1‑250 (MQ=255)
                                                                                                                                                                                           cAGCAATTAAAA‑GGGTGAGTGGGATAAAAAT‑CGGGCTTGCCTAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                                                 <  2:71678/251‑1 (MQ=255)
                                                                                                                                                                                               aaTTAAAA‑GGGTGAGGGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                                             <  1:550205/251‑1 (MQ=255)
                                                                                                                                                                                                    aaaaggggaggaggaaaaaaaaT‑CAGGTTGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                         <  2:508362/234‑1 (MQ=255)
                                                                                                                                                                                                    aaa‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGAATGTCCCAGCAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCACACTCCGCATGTAAACGACTACGAGCTTGCTTAGCGCTGGGCGTTAACCCATCCATGATGGGTTTAATCACCCGAGGCATTAGCAGCATCACTGCCGCAGTTTTTACCGCAAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGACGCCg                                        >  1:700033/1‑251 (MQ=255)
                                                                                                                                                                                                         ggTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGATTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATa                                     >  1:623198/1‑250 (MQ=255)
                                                                                                                                                                                                          gTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAAt                                   >  1:319418/1‑251 (MQ=255)
                                                                                                                                                                                                            gAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                  >  2:619092/1‑250 (MQ=255)
                                                                                                                                                                                                              gTGGGATAAAAAT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGc                               >  2:50959/1‑251 (MQ=255)
                                                                                                                                                                                                                     aaaaaT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGCTAATCAGCCCCATc                         >  2:639837/1‑250 (MQ=255)
                                                                                                                                                                                                                        aaT‑CAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca                      >  2:658576/1‑250 (MQ=255)
                                                                                                                                                                                                                                cTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAg                                >  1:117962/1‑233 (MQ=255)
                                                                                                                                                                                                                                       cgatacccccGCCGGATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCTCTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTCGACATCGTAACCCGCGAGGAAGCCGATAATCAGCCCCATCACAAAACCCACGGTGa         <  1:196613/245‑1 (MQ=255)
                                                                                                                                                                                                                                           cccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAg      <  1:375996/246‑1 (MQ=255)
                                                                                                                                                                                                                                           accaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAGAATCAGCCCCATCACAAAACCCACGGTGACAg      <  1:231481/246‑1 (MQ=255)
                                                                                                                                                                                                                                              cccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  <  1:614237/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                        
TAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGA‑TGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAAT‑CAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173127‑2173612

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