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.5% 78.9 / NA 30K305K (AAG→AAAyfbKIPR002035/DUF3520 domain‑containing protein YfbK
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (16/12);  minor base G (0/1);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACC  >  NC_000913/2384564‑2385050
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
cctaTACTTTCTGCGCTTAAGAGAGTGAAACAATTTATCTTAAGTGTCCTTGACCAACATCGACAATTCGCACCATCATTTCCTCGTTGTAATTTCTATTCACCTCTCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAAGCGATTTTGGATCGTCAATGCCAACGTTACAGTCATCGTCAGTGGCTAATAAAATTCGATTGATGCCGCCCTTAATAAACCCGTtcgtc                                                                                                                                                                                                                                              <  2:400109/248‑1 (MQ=255)
     aCTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCctg                                                                                                                                                                                                                                         <  1:397760/251‑1 (MQ=255)
     aCTTTCTGCGCGTTAGAGAGGGTATCAAGGTAGTTGTAGTTGCCGTAACCACCACCGGCAATTAGCACTATCATTGCCTCGCTGTAATTGCTATTCCCCACGCCTAAAGTCGACAGAGTAAAACCAGACTATCGTTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGACAACGTTAAAGTATCCGTCAGTGGCTAATAAAATGCGATTGATGCAACCCTTAATAAACCCTGTCGTCGTctg                                                                                                                                                                                                                                         <  2:36056/251‑1 (MQ=255)
                      gagGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCa                                                                                                                                                                                                                        <  1:320454/251‑1 (MQ=255)
                      gagGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACAATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCATTTCGTAGCCTGCTGATAAGACAGTTcc                                                                                                                                                                                                                         >  1:292808/1‑250 (MQ=255)
                            aTCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCgg                                                                                                                                                                                                                  >  1:72010/1‑251 (MQ=255)
                                  gTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGc                                                                                                                                                                                                            >  1:29332/1‑251 (MQ=255)
                                  gTAGCTGTAGTTGCCGTTACCAACAGCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCGCCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCg                                                                                                                                                                                                             <  1:189569/250‑1 (MQ=255)
                                   caGCTTTAGTTGCCGTTACCACCATCTGCAATTCGCGCCATCATTGCCTCTTTGTAATTGCTATGCCCCACGCCAAACGTCGTCAGAGTAACAGCAGACGCCCGCTGTTTTTTGACCATTGATTCAAACGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGGGATTGATGCCGCCCTGAATAACCCCTTTCGACGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGcc                                                                                                                                                                                                           <  2:393722/250‑1 (MQ=255)
                                         tAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCCCCCCCATTGGt                                                                                                                                                                                                     >  2:91373/1‑251 (MQ=255)
                                            ttGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCAAGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACt                                                                                                                                                                                                  >  1:191457/1‑251 (MQ=255)
                                                          aTCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATcc                                                                                                                                                                                     >  1:366587/1‑250 (MQ=255)
                                                                    cGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCACTTGGTAATGCCTTCGGCATCCAGCGAATCAAt                                                                                                                                                                          >  2:28593/1‑251 (MQ=255)
                                                                               tgcctcttgttAATTGCTATCCCCCCCGCCAAACGTCGACAGAGTAATACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACGGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCGTAAAAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTa                                                                                                                                                                <  1:28659/240‑1 (MQ=255)
                                                                               tGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGCCAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTACCTTCGGCATCCAGCGAATCAATTGCGGCATTa                                                                                                                                                                >  2:214141/1‑250 (MQ=255)
                                                                                           aTTGCTATTCCCCACTCCAACCGACTATAAAGAAACCCAGGACCCCCGTTGTTTTTTGACAATTGATTAAATCGATTTGGGATCGTAAATGCAAACGTAAAAGTCACCTTAAGTGGCTAAAAAAATGCGATTGATGCCGCCCTTAATAACCCTTTTCGTCGCCTGCTGATAAGCCATTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTAt                                                                                                                                                   <  2:137708/251‑1 (MQ=255)
                                                                                                     cccaccccaaccgtcaaatgagTAACCCCAGTCTCCCGTTTTTTTTTGCCCATTGATTCAACCGATTTTGGATCGAAAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATACCCCCTTTCGTCGCCTGCTGATAAGCCATTCCCACCCCGGCACCGCCATTGTTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGga                                                                                                                                         <  2:174801/232‑1 (MQ=255)
                                                                                                             aaaCGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCCCTCGCTTTTGGATAGTCAATGCAAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGGTGCCGCCCTTAATAAACCCTTTCGACGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCCGCGCATCAATTCCGGCATAAATTTCCGCCTTATGACTCCCGGAGATAGa                                                                                                                                   >  1:292525/1‑249 (MQ=255)
                                                                                                                    aacatagTAACCCCGGACCCCCGTTGTTTTTTGACCATGGATTAAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCCCCGTCAGTGGCTAATAAAAGGCGATTGATCCCGCCCTTAATAACCCCTTCCGTCGCCTGCGGATAAGCCAGTTCCACCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatgc                                                                                                                          <  2:443031/246‑1 (MQ=255)
                                                                                                                         agTAACCCCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAAGGCCAACGGTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTTCGTCGCCTGCTGATAAGCCAGTGCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATAc                                                                                                                     <  1:354137/251‑1 (MQ=255)
                                                                                                                               acCAGACTCCCGCTGGTTTTGGACCTTTGATTAAATCGATTTAGGATCGTCAATGCCAACGGTAAAGTCCCCGTCAGTGGCTAATAAAATGCGATGGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGt                                                                                                               <  1:387462/251‑1 (MQ=255)
                                                                                                                                  aGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAACCCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGc                                                                                                            <  1:419737/251‑1 (MQ=255)
                                                                                                                                          cGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGt                                                                                                      >  1:316719/1‑250 (MQ=255)
                                                                                                                                                             aTTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGtt                                                                                   >  1:437637/1‑250 (MQ=255)
                                                                                                                                                                                   aaTGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCTTTCGTCGCCTGATGATAAGCCAGTTCCAGCCCGGAACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCCTTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCTCGATGGCAATGTTATCCTGCTCACGAAGTTCTTTgc                                                            >  1:421654/1‑250 (MQ=255)
                                                                                                                                                                                               aaaatcaccagcAGTGGCTAATAAAATGCGATTGATGCCGCCATTAATCAACCCTTTCGTAGCCTGCTGATAAGACAGTTCCAGCACGGCACCGCAATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCCTTCATTTCCGCCTTATGACTCCCGTAGATAGAAGGCAATGCCATACGGGAGTCGCCAGCGTGGGACCAGATGGCAATGTTATCCTGCTCAAGAATTTCTTTGACCCATAATTTc                                                 >  2:321997/12‑250 (MQ=255)
                                                                                                                                                                                                                                           taataaACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAACCGTTCATCAGAAATCATTGa     >  2:354520/1‑250 (MQ=255)
                                                                                                                                                                                                                                             ataaACCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTg                                                             >  2:143701/1‑192 (MQ=255)
                                                                                                                                                                                                                                             ataaACCCTTCCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTg                                                             <  1:143701/192‑1 (MQ=255)
                                                                                                                                                                                                                                                aaCCCTTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAAcc  >  2:169148/1‑248 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
TCAATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTGACCATTGATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGATTGATGCCGCCCTTAATAAACCCCTTCGTCGCCTGCTGATAAGCCAGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAATCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACC  >  NC_000913/2384564‑2385050

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