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→AAAyfbKIPR002035/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

ATTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATC  >  NC_000913/2384562‑2385042
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
aTTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATACACCCTTtcg                                                                                                                                                                                                                                        >  2:118853/1‑251 (MQ=255)
      tACTTTCCGCGATTCCGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCACCATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACCGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCATTTCGTCGCct                                                                                                                                                                                                                                  >  2:139559/1‑251 (MQ=255)
                     agagagTGTATCAATGTAGCGGTAGTTGCCGTGACCAACATCGGCAATCCGCACCATCATTGCCTCGTTGTAAGTGCTATTCCCCACGCCAATCGTCGACAGAGTAACACCAGACTCCCGCTGTTTGTTGACCATTGATTCAATCGATTTTGGGTCGTCAATGCCAACTTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCTCTTAATAAACCCTTTCGTCGCCTGCTGAGAATCCAGtt                                                                                                                                                                                                                   <  2:192788/251‑1 (MQ=255)
                                                     tACCAACATCGGCAATTCGCCCCAACATTGCCTCGTTTTAATGGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGGCTCCCGCTGTTTTTTGACCATTGATTCAACCAATTTTGGATCGCCAATGCCAACGTTAAAGCCACCGTCAGTGGCTAATAAAATGCGATGGATGCGGCCCTGAATAAACCCTTTCGTCGCCGTCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTc                                                                                                                                                                                    <  1:384157/250‑1 (MQ=255)
                                                                        cACCATCATTGCCTCGTTGTAATTGCTATCCCCCACGCCAAACGGCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGTTCGTCAATGCCAACGTTAAGGTCACCGTCAGTGGCTAATAAAAGGCGATGGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTTATAAGCCAGTGCCACCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTg                                                                                                                                                                <  1:44114/251‑1 (MQ=255)
                                                                                            aaTTGTTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGTTGTTTTTTGACCATTGTTTCAATCGATTTTGGATCTTCAATGCCAACTTTAAAGTCCCCGTCAGTGGCTAATAAAAGGCGATGGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCtt                                                                                                                                             <  1:227179/250‑1 (MQ=255)
                                                                                               tgctaatcaacACGCCAAACGTAGACAGAGAAACACCAGACTCCCGCTGTCTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGAGGCATTAATTTCCGCCTTATg                                                                                                                                          >  2:383123/11‑250 (MQ=255)
                                                                                                 cTATTCCCCACGCCAAACGTCGCCAGAGTAACCCCAGACTCACGCTGTTTTTTGACCATTGATTCAATCGATTTTGGCTCATCAATGCCAACGTTAAAGTCAACGTCAGTGGATAATAGAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCGGCTGATCAGCCAGTTCCTGCCCGGTACCGCCATTTGTACTGCCTTCGGCATCCAGCGAATCAATTGCCGCATTAATTTACGCCTTATGACt                                                                                                                                       >  2:79118/1‑251 (MQ=255)
                                                                                                         cACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATAGATTCAATCGATTCTGGCTCGACAATGCCAACGTTAACGTCACCGTCAGTGGCTAATAACATGCGATTGATGCCGCCCTTACTAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGACTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCgg                                                                                                                                  >  1:62581/1‑248 (MQ=255)
                                                                                                            gCCAAACGTCGACAGCGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGCTTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGCTAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGAt                                                                                                                              >  2:329216/1‑249 (MQ=255)
                                                                                                                                                  ttttGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCTTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCAc                                                                                       <  1:63927/250‑1 (MQ=255)
                                                                                                                                                    ttGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGAt                                                                                    >  1:179421/1‑251 (MQ=255)
                                                                                                                                                       aCCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGc                                                                                 >  1:167463/1‑251 (MQ=255)
                                                                                                                                                          aTTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGAACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATATGGGAGTCGCCAGCGTAGGTCACGATAGCAAt                                                                              >  1:222655/1‑251 (MQ=255)
                                                                                                                                                            ggatccAATCTATTTTGGTTCGGCAATGCCACCTTTAAAGTCCCCGTCAGTGTCTAATAAAATGCGATTGTTGCCGCTTTTATTAAACCCTTCCGTCTCCTGCTTATAAGCCGGTCCCTGCCCGGCACCGCCATTGGTACGGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCTCCTTATGACTCCCTGAGATAGAAGGTTATGCAATACCGGAGTCGCCAGCGTAGGTCACGATGGCAATGt                                                                            <  1:178008/246‑1 (MQ=255)
                                                                                                                                                                    cctaTTTTGGATAGTCAATGCCAACGTTAAAGTACCCGCCAGTTGCTAATAAAATGCAATGGATCCCGCCCTTACTACCCCCTGTCGTCGCCTGCTGATAACCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGTCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCTCCAGCGTAGGTAACGATGGCAATGTTATCCttc                                                                    <  2:93271/248‑3 (MQ=255)
                                                                                                                                                                             gATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGt                                                           <  1:118853/251‑1 (MQ=255)
                                                                                                                                                                                                aaaGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTc                                         >  1:197556/1‑250 (MQ=255)
                                                                                                                                                                                                 aaGTCACCTTCAGTGGCTAATAAAATGCGATTGATGCCGCCTTGAATAACCGCTTGCGCCGCCGGCGTATAAGCCAGTTCCAGCCCGGCGCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATGGCGGCATTAAGTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGAGGGGAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCaa                                       <  1:165532/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gccCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATc  >  1:163179/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            
ATTCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATC  >  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: