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,8080CT92.6% 74.5 / ‑4.2 27K305K (AAG→AAA) yfbKIPR002035/DUF3520 domain‑containing protein YfbK
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (12/13);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTG  >  NC_000913/2384582‑2385046
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
cAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCTGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGt                                                                                                                                                                                                                        >  2:14139/1‑251 (MQ=255)
       ggTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGcc                                                                                                                                                                                                                 <  2:343811/251‑1 (MQ=255)
         tATCAATGTAGCTGTAGTTGCCGCTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAAACCTTT‑C‑GTCGCCGGCTGATAAGCACGTTCCAGCCCg                                                                                                                                                                                                               >  2:11024/1‑251 (MQ=255)
              aTGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GGCGCCTGCTGATAAGCCAGTTCCAGCCCGGCAcc                                                                                                                                                                                                          <  1:359158/251‑1 (MQ=255)
                  aGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGcc                                                                                                                                                                                                       <  2:151208/250‑1 (MQ=255)
                                 tACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTc                                                                                                                                                                                        <  1:200666/250‑1 (MQ=255)
                                 tACCAACATCGGCAATCCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACAATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTAAAAGCCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTT‑C‑GTCGCCTGTTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTc                                                                                                                                                                                        <  2:347849/250‑1 (MQ=255)
                                           ggCAATTCGCCCCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGCCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGTTTTTGGATCGCCAATGCCAACGTTAAAGTCACCGTCAGTGTCTAATAAAAGGCGATTGATGCCGCCCTTAATAAACCCTTG‑C‑GTCGCCTGCTGATAAGCAAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCg                                                                                                                                                                             <  2:34400/251‑1 (MQ=255)
                                                                  ggTTGTAATTGCTATGCCCCACGCCAAACGTCGACAGCGTAACACAAGCCTCCCGCTGTTTTTTTACCATTGATTCAATCGATTGTGATTCGTCAAGGCCAACGTTAAAGTCACCGTCAGTGGCGAATAAAATGGAATGGAGGCCGCCCTTAATAATCCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTcc                                                                                                                                                      <  1:385763/250‑1 (MQ=255)
                                                                     tGTAATTGATATTCCACACGCCAAACGTCGACAGAGTAACACCAGACTCGCGCTGTTTTTTGACCATTGCTTCAATCGATTTAGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGGGGCTAATAAAATACGATTGATGCCGCCCTTAATAAAACCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCAACGCCATTGGTACTGCCTTCGGCATCCAGCGAATAACTTGCGGCATTAATTTccacc                                                                                                                                                   >  2:314982/1‑248 (MQ=255)
                                                                            gCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGAc                                                                                                                                            >  1:229177/1‑251 (MQ=255)
                                                                              tATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTc                                                                                                                                          >  1:108890/1‑251 (MQ=255)
                                                                                 tCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTccc                                                                                                                                        <  1:344989/250‑1 (MQ=255)
                                                                                                      gagTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatgcaat                                                                                                                   >  2:136182/1‑250 (MQ=255)
                                                                                                      gagTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTAGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatgcaat                                                                                                                   >  1:309186/1‑250 (MQ=255)
                                                                                                          aaCACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACggg                                                                                                              <  1:120883/251‑1 (MQ=255)
                                                                                                                     cccGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCg                                                                                                   >  1:81269/1‑251 (MQ=255)
                                                                                                                        ccTGTTTTTTGACCATGAATTCAATCGATTTGGGATCGTAAATGCCAACGTAAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCTCCCTAAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTAGGACTCCCGGAGATAGAAGGCGATGCAATACGGTAGGCGCCGGCGTAg                                                                                                <  2:46370/250‑1 (MQ=255)
                                                                                                                                ttGACCATTGATTCAATCGATTTTGGATCGTAAATGCCAACGTTAAAGTCACCTTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑TTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGa                                                                                         <  2:280502/250‑1 (MQ=255)
                                                                                                                                                          aTCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGCATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGtt                                                              >  1:137228/1‑251 (MQ=255)
                                                                                                                                                                    ccAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCa                                                     >  2:395257/1‑250 (MQ=255)
                                                                                                                                                                                         gtgtcTAATAAAAGCCGATTGATGCCGCCCTTAA‑AAACCCCTTTC‑GCCGCCTGCTGATAAGCCAGTTCCACCCCGGCTCCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGAt                               <  2:199346/247‑1 (MQ=255)
                                                                                                                                                                                                      tGCGATTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGt                  >  2:97517/1‑251 (MQ=255)
                                                                                                                                                                                                          aTTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtcat              >  1:275929/1‑251 (MQ=255)
                                                                                                                                                                                                                   ccccTTAATAAACCCTTT‑C‑GGCGCCTCCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCa     <  2:108890/250‑1 (MQ=255)
                                                                                                                                                                                                                     ccTTAATAAACCCTTG‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAAGTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCAtt   <  1:368297/251‑1 (MQ=255)
                                                                                                                                                                                                                         taataaCCCCTT‑CCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTg  <  1:395257/245‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
CAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTT‑C‑GTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTG  >  NC_000913/2384582‑2385046

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