Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 518,989 C→T R54H (CGC→CAC)  ybbO ← NADP(+)‑dependent aldehyde reductase YbbO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913518,9890CT95.2% 52.1 / ‑6.0 21R54H (CGC→CAC) ybbONADP(+)‑dependent aldehyde reductase YbbO
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (8/12);  total (8/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

AACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAATTAAGACCGATTTTTGCATAACTTTACCTGTCAGGATCTCCGTTGCTTTATGAGTCATGATTTACTAAAGGCTGCAACTGCTTCGCCATCCAGTCGGCAATAAACGGCTGGGCGTCGCGGTTGGGATGAATACCGTCATCCTG  >  NC_000913/518749‑519234
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
aaCATCGCGGGTAACAGGCGCAGGGTGAGCTGTGGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAGGGGGGCCATACATGCCGAAGCCGGCATGGTTAAATATCCCATACAGACAATTATCGGTCATGGCGATCACCTCGTCGGCTGCTGTATCAACTCTTTCTGTTGAATCCAGATCGATAAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAAcatc                                                                                                                                                                                                                                             <  1:421669/251‑1 (MQ=255)
                  cgcATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGAGGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCGCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAAcc                                                                                                                                                                                                                           <  1:169152/251‑1 (MQ=255)
                              gggTGGGCGCCGAAAAAGTGGGAGGAAAACTGCTGTTCCATCTGTGCAGGGCGGATGGTGTAAAGGGGGCCATACATGCCGAAGCCGGCAGGGGTAAAGATCCCATAAAGACAATTAGCGGCCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTGTCGGGTGAAGCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTGCACGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATg                                                                                                                                                                                                               <  1:107096/250‑1 (MQ=255)
                                  gtgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGACAGGGGGCCATACATGCAGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCCGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGaaa                                                                                                                                                                                                            >  1:62510/1‑250 (MQ=255)
                                   tgCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGCGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAAc                                                                                                                                                                                                          >  1:186105/1‑251 (MQ=255)
                                       ccGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGACGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTg                                                                                                                                                                                                      >  1:10568/1‑251 (MQ=255)
                                                    cGGAAAACTGCTGTTCCCTCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCg                                                                                                                                                                                         >  1:354145/1‑251 (MQ=255)
                                                                   ccATCTGCGCACGGCTGATGGTGGAAAGGGGCCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTGCTGGTGAATCCAGATCGATCAACACGCCGTTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTcn                                                                                                                                                                                         <  1:53325/236‑2 (MQ=255)
                                                                    cATCTGCGCACGGCGGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTTTCGTTCAGGGCGATCACCTCGTCTGTTGCTCGATGAACACTTTCTGGTGAGGCCTGATCCAGTAACACGCGTGTTAATCCCTGGCTTTTCATGCTGTCGACAGCACCGGGTTTCCGGCAACCTGCCAGCATGTATAAACCCTTGCGTGTTAATTCGAGCGCGCTTTCCAGGc                                                                                                                                                                         <  1:418463/251‑1 (MQ=255)
                                                                         gcgcACGGCGGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTGGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGcc                                                                                                                                                                        <  1:186348/247‑1 (MQ=255)
                                                                                                   ccATACATGCCGAATCCGGCATTGTTAACGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGCTCAACACTTTCTGGTGAATCCAGCTCGATCAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTACTTCGCGCGCGCTTTCCAGACCAATTCCACTGGAACATCCGGTAATTAAGAc                                                                                                                                          >  1:163863/1‑251 (MQ=255)
                                                                                                   ccATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTTCGCGATCAACACGTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCGATTCTGTTCATGTGCTCATCATCATCCGGTTTcc                                                                                                                                                                                                                                  <  1:374221/163‑1 (MQ=255)
                                                                                                                                     cATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAgg                                                                                                                                                                          <  1:126981/185‑1 (MQ=255)
                                                                                                                                         aaGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAAGCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAATTAAGACCGATTTTTGCATAACTTTACCTGTCAGGATCTCCGttn                                                                                                    <  1:279260/250‑2 (MQ=255)
                                                                                                                                                                  aCCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCACCACGCCGGTAAATCCCATGCTGTTCATGTGCTCACCATAATCCGGTTTCCGGCAACCTGCCAGAACATGAAAACCCTGGCGTTTTCATTCGAGCGCGATTTCCAGGCCAATTCCACTGGAACATCCGGTAATTAAGACCGATTTTTGCATAACTTTACCTGTCCGGATCTCCGTTGCTTTATGAGTCATGATTTACTaa                                                                             >  1:284833/1‑249 (MQ=255)
                                                                                                                                                                           gCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTCAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCgg                                                                                                                                                                                                                                >  1:304378/1‑93 (MQ=255)
                                                                                                                                                                              gcgcGATCAACTCTTTCGGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGGTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCTCTTTCCAGGCCAATTCCACTGGAACATCCGGTAATTAAGACCGATTTTTGCATAACTTTAcc                                                                                                                     <  1:57708/197‑1 (MQ=255)
                                                                                                                                                                                                         aGATCGNNCAACACGACGGTAAATCCCATGCTGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGtttt                                                                                                                                                                                               <  2:10568/96‑1 (MQ=255)
                                                                                                                                                                                                                                        tGTTCATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACt                                                                                                                                                               <  2:266993/97‑1 (MQ=255)
                                                                                                                                                                                                                                        ngTTCATGTGCTCAACATCATCCGGGTTCCGGCAACCTGCCAGCACATGAAAAACCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACt                                                                                                                                                               >  1:266993/2‑97 (MQ=255)
                                                                                                                                                                                                                                            cATGTGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGAAACATCCGGTAATTAAGACCGATTTTTGCATAACTTTACCTGTCAGGATCTCCGTTGCTTTATGAGTCATGATTTACTAAAGGCTGCAACTGCTTCGCCATCCAGTCGGCAATAAACGGCTGGGCGTCGCGGTTGGGATGAATACCGTCATCCTg  <  1:197717/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
AACATCGCGGGTAACAGGCGCATGGTGAGCTGGTGTGCGCCGAAAAAGTTGGCGGAAAACTGCTGTTCCATCTGCGCACGGCTGATGGTGGAAAGGGGGCCATACATGCCGAATCCGGCATTGTTAAAGATCCCATACAGACAATTATCGGTCAGGGCGATCACCTCGTCGGCTGCGCGATCAACACTTTCTGGTGAATCCAGATCGATCAACACGCCGGTAAATCCCATGCTGTTCATGCGCTCAACATCATCCGGTTTCCGGCAACCTGCCAGCACATGAAAACCCTGGCGTTTTAATTCGAGCGCGCTTTCCAGGCCAATTCCACTGGAACATCCGGTAATTAAGACCGATTTTTGCATAACTTTACCTGTCAGGATCTCCGTTGCTTTATGAGTCATGATTTACTAAAGGCTGCAACTGCTTCGCCATCCAGTCGGCAATAAACGGCTGGGCGTCGCGGTTGGGATGAATACCGTCATCCTG  >  NC_000913/518749‑519234

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