Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,789,643 C→T P11S (CCT→TCT)  ydiE → PF10636 family protein YdiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,789,6430CT92.3% 73.5 / ‑0.4 26P11S (CCT→TCT) ydiEPF10636 family protein YdiE
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (16/8);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

CATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGT  >  NC_000913/1789400‑1789866
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                               
cATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAAc                                                                                                                                                                                                                          <  1:241199/251‑1 (MQ=255)
     ccGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATg                                                                                                                                                                                                                     <  2:175442/251‑1 (MQ=255)
                                   aaCGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCt                                                                                                                                                                                        >  1:257748/1‑250 (MQ=255)
                                          ggTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCg                                                                                                                                                                                >  2:343557/1‑251 (MQ=255)
                                                    cTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTccagc                                                                                                                                                                      >  2:254937/1‑251 (MQ=255)
                                                                                 aaTGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCATAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCTAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACt                                                                                                                                         <  1:380476/251‑1 (MQ=255)
                                                                                               ccTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATg                                                                                                                           <  2:118040/251‑1 (MQ=255)
                                                                                                      gAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCa                                                                                                                    <  1:64078/251‑1 (MQ=255)
                                                                                                                   tGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGaa                                                                                                                                                                              <  1:233677/180‑1 (MQ=255)
                                                                                                                   tGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGaa                                                                                                                                                                              >  2:233677/1‑180 (MQ=255)
                                                                                                                   tGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTc                                                                                                       >  2:365551/1‑251 (MQ=255)
                                                                                                                       cAACTTGTTGTCTTATCAACACATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGt                                                                                                    >  2:170853/1‑250 (MQ=255)
                                                                                                                                   ttctCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGCTGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTg                                                                                       >  1:111084/1‑251 (MQ=255)
                                                                                                                                      tcAACAAATTACTACTTGCTCAGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTACCCGCAGCCTATTCGGCGAATTTCCAGCCAGAAACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGCCGGGCAAGAATATCTGCTCCGTAAAACTAAGGCTgcc                                                                                     >  2:411642/1‑248 (MQ=255)
                                                                                                                                                                          ttGATCATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAAAGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTAGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGATGCTGATGACCAAGTAGCCTTTAACTCGatc                                                  >  2:396443/1‑247 (MQ=255)
                                                                                                                                                                             ataataAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGag                                                   >  2:323354/1‑245 (MQ=255)
                                                                                                                                                                             ataataAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGag                                                   <  1:323354/245‑1 (MQ=255)
                                                                                                                                                                                   aGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTcc                                       >  1:40819/1‑251 (MQ=255)
                                                                                                                                                                                             ttATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAc                             >  1:374197/1‑251 (MQ=255)
                                                                                                                                                                                                tatCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAcacc                          >  2:192330/1‑248 (MQ=255)
                                                                                                                                                                                                        attattaATTTTTATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTACAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGt                  >  1:391515/1‑251 (MQ=255)
                                                                                                                                                                                                                  tttATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGATGGCAAGCAGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTg         >  1:268707/1‑250 (MQ=255)
                                                                                                                                                                                                                    tATGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGAc       <  2:257748/250‑1 (MQ=255)
                                                                                                                                                                                                                      tGCGTTATACGGATAGCAGAAAACTCACGTCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACAt     >  1:203623/1‑250 (MQ=255)
                                                                                                                                                                                                                         gTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGt  <  1:340709/250‑1 (MQ=255)
                                                                                                                                                                                                                         gTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGt  >  2:340709/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                               
CATTACCGACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGT  >  NC_000913/1789400‑1789866

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