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,8080CT90.0% 50.3 / ‑3.7 20K305K (AAG→AAA) yfbKIPR002035/DUF3520 domain‑containing protein YfbK
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (11/7);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATC  >  NC_000913/2384562‑2385042
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
aTTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATACACCCTT‑t‑cg                                                                                                                                                                                                                                        >  2:118853/1‑251 (MQ=255)
      tACTTTCCGCGATTCCGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCACCATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACCGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCATT‑T‑CGTCGCct                                                                                                                                                                                                                                  >  2:139559/1‑251 (MQ=255)
                     agagagTGTATCAATGTAGCGGTAGTTGCCGTGACCAACATCGGCAATCCGCACCATCATTGCCTCGTTGTAAGTGCTATTCCCCACGCCAATCGTCGACAGAGTAACACCAGACTCCCGCTGTTTGTTGACCATTGATTCAATCGATTTTGGGTCGTCAATGCCAACTTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCTCTTAATAAACCCTT‑T‑CGTCGCCTGCTGAGAATCCAGtt                                                                                                                                                                                                                   <  2:192788/251‑1 (MQ=255)
                                                     tACCAACATCGGCAATTCGCCCCAACATTGCCTCGTTTTAATGGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGGCTCCCGCTGTTTTTTGACCATTGATTCAACCAATTTTGGATCGCCAATGCCAACGTTAAAGCCACCGTCAGTGGCTAATAAAATGCGATGGATGCGGCCCTGAATAAACCCTT‑T‑CGTCGCCGTCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTc                                                                                                                                                                                    <  1:384157/250‑1 (MQ=255)
                                                                        cACCATCATTGCCTCGTTGTAATTGCTATCCCCCACGCCAAACGGCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGTTCGTCAATGCCAACGTTAAGGTCACCGTCAGTGGCTAATAAAAGGCGATGGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTTATAAGCCAGTGCCACCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTg                                                                                                                                                                <  1:44114/251‑1 (MQ=255)
                                                                                            aaTTGTTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGTTGTTTTTTGACCATTGTTTCAATCGATTTTGGATCTTCAATGCCAACTTTAAAGTCCCCGTCAGTGGCTAATAAAAGGCGATGGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCtt                                                                                                                                             <  1:227179/250‑1 (MQ=255)
                                                                                               tgctaatcaacACGCCAAACGTAGACAGAGAAACACCAGACTCCCGCTGTCTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGAGGCATTAATTTCCGCCTTATg                                                                                                                                          >  2:383123/11‑250 (MQ=255)
                                                                                                 cTATTCCCCACGCCAAACGTCGCCAGAGTAACCCCAGACTCACGCTGTTTTTTGACCATTGATTCAATCGATTTTGGCTCATCAATGCCAACGTTAAAGTCAACGTCAGTGGATAATAGAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCGGCTGATCAGCCAGTTCCTGCCCGGTACCGCCATTTGTACTGCCTTCGGCATCCAGCGAATCAATTGCCGCATTAATTTACGCCTTATGACt                                                                                                                                       >  2:79118/1‑251 (MQ=255)
                                                                                                         cACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATAGATTCAATCGATTCTGGCTCGACAATGCCAACGTTAACGTCACCGTCAGTGGCTAATAACATGCGATTGATGCCGCCCTTACTAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGACTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCgg                                                                                                                                  >  1:62581/1‑248 (MQ=255)
                                                                                                            gCCAAACGTCGACAGCGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGCTTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGCTAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGAt                                                                                                                              >  2:329216/1‑249 (MQ=255)
                                                                                                                                                  ttttGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCTTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCAc                                                                                       <  1:63927/250‑1 (MQ=255)
                                                                                                                                                    ttGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGAt                                                                                    >  1:179421/1‑251 (MQ=255)
                                                                                                                                                       aCCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGc                                                                                 >  1:167463/1‑251 (MQ=255)
                                                                                                                                                          aTTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGAACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATATGGGAGTCGCCAGCGTAGGTCACGATAGCAAt                                                                              >  1:222655/1‑251 (MQ=255)
                                                                                                                                                            ggatccAATCTATTTTGGTTCGGCAATGCCACCTTTAAAGTCCCCGTCAGTGTCTAATAAAATGCGATTGTTGCCGCTTTTATTAAACCCTT‑C‑CGTCTCCTGCTTATAAGCCGGTCCCTGCCCGGCACCGCCATTGGTACGGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCTCCTTATGACTCCCTGAGATAGAAGGTTATGCAATACCGGAGTCGCCAGCGTAGGTCACGATGGCAATGt                                                                            <  1:178008/246‑1 (MQ=255)
                                                                                                                                                                    cctaTTTTGGATAGTCAATGCCAACGTTAAAGTACCCGCCAGTTGCTAATAAAATGCAATGGATCCCGCCCTTACTAC‑CCCCTGT‑CGTCGCCTGCTGATAACCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGTCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCTCCAGCGTAGGTAACGATGGCAATGTTATCCttc                                                                    <  2:93271/248‑3 (MQ=255)
                                                                                                                                                                             gATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGt                                                           <  1:118853/251‑1 (MQ=255)
                                                                                                                                                                                                aaaGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTc                                         >  1:197556/1‑250 (MQ=255)
                                                                                                                                                                                                 aaGTCACCTTCAGTGGCTAATAAAATGCGATTGATGCCGCCTTGAATAA‑CCGCT‑TGCGCCGCCGGCGTATAAGCCAGTTCCAGCCCGGCGCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATGGCGGCATTAAGTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGAGGGGAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCaa                                       <  1:165532/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gccCTTAATAAACCCTT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATc  >  1:163179/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
ATTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCT‑T‑CGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATC  >  NC_000913/2384562‑2385042

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