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.93.8% 197.4 / 6.7 48intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (29/16);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 5.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01
*NC_0009132,173,3640C.93.8% 198.3 / 3.5 48intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (29/16);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 5.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173131‑2173612
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                        
cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                        >  1:485686/1‑250 (MQ=255)
 gccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAgcagca                                                                                                                                                                                                                                       >  2:370261/1‑250 (MQ=255)
 gccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCAACAGGGTGATGCTCATAATGATGACACACGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATGCCACC‑‑GCCGTATGTCCCAgcagca                                                                                                                                                                                                                                       >  1:285442/1‑250 (MQ=255)
        gggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                                                                                                                                               <  1:415447/251‑1 (MQ=255)
         ggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTAATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                                              <  1:560957/251‑1 (MQ=255)
          gtgtGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                             >  2:461096/1‑251 (MQ=255)
                   ccAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                                                    >  2:63728/1‑251 (MQ=255)
                          tttGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCTAATCAGGAAc                                                                                                                                                                                                             >  2:60532/1‑251 (MQ=255)
                              cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCt                                                                                                                                                                                                         >  2:401537/1‑251 (MQ=255)
                                 cgcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                                                                                                                                      <  1:277970/251‑1 (MQ=255)
                                     aTCCACAGGTTGATGCTCATAATGATGACCCCCGAGATTAAGGGGCGGAAAAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                                                  <  2:271044/251‑1 (MQ=255)
                                         acaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                                                                                              >  2:468823/1‑251 (MQ=255)
                                             ggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                                          >  1:487442/1‑251 (MQ=255)
                                                gATGCTCATAATGATGACACACGAGATTACGGTGCGGAACAGATTTCCACGATGCACGGCGACCGACATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCt                                                                                                                                                                                       >  2:190915/1‑251 (MQ=255)
                                                                  acCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGctt                                                                                                                                                                     >  1:69398/1‑249 (MQ=255)
                                                                    ccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGc                                                                                                                                                                   >  2:173472/1‑251 (MQ=255)
                                                                    ccGAGATTAAGGTGCGGAACAGATTTCAACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATTGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGAACGCCGAACTTCGCCTGTAAAGGACTACGAGCTTGc                                                                                                                                                                   >  2:173454/1‑251 (MQ=255)
                                                                      gagaTTAAGGTGCGGAACAGATTTCCCCGAGGCAGGGCGACCGCCATCGCCACAAAGAAGCCGAGGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGT‑TGAGGGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                                 <  2:24883/251‑1 (MQ=255)
                                                                       agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCTTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTa                                                                                                                                                                <  2:167765/251‑1 (MQ=255)
                                                                       agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTa                                                                                                                                                                <  1:28098/251‑1 (MQ=255)
                                                                                         gATTTCCACGATGCACGGCGACCGCCATCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GTCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                               >  1:331720/1‑250 (MQ=255)
                                                                                               cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                        >  2:556928/1‑251 (MQ=255)
                                                                                               cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGAGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                        >  1:457987/1‑251 (MQ=255)
                                                                                                      cACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCTAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                                                                                                                                 <  1:60532/251‑1 (MQ=255)
                                                                                                       cgggcacccgccaccgccaCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAAAAAAAGAAATAAAAATGG‑TGAGTGGGAAAAAAACCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                <  2:578823/244‑1 (MQ=255)
                                                                                                         ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTACACCATCCATGATGGGtttc                                                                                                                               >  2:253320/1‑249 (MQ=255)
                                                                                                              cccccaccGCCACAAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAAGGG‑TGAGTGGGAAAAAAATCAGGCTTG‑CCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                                                         <  2:450998/244‑1 (MQ=255)
                                                                                                                           gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCCGCACCTGATTACCCGGCACACAAACAGCAATTACAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCAACGACTCCGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcc                                                                                                            >  1:270868/1‑250 (MQ=255)
                                                                                                                                            cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                           >  2:157069/1‑251 (MQ=255)
                                                                                                                                            cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                           >  2:447368/1‑251 (MQ=255)
                                                                                                                                                  cGCCAAAGGGCAGCACCTGATTACCCGCCACACAAACAGCAATAAAATTGT‑TGAGTGGGATAAAAACCAGGCTTG‑CCGTTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATTATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTGTTTAcc                                                                                     <  1:190915/251‑1 (MQ=255)
                                                                                                                                                       aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                  >  1:553579/1‑249 (MQ=255)
                                                                                                                                                                 cctgattcccgggcaacaaaaCAGCAATAAAAAGGG‑GGAGTGGAAAAAAAAACAGGCTTG‑CCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATaca                                                                      <  2:123884/233‑1 (MQ=255)
                                                                                                                                                                  cTGATTACCCGGCACACAAACAGCAATTAACATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                                     >  2:239587/1‑251 (MQ=255)
                                                                                                                                                                      ttccccGGCACAAAAACAGCAAT‑AAAATGGGTGGGTGGAAAAAAAATC‑GGCTGGGCCGATCCCACC‑‑GCCGTATTGCCCAGCAGCAACCCCGGAGCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCGCTGCCGCTGTTTTTGCCGCCAGCTGCAATACACCtt                                                                  <  2:547968/247‑1 (MQ=255)
                                                                                                                                                                            ggCAAACAAACAGCAATAAAAAGGG‑GGAGGGGGAAAAAAACCAGCCTTG‑CCAATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                                           <  2:3220/251‑1 (MQ=255)
                                                                                                                                                                              cacacaAACAGCAATTAAAATGG‑TGAGTGGGATTAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                         >  1:151803/1‑251 (MQ=255)
                                                                                                                                                                                cacaGACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                                        >  2:35479/1‑250 (MQ=255)
                                                                                                                                                                                cacaAACAGCAATTAAAATGC‑TGAGTGGGAAAAAAATCAGGCTTG‑CCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTaa                                                       <  1:277085/251‑1 (MQ=255)
                                                                                                                                                                                    aaaaGCAATTAAAATGG‑GGAGTGTGATAAAAATCAGGTTTA‑CCGAAACCACC‑‑GCCGTAGGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAGCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCTCCAGCTGCAATACACCTTTGACATCGTAACCcg                                                   <  1:40637/248‑1 (MQ=255)
                                                                                                                                                                                      cAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                                                 >  1:187094/1‑251 (MQ=255)
                                                                                                                                                                                           ataaaaaGGG‑TGAGGGGAATAAAAATCGGGCTGG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                                            <  1:180084/248‑1 (MQ=255)
                                                                                                                                                                                                               tAAAAATCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCCATACACCTTTGACATCGTAACCCGCGAGGATGCCTATAATCAGCCCCATc                         >  1:427235/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaaaacagggctG‑CCGGTACC‑CCCCGCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATca                        <  2:528120/240‑1 (MQ=255)
                                                                                                                                                                                                                     tCAGGCTTG‑CCGATACCACC‑‑GCCGTATGTTACAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGACCTTCGCCTGTCAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACaaa                    >  2:542234/1‑250 (MQ=255)
                                                                                                                                                                                                                      cAGGCTTG‑CCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAc                  >  1:204550/1‑251 (MQ=255)
                                                                                                                                                                                                                        ccaggctggcagacacccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACaaaa                   <  2:395316/235‑1 (MQ=255)
                                                                                                                                                                                                                            tG‑CCGATACCATC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATAAGGAACTCCTGACCGACGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGACAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCACATCACAAAACCCaagg            >  2:488244/1‑248 (MQ=255)
                                                                                                                                                                                                                               ccGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGAGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTg          <  1:370261/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gatacccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc        >  1:383854/6‑251 (MQ=255)
                                                                                                                                                                                                                                       ccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCTCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGc    <  1:570474/248‑1 (MQ=255)
                                                                                                                                                                                                                                         cccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc  <  1:461096/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                        
CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTG‑CCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTC  >  NC_000913/2173131‑2173612

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