Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,311,002 2 bp→GT intergenic (‑22/+104) alsA ← / ← alsB D‑allose ABC transporter ATP binding subunit/D‑allose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,311,0020TG96.8% 87.8 / ‑7.3 31intergenic (‑22/+105)alsA/alsBD‑allose ABC transporter ATP binding subunit/D‑allose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (20/10);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,311,0030AT100.0% 99.1 / NA 31intergenic (‑23/+104)alsA/alsBD‑allose ABC transporter ATP binding subunit/D‑allose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (21/10);  total (21/10)

AATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCGGATAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACC  >  NC_000913/4310780‑4311242
                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                               
aaTGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCGGATAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACAt                                                                                                                                                                                                                        <  1:148451/249‑1 (MQ=255)
          atggtactggtGCCTTTGGTCGGCTCATGTATTCCGGATAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTACCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGCAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCt                                                                                                                                                                                                             >  2:319569/1‑250 (MQ=255)
                                                     cTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTCGTAATGCCTGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTAAGGTTATTTTGTTGCCCGGATAAGACACTTCATTAGc                                                                                                                                                                  >  1:326474/1‑250 (MQ=255)
                                                               cGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCATGTCTACGGTCAGGTTATTTTGTTGGCCGGATAAGTCACTTCATTAGCATCATCATCCg                                                                                                                                                       <  2:152949/251‑1 (MQ=255)
                                                                    aTTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGcatc                                                                                                                                                  >  1:171288/1‑251 (MQ=255)
                                                                        aCCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCAtt                                                                                                                                               <  1:321904/250‑1 (MQ=255)
                                                                         cccGCGCCATTTTCTCCTAGTAATGCATGTATTTCTCCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCTATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCAACATTCa                                                                                                                                             <  1:93650/251‑1 (MQ=255)
                                                                                cATTTTATCCTAGTCAGGCATGTATTTCACCACGATAAACCGATAACTTAACCGACTTTAATGAGTGAACCGGACCAAAGGACTTGCCGATTCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAGATTACGTCTTATTGACTGATGCCTACATTCTTCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGcatcatcacacgggctacttcagttttgt                                                                                                                                      >  2:139175/1‑230 (MQ=255)
                                                                                   tttCTCCTAGTAAAGCATGTATTTCACCAGGATAAACCGTTAACTTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTgt                                                                                                                                   >  2:359212/1‑250 (MQ=255)
                                                                                               aTGCATGTATTTCACCAGGATAAACCGTTAAATGAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGAtt                                                                                                                       <  1:107857/251‑1 (MQ=255)
                                                                                                     gTATTTCACAAGGATCAACCGTTAAATTACCGGCCTTTAATACGTGAACCGGACCAAAGAACATGCCGAACCCCGCCAACGATCTCTATGGCGCCACCATAAATAACCTCTCTGCACGACAGTTCGTATTATTCACTGATGGCTACATTTATCAGATCTACGTTCAGGTTATTTTGTTGGACGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTAACCAGGATTGAATCg                                                                                                                 >  2:51931/1‑251 (MQ=255)
                                                                                                            aCCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGt                                                                                                                                           <  1:220068/218‑1 (MQ=255)
                                                                                                            aCCAGGATAAACCGTTAAATTAAACGACTTTAATGAGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGt                                                                                                                                           >  2:220068/1‑218 (MQ=255)
                                                                                                                ggATAAACCGTTAAATTAACCGAATTTAATGCGTGAACCGGACCAAAGGACTTGCCGATACCCACCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATTAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGGATCATCATCCGACATCATTCAGTTTTATTGAGTGACCATGATTGAATCGCCCAGTTTaa                                                                                                       >  2:333913/1‑250 (MQ=255)
                                                                                                                  aTAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGACAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCCTCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAAtt                                                                                                    >  2:326842/1‑251 (MQ=255)
                                                                                                                      aCCGTTAAATTCACCGACTTTACTGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCTATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGCAAGTTCGTATTATTGCCTGATGGCTAAATTTATCAGGTCTACGATCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCCGGATTGAATCGACCAGGTTAAATTCCgg                                                                                                >  2:46746/1‑251 (MQ=255)
                                                                                                                                  aCCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATcc                                                                                      >  1:301804/1‑249 (MQ=255)
                                                                                                                                         ttAATGCGTGCACCGAACCACAGGACTTGCCGATCCCCGCCATCGAAATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTAAGTTCAAGTTATTTTGTTGGCCGGATAAGACCATTCATTAGCATCATCATCCGGCCTCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGAt                                                                              >  2:137447/1‑250 (MQ=255)
                                                                                                                                          tAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCAATCGCTATATATGACGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATa                                                                             >  1:23462/1‑250 (MQ=255)
                                                                                                                                           aaTGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAAc                                                                           >  1:137501/1‑251 (MQ=255)
                                                                                                                                                       ggtccAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGAGATCCTTGCCGGAttt                                                               <  2:301804/248‑1 (MQ=255)
                                                                                                                                                         aCCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTtct                                                             >  1:305595/1‑251 (MQ=255)
                                                                                                                                                           cAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTtctc                                                            >  2:216093/1‑250 (MQ=255)
                                                                                                                                                                         aTCCCCGCCATCGATATATATAGCATGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGCCACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATACCCTTGCCGGTTTTCTCAGCGTCAACCATCac                                             >  2:24332/1‑250 (MQ=255)
                                                                                                                                                                         aTCCCCGCCATAGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGCCCAGGATTGAATCGACCAGTTTACATTCCGGTGCTTTATACAGCGGGGTTACCTTGCCGGATTTCTCAGCGTCAACCATc                                               >  2:130825/1‑249 (MQ=255)
                                                                                                                                                                                  aTCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTTCGTTCAGGTTATTTTGTTGGCCGGCTAAGACACTTCGTTAACATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAAACTTGCAGGATTTCTCAGCGTCAACCATCAGCTTCAGAcc                                    >  1:118553/1‑251 (MQ=255)
                                                                                                                                                                                                           ataaCCTCTCTGCACGAAAGTTCGTATCATTGCATGATGGCAACATTTATCAGGTATACGTTCAGGTTATTTTGTTGGCCGGATACGACACTTCATTAGCATCATCATCCGTCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGCCCAGTTTAAGTTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGttt           >  2:16772/1‑250 (MQ=255)
                                                                                                                                                                                                              aCCTCTCTGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTgg        <  1:183041/251‑1 (MQ=255)
                                                                                                                                                                                                                    ctGCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCa     <  1:216093/248‑1 (MQ=255)
                                                                                                                                                                                                                     tGCACGAAAGTTCGTATTATTGCCTGATGGATACATTTATCAGGTATACGTTAAGGCTTTTTTGTAGGACGGATAAGACACTTCATTAGCATCACCCACCGGCAAAATTCAGTTTTATTGAGTGAACAGAATTGAATCGACCAGTTTAAATTACGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATCTCCGCCGGGTTCTGGGCAAc   >  1:160971/1‑249 (MQ=255)
                                                                                                                                                                                                                      gCACGAAAGTTCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAAcc  <  2:171288/249‑1 (MQ=255)
                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                               
AATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCGGATAAAACTTTCATTAGCGTGGATTTACCCGCGCCATTTTCTCCTAGTAATGCATGTATTTCACCAGGATAAACCGTTAAATTAACCGACTTTAATGCGTGAACCGGACCAAAGGACTTGCCGATCCCCGCCATCGATATATATGGCGTGGCCATAAATAACCTCTCTGCACGAAATATCGTATTATTGCCTGATGGCTACATTTATCAGGTCTACGTTCAGGTTATTTTGTTGGCCGGATAAGACACTTCATTAGCATCATCATCCGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACC  >  NC_000913/4310780‑4311242

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