Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 246,747 G→T S12S (TCG→TCT)  yafL → NlpC/P60 family protein YafL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913246,7470GT95.8% 60.7 / ‑3.2 24S12S (TCG→TCT) yafLNlpC/P60 family protein YafL
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (15/8);  total (16/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

ATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAA  >  NC_000913/246509‑246958
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                   
aTCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGATTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCCGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGtt                                                                                                                                                                                                          <  1:211427/250‑1 (MQ=255)
aTCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGtt                                                                                                                                                                                                          <  1:272005/250‑1 (MQ=255)
       gtgCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGAGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCAATGATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTgttatgtt                                                                                                                                                                                                  >  2:77124/1‑244 (MQ=255)
         gCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCGGTCGATTCCATCCTGTGT‑ATTGTCTGGATTACTGTTGAtttgttt                                                                                                                                                                                                 <  1:231282/250‑1 (MQ=255)
                              aaTTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTACTTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACAGCAACATTAAAGCATGAACAGCCATTATTGGGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAgc                                                                                                                                                                           >  1:274576/1‑251 (MQ=255)
                                       tCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCcaccac                                                                                                                                                                   >  1:112920/1‑250 (MQ=255)
                                                                                                           tgtatatGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGAAACATTAAAGCATGACCAGCCATTATTGAGTATGACATTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTTTTGATTTGTTTGACGTTTTCTTCTTTTGCCTGCGCCACCACATCATATATCTCTTTCTGATACTCCGCCCGCCAGCGGACGCAAAACAGCGCGCTATTATTAAAACAGt                                                                                               >  2:83905/3‑250 (MQ=255)
                                                                                                                 tGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAAc                                                                                        >  2:24632/1‑251 (MQ=255)
                                                                                                                   ttAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACt                                                                                       <  1:151820/250‑1 (MQ=255)
                                                                                                                                      ttattaAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGcc                                                                     >  1:260694/1‑249 (MQ=255)
                                                                                                                                               ttggtCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACTCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCACACTCATCCGAAAAAGCAGGCCAGCTAAATtgt                                                          >  2:570/5‑251 (MQ=255)
                                                                                                                                                tGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGAATGACCAGCCATTATTGAGCATGTCCTTGCCGTCGATTCCATCCTTTCT‑ATTGTCTGGATTACTGTTGATTTGTTAGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGACAGCCAGCGGATGCAAAACCGTGCGCGTTTATTCACCCAGAACCAAACTCATCTGACAATGCAGACCAGCTATATtgtg                                                         >  1:141219/1‑251 (MQ=255)
                                                                                                                                                           gATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGCGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTTATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTACAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGc                                              >  2:156842/1‑251 (MQ=255)
                                                                                                                                                              ccGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTCCCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGc                                              >  1:150909/1‑248 (MQ=255)
                                                                                                                                                               cGGGATACACTGCAACTTTAAAGCATGCCCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑TTTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGaa                                          <  1:228053/251‑1 (MQ=255)
                                                                                                                                                                  gatacaaggCAAAATTAATGCATGACCTGCCATTATAGAGTCTGTCCTTGCCGACGATTCCATCCTTTGT‑ATTGTCTGGATTCCTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACTACATCACATATCTCTTTCAGCTTCGCCGCCCGACAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGACCCAAACTCATCTGAAAAAGCAGGCCACCAATATTGTGAAAGGAAATGCCAAAAGc                                       >  2:252191/9‑251 (MQ=255)
                                                                                                                                                                               cctatctaatccccAGCCATTATTTAGTTTTTCCTTGCCGTCGATTCCATCCTT‑GTGATTGTCTGGATTACTTTTGATTTGTTTGCCGTTTTCTTCATTGCCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATAGTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTa                           <  2:191187/238‑1 (MQ=255)
                                                                                                                                                                                     gCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCCTTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACCCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAAGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGc                    >  2:285611/1‑251 (MQ=255)
                                                                                                                                                                                      cATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTa                           >  1:191187/1‑243 (MQ=255)
                                                                                                                                                                                         acccatctattatttagtgttcccTTGCCGTCTATTCAACCCTTTGT‑ATTGTCTGTCTTACTGTTGATTGTTTTGCTGTTTTTTTAATTTGCCAGCGCCCCCACATCACATATCTCTTTCAGCTACGCAGCTCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACaa                <  2:9268/229‑1 (MQ=255)
                                                                                                                                                                                           ccagccagTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAGGAAAAAGGGCGCTCCGCCAGCACAAcc              >  2:105601/1‑251 (MQ=255)
                                                                                                                                                                                                cATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTTT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAg         <  1:254006/251‑1 (MQ=255)
                                                                                                                                                                                                    attGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGc        >  1:264938/1‑248 (MQ=255)
                                                                                                                                                                                                        agaaTGTCCTTGCCGTCGATTCCATCCTTTGA‑ATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCAAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGAGCGCGTTTATTAAAAAAGTACCAAACTCATCTGAAAAAGCAGGACAGCTATATTGTGGAAGGCAATGACGAAAGCAAAAGGGCGCTACGCCCGCACAACCGGGCGCAGATaa  >  1:238061/4‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                   
ATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGT‑ATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAA  >  NC_000913/246509‑246958

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