Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,630,221 C→T E128K (GAA→AAA)  yhiI ← putative membrane fusion protein YhiI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,630,2210CT96.7% 83.6 / ‑4.1 30E128K (GAA→AAA) yhiIputative membrane fusion protein YhiI
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (13/16);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

ATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCG  >  NC_000913/3629979‑3630438
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                         
atCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTgg                                                                                                                                                                                                                    <  1:122138/250‑1 (MQ=255)
                cGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGa                                                                                                                                                                                                   <  1:310292/251‑1 (MQ=255)
                       gggCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCg                                                                                                                                                                                            <  1:182333/250‑1 (MQ=255)
                            tgtgCCGCTTCGACGCGGTTGTGCGCCGTCATGATTTTGGTGCGTTCCGCTTCTCTCGCCGCTTTAGAAGCTGATACCTGATCTTTCGCCGATTCCAGCGCAGCTCGGTCGCTCTCAGCGTCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCTGTGATCGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCat                                                                                                                                                                                        <  1:356680/250‑2 (MQ=255)
                             gtgCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATCGCCGCTTTAGAAGCCGATCCCTGTGCTTTCGCCGATTCCCGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCCGCTGTTGCGCAGAAATAGCCCCTCGTTGGACCAGTGAACGGGAACGCGTATGAAGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACg                                                                                                                                                                                       >  1:352275/1‑250 (MQ=255)
                             gtgCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCAGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACg                                                                                                                                                                                       >  1:239199/1‑250 (MQ=255)
                                  gCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATAACTGCGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTc                                                                                                                                                                                 >  1:250959/1‑251 (MQ=255)
                                                  gcgcCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTcc                                                                                                                                                                 <  1:110476/251‑1 (MQ=255)
                                                                                 tataGCCGCTTTAGATGCCGATACCTGAGCTTTCGCCGATTCCAGCGCCGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCCGAAATAGCCCCTCGTTGGGCCAGTGTACGGGAACGCGTATGACGTTTTGCTACGGATTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAATTTCGCTTTGTCGTTCACCCAGCAAAGCCTGCCCGGAAGCAACGGCGCttt                                                                                                                                  >  2:144906/1‑251 (MQ=255)
                                                                                  attgCCGCTGTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTtg                                                                                                                                 <  1:175684/248‑1 (MQ=255)
                                                                                               gaggcgataccggggatttCGCCGATTCCAGCGCAGCTCGGTCGGTCTCAGCTGCGTCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCCGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGTTCGAGTTTCGCGTTGTCGTGGCGCCAGCAAAGCCTGCGCTGCAGCAACGGCGCTGTGTGCCTCTGTGAtt                                                                                                                    <  2:351110/235‑1 (MQ=255)
                                                                                               aaGCCTATATCTGAGTTGTCGCCTATTCAGGCGAAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCCCTTGGGGCCAGTGAACGGGAATGCGTATGACTCTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCATCAACGGCGCTTTGTGCCTCTTTGAtt                                                                                                                    <  2:193294/251‑1 (MQ=255)
                                                                                                                     cGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCt                                                                                              <  1:48811/251‑1 (MQ=255)
                                                                                                                     cGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTACAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGATTCGATTTGTCATTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCt                                                                                              >  2:334286/1‑251 (MQ=255)
                                                                                                                                         cgcTCTTAGCGGCGTCGCGATCGTCAGCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAAGGCGTTTGACGTTTTGCTACGGAGTCCAGTTTGGCCTGGCGTTGATTACCCAGCGACTGGGCTGCACGAGGTTCGGTTTGTCGTTGCTCCAGCAAAGCCTGCGCGTCAGCAACGGCTCTTTGTTCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTa                                                                          <  1:183033/251‑1 (MQ=255)
                                                                                                                                              tcAGCGGCGGCGCGCTCGTCATCCAGCTGTTGCGCAGAAATAGCCCATCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAATCGCTGTTCCTGCAACACGCGAGTATCCAt                                                                     >  1:316289/1‑251 (MQ=255)
                                                                                                                                               cAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCAt                                                                     <  1:9139/250‑1 (MQ=255)
                                                                                                                                                 gcggcggcgCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGAATCCATCtt                                                                  >  1:294004/1‑251 (MQ=255)
                                                                                                                                                   ggcggcgCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCg                                                                <  1:68074/251‑1 (MQ=255)
                                                                                                                                                     cggcgCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGcc                                                              >  1:72934/1‑251 (MQ=255)
                                                                                                                                                      ggcgCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGcc                                                              >  1:298789/1‑250 (MQ=255)
                                                                                                                                                                                cAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTg                                   <  1:42012/251‑1 (MQ=255)
                                                                                                                                                                                 aGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAACGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGCTGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCAATTCACCTTCGCGCACAAACTTGc                                  >  2:200067/1‑251 (MQ=255)
                                                                                                                                                                                  gAAATAGCCCCTTGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGcc                                 <  1:155140/251‑1 (MQ=255)
                                                                                                                                                                                     aTAGCCCCTCGATGGGCCAGTACACGGGAACGAGTATGACGATTTGCTAAGGAGTCCAATTTTGCCGGACGTTGAGTAACAAGCGCCTGAGCGGCACGCGTTTCGCTTTGTCGTAGCTCCAACAAAGCCTGCGCGACCGCAACGGCGCTATGTGCATCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCGGCAACACGCGAGTATACATATTAGCCCGCACTTCACCTTCGCGAAACAACTTGCcttc                              >  2:233054/1‑251 (MQ=255)
                                                                                                                                                                                      gagcaccTCGTTGGGCCAGTGAACGGGAAAGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTGCGCCAGCACTTCACCTTCGCGAACAAACTTGCcttct                             <  1:281434/246‑1 (MQ=255)
                                                                                                                                                                                                   ggCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATc                                                                        >  2:150336/1‑195 (MQ=255)
                                                                                                                                                                                                   ggCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATc                                                                        <  1:150336/195‑1 (MQ=255)
                                                                                                                                                                                                            aCGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAACGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACACACTTGCCTTCTTTCACCAGAATGGTGTCGATAc       >  1:53228/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaCGCGTATGACGTTTTGCTACGGAGACCAGTTTTACCTGGAGTTTATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCTTCTTTGATTTGCGCAATGGCTTCCAGTCGCTGTTACTGCAACACGAGAGCACCCATCTTCGCCACCCCCTCACCTTCGCGACCAAACTTGCCTTCTTTCACCACAATGGTGTCGATACGCCCg  >  2:269486/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                         
ATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCG  >  NC_000913/3629979‑3630438

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