Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,721,136 C→T L38F (CTT→TTT)  ydhI → DUF1656 domain‑containing protein YdhI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,721,1360CT94.4% 40.7 / ‑6.6 18L38F (CTT→TTT) ydhIDUF1656 domain‑containing protein YdhI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/12);  minor base A (1/0);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGC  >  NC_000913/1720964‑1721354
                                                                                                                                                                            |                                                                                                                                                                                                                          
cGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCTGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTg                                                                                                                                               <  1:400309/250‑1 (MQ=255)
  cACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGttt                                                                                                                                            >  1:34712/1‑251 (MQ=255)
    cTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTacac                                                                                                                                                                                                                 >  1:143394/1‑180 (MQ=255)
          naaTATTTATTACTCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGc                                                                                                                                                                                       >  1:131350/2‑200 (MQ=255)
                           aCGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGATTGGATTCGTCATCTGGTTTGTCGTa                                                                                                                                                                                                                    <  1:5515/154‑1 (MQ=255)
                                       cGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTg                                                                                                                                           <  1:325254/215‑1 (MQ=255)
                                                                 aGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGNTGGNNNNNCg                                                                                                                                                                                        >  2:305836/1‑144 (MQ=255)
                                                                 aGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTggn                                                                                                                  <  1:305836/214‑2 (MQ=255)
                                                                  ggttATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGAAAGACCGAGTCCTTTACTCCGCTTTTCAACGCATTCGCGTTAGGATTCGTCATCTGGATTGTCGTACACCGACTGATTCGTGGCTGGATCTACGCCGGTGAAATCTGGCATCACTTGTTAATGGATTTATCGATGTTTGCGATTTGCGTTTGACTTGCTCTGGCAATACTGATTGTGTGGTAACTATGTCAATTAAAACAAt                                                                            >  1:457238/3‑251 (MQ=255)
                                                                         cTCAAGGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCGCCGTTGTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTTTCGTACGCCGCTCGCTGCTTTGATGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGttt                                                                                                                                            <  1:474825/180‑1 (MQ=255)
                                                                           cAATGCAANANNNTNNCCCTTACAAGATCTGGTGTTCAGTGCGTCCGTCTACTTTCCTCCNTTTTTCAAAGCATTCNCGNTTGGNTTCNTCATCTGGTTTGt                                                                                                                                                                                                                        >  2:358223/1‑102 (MQ=255)
                                                                           cAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCAGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGt                                                                                                                                                                                                                        <  1:358223/102‑1 (MQ=255)
                                                                                                                    cgtccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTc                                                                                                                                                                                                       <  1:119890/78‑1 (MQ=255)
                                                                                                                      tccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTgg                                                                                                                                                                                               <  1:407652/84‑1 (MQ=255)
                                                                                                                      tccgtcTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTcaan                                                             <  1:57738/214‑2 (MQ=255)
                                                                                                                                         ttttCAAAGCATTCGCGTTTGGATTCGTCATCTGGTTTGTCGTACTCCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATg                                                                                                                                                           <  1:71684/101‑1 (MQ=255)
                                                                                                                                                                     cATCTGGTTTGTCGTACACCGCCTGCTTCGTGTCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTTTAGCGGTAGTTGCTGTGCTTGCCGGATGGTGGCTGTGGAATTATTACATGc  <  1:195928/226‑1 (MQ=255)
                                                                                                                                                                       tCTGGTTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTggn                                                                                                                                                                         <  1:94126/57‑2 (MQ=255)
                                                                                                                                                                            |                                                                                                                                                                                                                          
CGCACTTAAAGAATATTTATTAATCTAACGCAATATATTCGGTCGTAAAGGAATCTACTTTGTGAAGTTTATGCTCAATGCAACAGGATTGCCCTTACAAGATCTGGTGTTCGGTGCGTCCGTCTACTTTCCTCCGTTTTTCAAAGCATTCGCGTTTGGATTCGTCATCTGGCTTGTCGTACACCGCCTGCTTCGTGGCTGGATCTACGCCGGTGACATCTGGCATCCCTTGTTAATGGATTTATCGCTGTTTGCGATTTGCGTTTGCCTTGCTCTGGCAATACTGATTGCGTGGTAACTATGTCAATTAAAACAATTAAGTATTTCTCAACAATCATTGTAGCGGTAGTTGCGGTTCTTGCCGGATGGTGGCTGTGGAATTATTACATGC  >  NC_000913/1720964‑1721354

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