Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,855,920 A→G L18L (TTG→CTG)  ydjG ← NADH‑dependent methylglyoxal reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,855,9200AG96.3% 82.1 / NA 27L18L (TTG→CTG) ydjGNADH‑dependent methylglyoxal reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (13/13);  minor base C (1/0);  total (14/13)
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

CCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAACCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTACGCCGGTGGTGGCACCG  >  NC_000913/1855799‑1856047
                                                                                                                         |                                                                                                                               
ccTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTcgcg                                                                                                                      <  1:8820/133‑1 (MQ=255)
 cTGGCGCAGTGTCAACCAGATTAACGCCCCAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGa                                                                                                                     <  2:146575/133‑1 (MQ=255)
    gcgcAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCCATGGCCCATGTCCCCACCCCCCGTCGCGAAAg                                                                                                                  >  2:239902/1‑133 (MQ=255)
    gcgcAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAg                                                                                                                  >  1:127698/1‑133 (MQ=255)
    gcgcAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAAGGGCCCATTTCCACATCCCCATCCGCGAAAg                                                                                                                  >  1:155107/1‑133 (MQ=255)
    gcgcAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAAGGGCCCATATCCGCAGCCCCATGCGCGAAAg                                                                                                                  >  1:133476/1‑133 (MQ=255)
         gTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTaa                                                                                                             >  2:352504/1‑133 (MQ=255)
                tcaGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGt                                                                                                      <  2:419141/133‑1 (MQ=255)
                    aTTAATGCCACAACGATGGGCTTCGAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTtgtg                                                                                                  >  1:345243/1‑133 (MQ=255)
                    aTTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGGCCCATTCTCGACAGAGTAATAGCCGTtgtg                                                                                                  >  1:159227/1‑133 (MQ=255)
                      tAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGcc                                                                                                <  1:32306/133‑1 (MQ=255)
                                   aTGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCttt                                                                                   >  1:288495/1‑133 (MQ=255)
                                         ttCAAGAATCTTATCCATACATATTTGCTGATCGAGATCGCTATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCAtt                                                                             <  2:133476/133‑1 (MQ=255)
                                            aaGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGc                                                                          >  1:497580/1‑133 (MQ=255)
                                            aaGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGc                                                                          <  2:288495/133‑1 (MQ=255)
                                                   gaaTCAATACATATTCGCCGATCGCGATTGCCCTTCCATGCCGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCAtt                                                                   <  2:234399/131‑1 (MQ=255)
                                                    tATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTa                                                                  >  2:484013/1‑133 (MQ=255)
                                                         aTACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCtt                                                             <  2:127698/133‑1 (MQ=255)
                                                               atTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTAtt                                                       >  2:275924/1‑133 (MQ=255)
                                                                   gCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTAttcttc                                                   >  2:270189/1‑133 (MQ=255)
                                                                                           cAGGACCGCCGCCAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCt                           <  1:484013/133‑1 (MQ=255)
                                                                                                   ccgccAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTa                   >  1:222922/1‑133 (MQ=255)
                                                                                                   ccgccAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGCTGTGCCTAAAGGAATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCGGCAATTGTTCGGCCGGCTTTCGGTTTTTTa                   >  1:144939/1‑133 (MQ=255)
                                                                                                     gccAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTACg                 <  2:222922/133‑1 (MQ=255)
                                                                                                     gccAATGGCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCATTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTACg                 <  2:497580/133‑1 (MQ=255)
                                                                                                            gCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTAca                 <  1:496599/126‑2 (MQ=255)
                                                                                                            gCCCATGTCCCCAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTAca                 >  2:496599/1‑125 (MQ=255)
                                                                                                                    ccccAGCCGCATTCGTGACAGCATAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTACGCCGGTGGTGGCACCg  <  1:239902/133‑1 (MQ=255)
                                                                                                                    ccccAGCCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTACGCCGGTGGTGGCACCg  <  1:10190/133‑1 (MQ=255)
                                                                                                                         |                                                                                                                               
CCTGGCGCAGTATCAATCAGATTAATGCCACAACGATGGGCTTCAAGAATCGTATCAATACATATTTGCCGATCGAGATCGCCATTCCATGCAGGACCGCCGCCAATGGCCCATGTCCCCAACCCCATTCGCGAAAGCGTAATATCCGTTGTGCCTAAAGGTATCTTTTTCATTTGCCTTCATTATCCTTCGTATTCTTCCAGCAATTGTTCGACCAGCTTTCTGTTTTTTACGCCGGTGGTGGCACCG  >  NC_000913/1855799‑1856047

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