Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,384,808 C→T K305K (AAG→AAA yfbK ← IPR002035/DUF3520 domain‑containing protein YfbK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,384,8080CT96.7% 96.8 / ‑5.0 30K305K (AAG→AAAyfbKIPR002035/DUF3520 domain‑containing protein YfbK
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/10);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACC  >  NC_000913/2384612‑2385050
                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
cGTTACCAACATCTGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTACTAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGAACTGCCt                                                                                                                                                                                              >  2:286150/1‑251 (MQ=255)
cGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCt                                                                                                                                                                                              <  1:76657/251‑1 (MQ=255)
                         catcatTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCgg                                                                                                                                                                     >  2:155297/1‑251 (MQ=255)
                                  cTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAAttt                                                                                                                                                            >  1:230044/1‑251 (MQ=255)
                                     gTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCg                                                                                                                                                         >  2:94008/1‑251 (MQ=255)
                                             tGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATg                                                                                                                                                  >  2:64699/1‑250 (MQ=255)
                                               cTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCAcc                                                                                                                                                                                                              >  2:312271/1‑188 (MQ=255)
                                               cTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCAcc                                                                                                                                                                                                              <  1:312271/188‑1 (MQ=255)
                                               cTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACt                                                                                                                                               >  1:28003/1‑251 (MQ=255)
                                                 aTTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACt                                                                                                                                               <  2:231564/249‑1 (MQ=255)
                                                                    gACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCgg                                                                                                                                                                     >  2:106861/1‑208 (MQ=255)
                                                                    gACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCgg                                                                                                                                                                     <  1:106861/208‑1 (MQ=255)
                                                                                               tttttGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCa                                                                                                <  1:358503/250‑1 (MQ=255)
                                                                                                   acaCCATTCATTCAATGCATTTTGGATCTTAAAGGCCACCTTTAAAGCCACCGCCGTTGGCTAATAAAAGGCTTTGGATCCCTCCCTTAATAAACCCTTTCGTCGCCTATTGATAGGCTATTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAAGTCTCGag                                                                                            <  2:87514/248‑2 (MQ=255)
                                                                                                       cATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCa                                                                                        >  1:307073/1‑250 (MQ=255)
                                                                                                         gtGATTCAATCGATTTTGGATCGTAAATGACAACGTTAAAGACACCGTCAGTGGATAATAAAATGCGATTGAGGCCGCCCTTAATAAACCCTGTCGTCCCCTGCGGATAAGCCAGTTACAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAAt                                                                                      <  2:29949/249‑1 (MQ=255)
                                                                                                                    gATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCg                                                                                                             >  1:151820/1‑216 (MQ=255)
                                                                                                                    gATTTTGGATCGTCAATGCAAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCg                                                                                                             <  2:151820/216‑1 (MQ=255)
                                                                                                                     aTTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTc                                                                          >  2:183382/1‑250 (MQ=255)
                                                                                                                                              aaaGTCACAGTCAGAGGCTAATAAAATGAGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGACGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTc                                                 >  1:320492/1‑250 (MQ=255)
                                                                                                                                                    aCCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGa                                            >  2:261548/1‑249 (MQ=255)
                                                                                                                                                       gTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACt                                       >  2:137276/1‑251 (MQ=255)
                                                                                                                                                                taataaCATGCGATTGATGACGCCCTTAATCAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTg                              >  2:150279/1‑251 (MQ=255)
                                                                                                                                                                  ataaAATGCGATTGCTGCCCCCCTTATAAACCCCTTTCGTCGCCTGCTGATAAGCCATTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGc                            <  2:361238/251‑1 (MQ=255)
                                                                                                                                                                    aaattGCGATTGATGACGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGACACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCATCGAAGGTCACGATGGCAATGTTATCCTGCTCACGATGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTgtca                           >  2:297248/5‑247 (MQ=255)
                                                                                                                                                                            aTTGATGCCGCCCTTAATAACCACTTTCGTCGCCTGATGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAAAGAAGACTGGATAAGTGTCAAACGTtcat                  <  2:345698/251‑1 (MQ=255)
                                                                                                                                                                               gATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcatca                <  2:43508/250‑1 (MQ=255)
                                                                                                                                                                                  gccgccCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGaa             >  1:146066/1‑250 (MQ=255)
                                                                                                                                                                                        cTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTg      >  1:20813/1‑251 (MQ=255)
                                                                                                                                                                                              aaaCCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAAcc  >  1:201558/1‑249 (MQ=255)
                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACC  >  NC_000913/2384612‑2385050

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