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.97.1% 148.1 / ‑3.4 35intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (23/11);  total (24/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01
*NC_0009132,173,3640C.97.1% 147.9 / ‑3.4 35intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (23/11);  total (24/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGC‑TCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173131‑2173613
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                         
cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGC‑TCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                         >  1:404691/1‑250 (MQ=255)
   cAGTTGGGTGAGCAAACCAATAGTTTGCGTCGCGCTCCACAGGGTGATGC‑TCATAATGATGACACCCGAGATTAAGGTGCGGCACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCAGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGAAACCACC‑‑GCCGTATGTCCCAGCAGCAACg                                                                                                                                                                                                                                     >  2:429899/1‑251 (MQ=255)
     gTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGGGATGC‑TCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCCGCAACGcc                                                                                                                                                                                                                                   >  2:388511/1‑251 (MQ=255)
                  acaaatatgttGCGTCGCCAACCCCAGGGGGATGCCTAATAA‑GAAGACCCCCGAAATTAAGGGGCGGAACAAATTTCCCCGATGCCAGGCAACCCCCAACGCCAAGAAAAAGCCGATGGGGGCAAAATCGCCAAACGGCAGCACCAGATTTCCCGGCCCCCAAACAGCAATTAAAATGGGGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                                                                                      <  2:363512/242‑1 (MQ=255)
                              cGTCGCGATCCACAGGGTGATGC‑TCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                          >  2:3052/1‑251 (MQ=255)
                                gcgcgACCCACGGGGGGATGC‑TCATAATGAGGACCCCCGACATGAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCAAGAAGAAGCCGATGGTGGCAAGATCGCAAACCGGCAGCCCCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCGTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCATGAACTCCTGa                                                                                                                                                                                                        <  2:280713/250‑1 (MQ=255)
                                                   gC‑TCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                                                      >  1:220695/1‑250 (MQ=255)
                                                   gC‑TCATAAGGATAACCCCCGAGATTAGGGTGGGGAACAGATTCCCACGAGCCCCGGCGCCCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGACCGCCAAAGGGAAGCACCTGATTACCCGGCACACAAACTGCAATTAAAATGGTGAGGGGGAAAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATTTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTTACCGCCGAACTTCGCCTGt                                                                                                                                                                                      <  2:404691/250‑1 (MQ=255)
                                                    ccttaTAATGATGACCCCCGGGATTAAGGTGCGGAACAGATTCACACTATGCCCGCCGACAGCCATCGCTACAAAGAAGCCTATGGTGGAAAGATCGCCAAACGGCGGCCCCTGAGTACCCGGCACACAAATAGCAAATAAAATGTTGAGGGGTATAAAAAACAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                                                     <  1:431165/247‑1 (MQ=255)
                                                      tCATAATGATGACACCCGAGATTAAGATGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑GCCGAATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCACCTGTaa                                                                                                                                                                                    >  1:38495/1‑250 (MQ=255)
                                                            tgatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                              >  1:374979/1‑250 (MQ=255)
                                                            tgatgaCACCCGAGATTAAGGGGCGGAACAGATTTACAAGAAGCACGGAGACCGCCAAAGACACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGAAGCACCTGAGTAACCGGCACAAAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGAAA‑‑ACACCGCCGTATGTAACAGAAAAAACGCCGGATCAAGGCCAATCCGGCAATCCTGACCGCAGAAAGTCGCCGGTAAACGACt                                                                                                                                                                              >  2:399210/1‑250 (MQ=255)
                                                            tgatgaCACCAGAGATTAAGGTGCGGAACAGATTTCCACGATGCACAGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACAGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGTGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                              >  1:41227/1‑250 (MQ=255)
                                                                  cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGcc                                                                                                                                                                       >  1:414181/1‑251 (MQ=255)
                                                                    cccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTg                                                                                                                                                                     >  2:112097/1‑251 (MQ=255)
                                                                       gagaTTAAGGTGCGGAACAGCTTTCCACGATGCACGGCGACCGCCATCGACCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCAACTGATTACCCGGCCCACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGCCTACGAACCTGCtt                                                                                                                                                                  >  1:146251/1‑251 (MQ=255)
                                                                       gaaaTTAAGATACGGAACAGATTTCCCCGCTCCACGGCGCCCGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAACGGCAGAACCTGCTTACCAGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAACAATCAGGCTTGACGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGAACGCCGAACTTCGCCTGTAAACGACTAAGCGCCTGCtt                                                                                                                                                                  >  2:291298/4‑251 (MQ=255)
                                                                        agtttagggTGCGGAACAGATTCCCACGATGCCCGGCGCCCCCCACCGCCACGAAGAAGCCGAGGGGGGCAAGACCGCCAAACGGCAGCACCTGATTCCCCGGCCCACAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                                                                                                                                 <  2:374979/244‑1 (MQ=255)
                                                                                gTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                                                                                                                                         >  2:110370/1‑251 (MQ=255)
                                                                                   cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                                                                                                                                      >  2:346669/1‑251 (MQ=255)
                                                                                       aCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                                                  <  1:388511/251‑1 (MQ=255)
                                                                                           ccgatcgcgaggCACGGCGACCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCTATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAATGCCAATCAGGATCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                              <  1:226011/240‑1 (MQ=255)
                                                                                                  cGATGAACGGCGACCGCCATCGCCCAGAAGAAGCCGATGGTGGCCAGATCGCCAAACGGCAGCACATGATTACCCGGCACACAAACAGCAATTAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                              >  2:226011/1‑244 (MQ=255)
                                                                                                    aTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                                                                                                                                     >  1:197742/1‑251 (MQ=255)
                                                                                                             gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa                                                                                                                            <  1:122536/251‑1 (MQ=255)
                                                                                                                            gaagaaGCCGATGGGGGAAAGATCCCCAAACGGCGGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagca                                                                                                             <  1:112570/251‑1 (MQ=255)
                                                                                                                                  gCCGGTGGTGGAAAGAGCGCCAAAAGGCACTACTCGATTACCCGGCACAAAAAATACAATAAAAATGGGGGGGGGGAAAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGGGGCATTAGCAGCATCa                                                                                                       <  2:135673/251‑1 (MQ=255)
                                                                                                                                                      cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                   >  1:198715/1‑251 (MQ=255)
                                                                                                                                                       aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                   >  2:127864/1‑250 (MQ=255)
                                                                                                                                                                             gccacaaaaaaccaaaTTAAAAGGTTGAGTGGAAAAAAAACCAGGCTGGCAGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCTTGGCATTAGCATCATCACTGCCGCCGTTTGTACCGCCAGCTGCAATACACCTTTGACATc                                                            <  2:414181/237‑1 (MQ=255)
                                                                                                                                                                                                aaaTGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                                         >  1:114435/1‑251 (MQ=255)
                                                                                                                                                                                                           gggATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGccc                              >  1:207701/1‑251 (MQ=255)
                                                                                                                                                                                                           gggATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGAAACACGCGAGGATGCCGATAATCAGccc                              >  2:269760/1‑251 (MQ=255)
                                                                                                                                                                                                                       aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGAAACCCGCGAGGATGCCGATAATCAGCCCCATCACAACAcc                  >  2:179413/1‑251 (MQ=255)
                                                                                                                                                                                                                         gCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGAAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAccc                 >  1:279578/1‑250 (MQ=255)
                                                                                                                                                                                                                              cgccgatcccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAGCTCCGGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCGGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc               <  1:110370/238‑1 (MQ=255)
                                                                                                                                                                                                                                         ccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg  <  1:173824/251‑1 (MQ=255)
                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                                         
CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGC‑TCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173131‑2173613

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