Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,769,478 (G)7→8 coding (405/1113 nt) ydiK → putative transporter YdiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,4711.G95.2% 69.6 / ‑3.2 21coding (398/1113 nt)ydiKputative transporter YdiK
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (8/12);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCC  >  NC_000913/1769347‑1769599
                                                                                                                             |                                                                                                                                
ggCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGa                                                                                                                    >  2:94835/1‑140 (MQ=255)
ggCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGa                                                                                                                    <  2:8435/140‑1 (MQ=255)
  cAGCGGCCCACTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATc                                                                                                                  >  1:71439/1‑140 (MQ=255)
             tAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGt                                                                                                       <  1:390000/140‑1 (MQ=255)
              aaTTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                      <  2:4490/140‑1 (MQ=255)
                       cATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTtat                                                                                             <  2:71439/140‑1 (MQ=255)
                         tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTatat                                                                                           >  2:207143/1‑140 (MQ=255)
                                               aCCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTc                                                                     >  1:198155/1‑140 (MQ=255)
                                                     tCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTggg                                                               >  1:34373/1‑140 (MQ=255)
                                                         gCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGGAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                           <  2:356875/140‑1 (MQ=255)
                                                                       cATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCggg                                             <  1:409952/140‑1 (MQ=255)
                                                                        aTTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGGCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGGATATCGGgcg                                           >  1:34845/1‑140 (MQ=255)
                                                                         ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGCGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgcg                                           >  2:125412/1‑140 (MQ=255)
                                                                         ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgcg                                           <  2:198155/140‑1 (MQ=255)
                                                                           ccGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                         <  2:34373/140‑1 (MQ=255)
                                                                           ccGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                         <  1:125412/140‑1 (MQ=255)
                                                                           ccGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                         <  1:63724/140‑1 (MQ=255)
                                                                                           aGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTgc                         >  2:17678/1‑140 (MQ=255)
                                                                                                                 cTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg     >  2:187451/1‑138 (MQ=255)
                                                                                                                 cTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg     <  1:187451/138‑1 (MQ=255)
                                                                                                                  ttGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACATGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGccc  <  1:207143/140‑1 (MQ=255)
                                                                                                                             |                                                                                                                                
GGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCC  >  NC_000913/1769347‑1769599

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