Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 234,876 T→C 10.7% W21R (TGG→CGG)  yafS → putative S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913234,8760TC10.7% 81.3 / 4.0 38W21R (TGG→CGG) yafSputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (10/23);  new base C (2/2);  total (12/26)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

TCAAAACCCAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTG  >  NC_000913/234731‑235023
                                                                                                                                                 |                                                                                                                                                   
tCAAAACCCAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTggg                                                                                                                                                  <  1:745696/149‑1 (MQ=255)
tCAAAACCCAGATGGAAATGACAGCAAAGGCGGTAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTggg                                                                                                                                                  <  1:213116/149‑1 (MQ=255)
    aaCCCAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGaaa                                                                                                                                              <  2:236730/149‑1 (MQ=255)
       ccAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCt                                                                                                                                           >  2:774864/1‑149 (MQ=255)
                              cGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGGCCCTCCGACTGACGTGGATAACGCTTGCAGAAGCGATTTCCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGGGCGcc                                                                                                                    >  1:149835/1‑149 (MQ=255)
                              cGGGAATACTGTTAAGATTCATAGATTAAACACGAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTCTATCGCAAGGCGCGGGATCGcc                                                                                                                    >  2:363872/1‑149 (MQ=255)
                                 gAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGc                                                                                                                 <  1:774864/149‑1 (MQ=255)
                                     aCTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCaa                                                                                                             <  1:110288/149‑1 (MQ=255)
                                        gTTAAGTTTCATAGATTACCTCTCGGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCGAACTGTCGAGGCTCCTGATTGGTGGGGCGATTTGCCCTGGGGAAAGCTTTATGGCAAGGCGCTGGAGCGCCAGCTCAAccc                                                                                                          <  2:226235/149‑1 (MQ=255)
                                           gaGATTCAGATGTTACCTCTCAGGGGGGAACGGAAGGTTGTGATGAAACCGGCAAGAGGCCCTCAAACGGTCGTGGCTCCTGATTGCTGGGGCGTTGGTCGCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTg                                                                                                       <  1:347392/148‑1 (MQ=255)
                                           aaGATTTAGGGATTACCTCTCACTGTGACACGGAAGGTTGTGACGAAACCGGCAAGCGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCTAGGCGCTGGAGCGCCAGCTCAACCCGTg                                                                                                       <  1:67506/149‑1 (MQ=255)
                                            aGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTgg                                                                                                      <  2:77086/149‑1 (MQ=255)
                                                tCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                  <  2:725760/149‑1 (MQ=255)
                                                tCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCCGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                  >  2:403489/1‑149 (MQ=255)
                                                tCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCCGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCa                                                                                                  >  2:305305/1‑149 (MQ=255)
                                                     gATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTaaa                                                                                             <  2:97739/149‑1 (MQ=255)
                                                             ctcAGTGTGAAACGGCAGGTTGTGATGAAACCGGCAAGAGTCCCTAAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGGGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATgg                                                                                     >  2:191156/1‑149 (MQ=255)
                                                             ctcAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATgg                                                                                     <  2:79677/149‑1 (MQ=255)
                                                             cgcAGTGTGAAACGGAAGGTTGTGATGAAACGGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTGGCCCGGGGGAATGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATgg                                                                                     <  2:710985/147‑1 (MQ=255)
                                                                 gtgtGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGtttt                                                                                 <  2:663948/149‑1 (MQ=255)
                                                                    tgAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGCGCGCCAGCTCAACCCGTGGTTCCCTACAATGTATGGCTTTCAt                                                                              >  1:547510/1‑149 (MQ=255)
                                                                     gAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGCGTTTCATc                                                                             >  2:608836/1‑149 (MQ=255)
                                                                         cGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCt                                                                         <  1:633801/149‑1 (MQ=255)
                                                                         cGGAACGTTGTCATGAAAGCCGCAAGCGCCCCTCGAACTGTCGTGGCTCCTGATTGCTGGGGCGATTCGCCCTGGGGAAAGCTTTATCTCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCt                                                                         <  1:86485/149‑1 (MQ=255)
                                                                            aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTTTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTaa                                                                      <  2:252558/149‑1 (MQ=255)
                                                                            aaGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAGGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTaa                                                                      <  1:307188/149‑1 (MQ=255)
                                                                              ggatgtgATGAAACCGGCAAGAGTCCCTCAACCTGTCGTCGCTCCTGATTGCAGGGGCGATTTGCCCTGGGGAAAGCTGTATCGCAAGGCGCTGGTGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGa                                                                    <  2:552210/146‑1 (MQ=255)
                                                                                  gtgATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTgg                                                                <  1:97259/149‑1 (MQ=255)
                                                                                         aaCCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCGGGAGCGCCAGCTCAACCCGTGGTTCACTAACATGGATGGTTCTCATCTGCTTAAGCTTGGCACTTTa                                                         >  2:11927/1‑149 (MQ=255)
                                                                                            cGGCGAGAGTCCCTCAAGCTGTCGTGTCTCCTGATTGCTGGGGCGGTTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAgc                                                      <  2:254230/149‑1 (MQ=255)
                                                                                                         tCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAACTGTATGGTTTTCATCTGCTTAAGATTGGCGATTTAAGCGCAGAAATCAAtt                                         >  2:586021/1‑149 (MQ=255)
                                                                                                         tCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAAtt                                         >  2:698806/1‑149 (MQ=255)
                                                                                                                       cTCCTGATTGCTGGGGCGATTTGCCCCGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGcgcg                           <  1:403489/149‑1 (MQ=255)
                                                                                                                       cTCCTGATTGCTGGGGCGATTTGCCCCGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGcgcg                           <  1:305305/149‑1 (MQ=255)
                                                                                                                                    gggCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAg              >  1:282029/1‑149 (MQ=255)
                                                                                                                                             gCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGttt     <  2:190785/149‑1 (MQ=255)
                                                                                                                                                cTGGGGAAAGCTTTATCGGAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTGAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTg  <  1:618497/149‑1 (MQ=255)
                                                                                                                                                cTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTg  <  1:431909/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCAAAACCCAGATGTAATTGTCATCAAAGGCGGGAATACTGTTAAGATTCATAGATTACCTCTCAGTGTGAAACGGAAGGTTGTGATGAAACCGGCAAGAGTCCCTCAAACTGTCGTGGCTCCTGATTGCTGGGGCGATTTGCCCTGGGGAAAGCTTTATCGCAAGGCGCTGGAGCGCCAGCTCAACCCGTGGTTCACTAAAATGTATGGTTTTCATCTGCTTAAGATTGGCAATTTAAGCGCAGAAATCAATTGCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTG  >  NC_000913/234731‑235023

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