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.87.2% 128.4 / 13.5 39intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base . (25/9);  total (26/13)
Fisher's exact test for biased strand distribution p-value = 3.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01
*NC_0009132,173,3640C.87.2% 127.4 / 14.5 39intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base . (25/9);  total (26/13)
Fisher's exact test for biased strand distribution p-value = 3.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

CCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGC  >  NC_000913/2173133‑2173614
                                                                                                                                                                                                                                       ||                                                                                                                                                                                                                                                          
cgaGTTGGGTGGGCAAACCAATAGTTTGCGGCGCGACCCACAGGGGGATGCTCATAATGATGACACCCGAGATTAAGGTGGGGAACAGATTTCCACGATGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAAc                                                                                                                                                                                                                                       <  2:66003/249‑1 (MQ=255)
  aGTTGGGAGGGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                                                     >  2:187589/1‑251 (MQ=255)
     tGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGAGCGGAACAGATTTCCACGATGCACGCCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGCATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCCGCAACGCCgg                                                                                                                                                                                                                                  >  1:536456/1‑251 (MQ=255)
                     tAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                                                  >  2:371782/1‑251 (MQ=255)
                         ttGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                                                              >  1:455153/1‑251 (MQ=255)
                          tGCGTAGCGATCCACAGGGTGCTGCTCATAATGATGACACCCGAGCTTAAGGCGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGACGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCAGGCACACAAACAGCAATTAAAA‑TGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                                                                                                                                             >  2:37720/1‑251 (MQ=255)
                           gCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACTAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                                                                                                                                             <  1:22898/250‑1 (MQ=255)
                           gCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGCGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTcc                                                                                                                                                                                                            >  2:588935/1‑251 (MQ=255)
                                    tCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGCTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCACGATCGCCACACGGCAGCACCTGATCACCCGGCACACAAACAGCAATTAAAA‑TGGTAAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGAAACGACGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                                   >  1:590203/1‑251 (MQ=255)
                                     ccACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                                   >  1:109156/1‑250 (MQ=255)
                                         aGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCAACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGAACACAAACAGCAATTAAAA‑TGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                                                                                              >  2:22195/1‑251 (MQ=255)
                                               aTGCTCATAATGATGAAACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGACATCGCCACGAAGAAGCCAATGGTGGCAAGATCGCCAAAAGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑CGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATAACGGCCAATCAGGAACTCCTGAACGCCGAACTTCGCCTg                                                                                                                                                                                        >  2:293123/1‑251 (MQ=255)
                                                     aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGCTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCAGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCACGCCGAACTTCGCCTGTAAACg                                                                                                                                                                                  >  2:477738/1‑251 (MQ=255)
                                                                              tGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                                                                                                                                         >  2:191018/1‑251 (MQ=255)
                                                                              gggggAACAGATTTCCACGAGGCCCGGCCACCCCCCACGCCCCGAAGAAGCGGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                                                                                                                                         <  2:371177/248‑1 (MQ=255)
                                                                                gtgaccaaattcccgcGATGCACGGCGACCCCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATGACCCGGCACACAAACATCAATTAAAA‑TGGGGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCATCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCTCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                                                                                                                                       <  1:300549/236‑1 (MQ=255)
                                                                                   aaCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGACAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGACATCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                                                                                                                                     >  2:401968/1‑250 (MQ=255)
                                                                                    aCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGCACt                                                                                                                                                                                               >  2:377278/1‑207 (MQ=255)
                                                                                          gcccACGATGCACGGCGACCGCCAGCGCCACGAAGAAGCCGATGTTGGCAAGATCGCCAAACGGAAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACGTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGTGCGTTAAAccatc                                                                                                                                              <  1:494967/248‑1 (MQ=255)
                                                                                             cACGATGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGAATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCCTGCTTAGCGATTGTCGTTAAACCATCCatg                                                                                                                                          <  1:487345/251‑1 (MQ=255)
                                                                                                       ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACCGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGTGCGTTCAACCATCCATGATGGGTTTa                                                                                                                                 >  1:564006/1‑250 (MQ=255)
                                                                                                            ccGCCCTCGCCACGAAGAAACCGATGGTGGCAAGATCGCCAAACGCCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAACATCAGGCTTGCCGATACCACC‑‑GACGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCAATGATGGGTTTAATCAcc                                                                                                                           >  2:577672/1‑251 (MQ=255)
                                                                                                                     ccACGAAGACGCCGGTGGAGGCAAGGTCTCCAAACGGCAGCACCTCATTCACCGGCACACAACAAGCACTTAAAA‑TCGTGAGTGGGCTAAAAATCAGGCTTGCCGATACCACC‑‑GACGTATGTACAAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCACCGAACTTCGCCTGTCACAGACTACGAGCATGCTTCGCCATGGGCGTTAAACCATCCCTTATGGGTTTAATCACCCGTGGCAtt                                                                                                                  >  2:554681/1‑251 (MQ=255)
                                                                                                                         gaagaaGCCGATGGTGGCAAGATCGCCAAACGGAAGCACCTGATTACCCGGCACACAACCAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAACCTGCTTAGCGCAGGGCGTTAAACAATCCATGATGGGTTTAATCACCCGGGGCATGagc                                                                                                               >  2:485008/1‑250 (MQ=255)
                                                                                                                                 cGATGGGGGAAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCAAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt                                                                                                      <  1:550426/251‑1 (MQ=255)
                                                                                                                                              accgccagccggaagcccttgataccccGGCACCCAAACAGCAATGAAAA‑GGTTGAGGGGGAAAAAAATCGGGCTTGCTGATCCCACC‑‑GCCGTGTGTCCCAGCAGCAACGCCGGATCAGGGCCAATCAGTAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATTATGGGTTTAATCACCCGTGGCATTAGCATCATCACGGCCGCCGTTTTTa                                                                                         <  2:345974/226‑1 (MQ=255)
                                                                                                                                                     agcGGCAGCCCCTGATTCCCGGGCACACAAAAAG‑AATTAAAAAGGGTGAGGGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg                                                                                  <  1:213883/249‑1 (MQ=255)
                                                                                                                                                         gcacccccTGATTCCCCGCCACACAAACAGCAATAAAAA‑TGGTGAGTGGTATAAAAATCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGc                                                                              <  1:485008/245‑1 (MQ=255)
                                                                                                                                                                   aTTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGGGGGAAAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATGAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCtt                                                                    <  1:509821/251‑1 (MQ=255)
                                                                                                                                                                         cGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                              >  2:594491/1‑251 (MQ=255)
                                                                                                                                                                          ggCACAAAAACAGCAATTAAAA‑TGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                                             <  1:37720/251‑1 (MQ=255)
                                                                                                                                                                           gCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                            >  2:185607/1‑251 (MQ=255)
                                                                                                                                                                                   aCAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGAccgccg                                                                                                                                                                                                   >  1:595246/1‑108 (MQ=255)
                                                                                                                                                                                                 tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                        >  1:377960/1‑250 (MQ=255)
                                                                                                                                                                                                 gggTGAGGGGGAAAAAACTCACGCTTGCCGA‑ACC‑CACCGCCGTATGTCCCAGCACCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTATTCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAGCGTAACCCGCGAGGATGCCGAt                                        <  2:417068/249‑1 (MQ=255)
                                                                                                                                                                                                    tGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                    >  1:564882/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                    >  1:385860/1‑249 (MQ=255)
                                                                                                                                                                                                                      ggCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGCTAATCAGCCCCATCACAAAAccc                  >  1:241870/1‑251 (MQ=255)
                                                                                                                                                                                                                         ttGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACCTCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg               >  2:4905/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGc      <  1:477738/248‑1 (MQ=255)
                                                                                                                                                                                                                                     caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCt     >  1:578618/3‑251 (MQ=255)
                                                                                                                                                                                                                                        cgccgTAAATCCCAACAGCAACGCCGGATCAAGGCCAATCAGGGACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACAATCCATGAAGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACATTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAACCCCACGGGGACAGGCACGc  >  2:373219/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ||                                                                                                                                                                                                                                                          
CCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGC  >  NC_000913/2173133‑2173614

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