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.88.0% 183.8 / 16.2 50intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base . (25/19);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 2.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3640C.90.0% 193.3 / 11.3 50intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (25/20);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 5.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGCCAGTTGGGTGTGCAAACCAA‑TAGTTTGCGTCGCGATCCA‑CAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173128‑2173612
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                                        
aGCCGCCAGTTGGGTGTGCAAACCAA‑TAGTTTGCGTCGCGATCCA‑CAGGGCGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                           >  2:554126/1‑250 (MQ=255)
     caaGTTGGGTGTGCAAACAAA‑TAGTTTGCGTCGCGATCCA‑CAGGGGGATGCTCAAAATGATGACCCCCGAGATTAAGATGCGAAACAGATTTCCACGATGCACGGGGACCGCCTGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCAGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAAc                                                                                                                                                                                                                                     <  1:515261/249‑1 (MQ=255)
        gTTGGGTGTGCAAACCAA‑TAGTTTGCGTCGCGATCCA‑CAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCAACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                                  >  1:573395/1‑251 (MQ=255)
              gtgtgaaaCCAAATAGTTGGCGTCGCG‑TCCAACAGGGTGATGCCCATAATGATGACCCCCGAAATTAAGGGGGGGAACAGATTTCCACGATGCACGGCGACAGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                             <  2:286931/246‑1 (MQ=255)
                 gCAAACCAA‑TAGTTTGCGTCGCGATCCA‑CAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGCGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                                         >  2:345175/1‑251 (MQ=255)
                         a‑TAGTTTGCGTCGCGATCCA‑CAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                   >  2:609796/1‑249 (MQ=255)
                         a‑TAGTTTGCGTCGCGATCCA‑CAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                                                                 <  1:132317/251‑1 (MQ=255)
                                tGCGTCGCGATCCA‑CAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACt                                                                                                                                                                                                            <  1:152395/250‑1 (MQ=255)
                                          tCCA‑CAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                                 >  2:522358/1‑251 (MQ=255)
                                           tccccaGGTTGATGCTCATAATGATGACTCCCGAAATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGTTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCATCGCCGGATCAAGGCCAATCAGTAACTCCTGAccgcc                                                                                                                                                                                                  <  1:486242/246‑1 (MQ=255)
                                                ggggtgagacTCATAAGGAGGCCACCCGAGATTAAGGTGCGGAACATAGTCCCCCAATGCACGCCGACCGCCATCGCCCCGAAGAACCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTAATTACCCGGCACACAAACAGCAATGAAAAGGGGGAGG‑GGGATAAAAACCA‑GGCTTG‑CCGAAACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                                                                                            <  2:150974/242‑1 (MQ=255)
                                                    tGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                                                                                                                                        >  1:177782/1‑251 (MQ=255)
                                                      aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                                                                                                                                      >  2:406715/1‑251 (MQ=255)
                                                      aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg                                                                                                                                                                                      >  1:241020/1‑251 (MQ=255)
                                                           aaaaaTGATGACCCCCGAGATTAAGGTGCGGAACAGATTCCCACTATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAGGGCAGCACCTGATTACCCGGCACAAAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaan                                                                                                                                                                                 <  2:269616/248‑2 (MQ=255)
                                                               atgatgACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                             <  1:559744/251‑1 (MQ=255)
                                                                tgatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGTTTACCCGGCACACAAACAGCAATTAACATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCCGAACTTCGTCTGTAAACGACtt                                                                                                                                                                            >  1:357579/1‑250 (MQ=255)
                                                                 gatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCAGGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                                                           <  1:148413/251‑1 (MQ=255)
                                                                 gatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                                                                           >  2:106073/1‑251 (MQ=255)
                                                                    gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                          >  2:568334/1‑249 (MQ=255)
                                                                              aTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                                                                                                                                              >  2:10012/1‑251 (MQ=255)
                                                                                    gTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTCGCGACggg                                                                                                                                                        >  2:101080/1‑251 (MQ=255)
                                                                                        ggAACAGATTTCCACGATCCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGCGGCAAGATCGCCAAACGGCCGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTCGCGATCGGCGtt                                                                                                                                                    >  2:287560/1‑251 (MQ=255)
                                                                                                          gcccGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                                                                                                                                  <  1:522358/248‑1 (MQ=255)
                                                                                                          ccccGGCGCCCGCCATCGCCACGAAGAGGCCGTTGGTGGAAAGATCGCCAAACGGCAGCCCCTGATTCCCCGCCACAAAAACAGCAATAAAAATGGGGAGG‑GGAATAAAAACCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                                                                                                                                  <  2:130461/248‑1 (MQ=255)
                                                                                                            ccGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                <  1:11108/250‑1 (MQ=255)
                                                                                                            aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                >  1:10326/1‑251 (MQ=255)
                                                                                                                       cgtCGCCCCGAAGAAGCGGAGGGGGGCAAACCCGCCAAACGGCAGCACCTGATAACCCGGCACACAAACGGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                                     <  2:83471/249‑1 (MQ=255)
                                                                                                                         ccGCCACGAAGAACCCGATGGTGGCAAGATCGCCAAACGCCAGCACCTGATTCCCCGGACCACAACCAGCAATAAAAATGGT‑AGTTGGAATAAAACTCA‑GCCTTG‑CCG‑TACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTggn                                                                                                                   <  2:277970/250‑2 (MQ=255)
                                                                                                                            ccACGAAGAAGCCGGTGGTGGAAAGACCGCCATACGGCAGCACCTGATTCCCCGGCCCACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GCCTTC‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                <  1:159243/251‑1 (MQ=255)
                                                                                                                                       ccGATGGCGGAAAGCTGGCCAAACGGCACCACCTGATTACCCGGCAAACAAACGGCAATTAAAATGGTGAGT‑GGGATAAAATTCA‑GG‑TTGTCGGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTTCTTAGCGATGGGCGTTAACCCATCCATGATTGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                                                     <  1:245388/251‑1 (MQ=255)
                                                                                                                                           gggtggCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCAAACAAACAGCAATTAAAATGGGGAGT‑GGGATAAAAATCA‑GGCTGG‑CCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                                                 <  1:93209/250‑1 (MQ=255)
                                                                                                                                           gggtggAAAGATCGCCAAACGGCAGCACTTGATTCCCCGGCACACAAACATCAATTAAAATGGTGAGT‑GGGATAAAAACCA‑GGCTTG‑CCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                                                 <  1:211318/250‑1 (MQ=255)
                                                                                                                                               ggCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAACCA‑GGCTTG‑CCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                              <  1:568334/250‑1 (MQ=255)
                                                                                                                                                               ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                             >  1:500409/1‑251 (MQ=255)
                                                                                                                                                                            ccccagcccaccaaccAGCAATTAAATTGGTGGGT‑GGGA‑AAAAAACACGTCTT‑‑CCCATCCCACCC‑GCGGTATTTCCCGGCAGCAACGCCGGACCAAGGCCAAGCAGGAACCCATGACCGCCGAACTTCGCCTGTAAACGACTACTAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                                 <  2:547805/235‑1 (MQ=255)
                                                                                                                                                                                cGGCACACAAACAGCAATTAAAAGGGGGAGG‑GGGATAAAAATCG‑GGCTTG‑CCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                            <  1:345175/251‑1 (MQ=255)
                                                                                                                                                                                 ggCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCACTACACCTTTGACATc                                                           >  2:9064/1‑251 (MQ=255)
                                                                                                                                                                                   cacacaAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                         <  1:406715/251‑1 (MQ=255)
                                                                                                                                                                                      acaAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                                      >  1:142740/1‑251 (MQ=255)
                                                                                                                                                                                              cAATTAAAATGGTGAGA‑GGGATACAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGAATAACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGacg                                              >  2:579054/1‑249 (MQ=255)
                                                                                                                                                                                                 ttAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATg                                           >  1:79050/1‑251 (MQ=255)
                                                                                                                                                                                                   aaaaGGGTGAGT‑GGGATAAAAATCG‑GGCTTG‑CCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                         <  1:609796/251‑1 (MQ=255)
                                                                                                                                                                                                       tggagcgT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGATTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACAGACAGCTGCAATACACCTTTGACATAGTAACCcgc                                                  >  2:132926/7‑238 (MQ=255)
                                                                                                                                                                                                         gTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAAt                                   >  1:459049/1‑251 (MQ=255)
                                                                                                                                                                                                            aGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                              >  1:18927/1‑221 (MQ=255)
                                                                                                                                                                                                                         aTCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACaa                     >  1:524659/1‑250 (MQ=255)
                                                                                                                                                                                                                         aTCA‑GGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGACGCCGCTAATCAGCCCCATCACaa                     >  2:427152/1‑250 (MQ=255)
                                                                                                                                                                                                                           ccgggCTTG‑CCGAT‑CC‑CCCCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACaaaa                   <  2:177782/248‑1 (MQ=255)
                                                                                                                                                                                                                                  tG‑CCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc              <  1:554126/249‑1 (MQ=255)
                                                                                                                                                                                                                                               cccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  <  1:287560/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                                        
AGCCGCCAGTTGGGTGTGCAAACCAA‑TAGTTTGCGTCGCGATCCA‑CAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGT‑GGGATAAAAATCA‑GGCTTG‑CCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173128‑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: ‑