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,7470GT97.8% 128.2 / NA 46S12S (TCG→TCT) yafLNlpC/P60 family protein YafL
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (22/23);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGC  >  NC_000913/246534‑246929
                                                                                                                                                                                                                     |                                                                                                                                                                                      
gCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTg                                                                                                                                                   <  1:355845/251‑1 (MQ=255)
                                             aTTTAATAAACGGATGTTTTTACATAATTAACTTAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCAGGACCAGCCATGATTGAGTATGTCCTGGCCTTCGATTCCATCCTTTGTATTGTCTGGATTACTTTTGATTTGTTGTCCTTTTCCTTCATTTGCCAGCGCCACCACATCCCATATCTCTTTCAGCTACGccgcccgac                                                                                                      <  1:474733/251‑3 (MQ=255)
                                             aTTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGTGTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCAATTACATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTGTCTTCATTTGCCAGCGCCACCACATCACATAGCTCTTTCAGCTACGCcgcccgcc                                                                                                      <  1:158457/251‑1 (MQ=255)
                                                        ggATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAACGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCa                                                                                            >  1:248144/1‑250 (MQ=255)
                                                                  aCATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTAcgccgc                                                                                                           <  1:353755/225‑1 (MQ=255)
                                                                    aTAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGcgcgn                                                                               <  1:404301/251‑2 (MQ=255)
                                                                                                     ataatGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGTGTATGTCCTTGCCGTCGCTTCCTACCTTTGTATTATCTGGATTAATGTTGATTTGTTTGCCGTTTTCTTAATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCAGATGCAAAACCGTGCGCGTTTTTTAAAACAGTACCAAACTCATCTGaaaaa                                               >  1:242033/1‑250 (MQ=255)
                                                                                                                           naGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTaaaa                                                                      >  1:314262/2‑205 (MQ=255)
                                                                                                                           naGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTAATTGCCGTCGATTCCATCCTTTGCATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCAACCACATCACATATCTCTTTCAGCTAcgccgc                                                                                                           >  1:44635/2‑168 (MQ=255)
                                                                                                                            aGAAACGAACCGGGATACGCTGCAACATTAAATCATGACCAGCCATTATTGAGTATGTCCTTTCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTg                                                                                                                                                                            <  1:506029/102‑1 (MQ=255)
                                                                                                                                ccnntctgngacacactctcaCATTAATGAATGNCCAGCCATTTTGGATTATTTCCTGGCCGTCTNTNCCATCCTTTGTATTGTCTGGATTCCTGTGGATGTGTTGGCCGTTTTCTNCATGTGCAAGNGCCNNCACANCNCNTATCTCTTTCATCTACGCCGCCCGCCAGCNGATGCAAAACCGTGCGCGTTTATNNAANNNGTNNCAAAc                                                           <  2:364028/191‑1 (MQ=255)
                                                                                                                                   aTCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTGTTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATCTCTCTTTCCGCTAcgccg                                                                                                            >  1:296776/1‑159 (MQ=255)
                                                                                                                                    tCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGc                                                                                                                   <  1:141172/151‑1 (MQ=255)
                                                                                                                                     ccGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCan                                                                                                                                                      <  1:310328/115‑2 (MQ=255)
                                                                                                                                     ccGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCAt                                                                                                                                                      >  2:310328/1‑115 (MQ=255)
                                                                                                                                      cGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGtt                                                                                                                                                                             >  2:318170/1‑91 (MQ=255)
                                                                                                                                      cGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGtn                                                                                                                                                                             <  1:318170/91‑2 (MQ=255)
                                                                                                                                       gggATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATACTTTGTATTGTCTGGCTTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAAc                                                           >  1:364028/1‑204 (MQ=255)
                                                                                                                                             cacTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGACAGAGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCaaa       >  1:202360/1‑250 (MQ=255)
                                                                                                                                             cacTGCAAAATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAg                        <  1:27264/233‑1 (MQ=255)
                                                                                                                                               cTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGAGTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTAcgccga                                                                                                           <  1:33799/148‑2 (MQ=255)
                                                                                                                                                 gCAACATTAAAGCATGACCAGCCATTATTGCGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTAcgccgc                                                                                                           >  1:424510/1‑146 (MQ=255)
                                                                                                                                                 gCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTgcgc                                                                                 >  1:10023/1‑172 (MQ=255)
                                                                                                                                                              aTGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATGTGTTTGCCGTTGTCTTCATTTGCCAGCGCCACcacatcac                                                                                                                                 <  1:486112/111‑1 (MQ=255)
                                                                                                                                                              aTGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTATGGATTACTGTTGCTTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTc                                                                                                                      >  1:252736/1‑122 (MQ=255)
                                                                                                                                                                gACCANCCATTNTNNAGTATGTCCTCGCCGTCGANTCCATCCTTTGTATTNTCTGNATTACNGt                                                                                                                                                                              >  2:196594/1‑64 (MQ=255)
                                                                                                                                                                gACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCGTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCGCTTTCAGCTACGCCTCCCGCCAGCGGATGCAAAACCGTGCGCGTttatta                                                                         <  1:196594/165‑1 (MQ=255)
                                                                                                                                                                   naGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCAGTTTTCTTCATTTGCCAGCGCCCCCACATCACATCTCTCTTTAAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGCAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGgc  >  1:52990/2‑233 (MQ=255)
                                                                                                                                                                   cAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTTCTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATtgtg                            <  1:21865/207‑1 (MQ=255)
                                                                                                                                                                   cAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGNNAGCGccnnc                                                                                                                                        >  2:21865/1‑99 (MQ=255)
                                                                                                                                                                    aGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTAATGTTGATTTGCTTGCCGTTTTCTTCATTTGCCAGCGCCACAACATCACATATCTCTTTCAGCTAAGCCGCACGCCAGCGGATGCAAAACCGTGCGCGTTTATTACAACAGTACCAAAc                                                           >  1:501914/1‑175 (MQ=255)
                                                                                                                                                                      ccATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTAcgc                                                                                                              <  1:113926/122‑1 (MQ=255)
                                                                                                                                                                      ccATTATTGAGAATGTCCTTGCCGTCNTTTGCGNNCTNTGTCTTGTCTGGATTAATGNTGNNNTNNNTNATGTTNTCNNCATTTGCCANCGCCACCACATCACATATCTCTTTCAGCTAcgc                                                                                                              >  2:113926/1‑122 (MQ=255)
                                                                                                                                                                              gAGTATGTCCTTGCCGTCGATTCCATCCTTTTTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCATCGCCACCCCATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATg                                                                                              <  1:407923/130‑1 (MQ=255)
                                                                                                                                                                              gAGTATGTCCTCGCCGTCGATTCCATCCTTTCTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTttcttcct                                                                                                                                                      >  1:180979/1‑72 (MQ=255)
                                                                                                                                                                                gTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGAtt                                                                                                                                                                                   <  1:4271/43‑1 (MQ=255)
                                                                                                                                                                                gTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGAtt                                                                                                                                                                                   >  2:4271/1‑43 (MQ=255)
                                                                                                                                                                                  nnnnnnttNNNNNCGANTCCATCCTTGTTATTGCCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTaaaa                                                                      <  2:314262/144‑1 (MQ=255)
                                                                                                                                                                                     tCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACCTCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTaaa                                                                       >  1:31048/1‑146 (MQ=255)
                                                                                                                                                                                           ccGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTg                                                                                                                                                   <  1:290160/64‑1 (MQ=255)
                                                                                                                                                                                             gTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAgcgcn                                                                                                                                           <  1:302785/70‑2 (MQ=255)
                                                                                                                                                                                             gTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGcc                                                                                                                                           >  2:302785/1‑70 (MQ=255)
                                                                                                                                                                                               ngATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTttcttc                                                                                                                                                        >  1:315359/2‑55 (MQ=255)
                                                                                                                                                                                               cGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTttcttc                                                                                                                                                        <  2:315359/55‑1 (MQ=255)
                                                                                                                                                                                                    ccATCCTTTGGATTGTCTGGATTACTGTTGATTTGTTTGCCGTTttcttc                                                                                                                                                        <  1:431568/50‑1 (MQ=255)
                                                                                                                                                                                                        ccTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACAGATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGc                                             <  1:237070/153‑1 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                      
GCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGC  >  NC_000913/246534‑246929

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