Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,016,166 A→G Y36C (TAC→TGC)  ygfK → putative Fe‑S subunit oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,016,1660AG98.0% 177.3 / ‑4.0 50Y36C (TAC→TGC) ygfKputative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (32/17);  total (33/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

CTCAATGATATCTGTATAAGCTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTACTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGAACTAAAAACC  >  NC_000913/3016026‑3016311
                                                                                                                                            |                                                                                                                                                 
cTCAATGATATCTGTATAAGCTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGt                                                                                                                                          >  2:24778/1‑150 (MQ=255)
 tCAATGATATCTATATAAGCTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTa                                                                                                                                         >  1:722393/1‑150 (MQ=255)
    aTGATATCTGTATAAGCTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTaaa                                                                                                                                       >  1:76358/1‑149 (MQ=255)
       atatCTGTATAAGCTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAgg                                                                                                                                    >  2:310176/1‑149 (MQ=255)
             gTATAAGCTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACt                                                                                                                             <  1:609417/150‑1 (MQ=255)
                  aGCTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGc                                                                                                                       >  2:55023/1‑151 (MQ=255)
                    cTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGt                                                                                                                     >  2:312832/1‑151 (MQ=255)
                          agagGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCgg                                                                                                               >  1:633907/1‑151 (MQ=255)
                               gTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGaaa                                                                                                          <  2:223157/151‑1 (MQ=255)
                               gTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGaa                                                                                                           >  2:554159/1‑150 (MQ=255)
                                ttATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAAc                                                                                                         >  2:411679/1‑151 (MQ=255)
                                 tATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAAcc                                                                                                        >  2:542230/1‑151 (MQ=255)
                                   tGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCtg                                                                                                      >  2:155237/1‑151 (MQ=255)
                                    gggggATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCtg                                                                                                      <  1:594944/150‑1 (MQ=255)
                                     ggggATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCtgtg                                                                                                    <  2:101960/151‑1 (MQ=255)
                                        gATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCa                                                                                                 <  2:91067/151‑1 (MQ=255)
                                         atatTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCAc                                                                                                >  2:457297/1‑151 (MQ=255)
                                            ttatGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTc                                                                                              >  1:701946/1‑150 (MQ=255)
                                                  gTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCgg                                                                                       >  2:235610/1‑151 (MQ=255)
                                                  gTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCgg                                                                                       >  1:719385/1‑151 (MQ=255)
                                                  gTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCgg                                                                                       >  1:313071/1‑151 (MQ=255)
                                                             tttGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCCCCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCt                                                                             <  2:524917/150‑1 (MQ=255)
                                                             tttGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCt                                                                             <  1:204577/150‑1 (MQ=255)
                                                             tttGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTACTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTg                                                                            >  2:651700/1‑151 (MQ=255)
                                                                     aCTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGcac                                                                    <  2:348673/151‑1 (MQ=255)
                                                                      cTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGcac                                                                    >  2:747974/1‑150 (MQ=255)
                                                                      cTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGcac                                                                    >  2:548538/1‑150 (MQ=255)
                                                                           gACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGc                                                                >  1:727367/1‑149 (MQ=255)
                                                                           gACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCa                                                               <  1:522019/150‑1 (MQ=255)
                                                                                    atTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCa                                                      >  1:744000/1‑150 (MQ=255)
                                                                                         atgaATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTcgc                                                         >  2:221052/1‑142 (MQ=255)
                                                                                         atgaATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTcgc                                                         <  1:221020/142‑1 (MQ=255)
                                                                                           gaATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATAtt                                               <  2:17777/150‑1 (MQ=255)
                                                                                           gaATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTg                                              <  1:329031/151‑1 (MQ=255)
                                                                                               aCCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAACCTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCa                                           <  1:344800/150‑1 (MQ=255)
                                                                                                 caacaacaaCGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACtt                                        >  2:717209/1‑151 (MQ=255)
                                                                                                        aacGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCt                                 >  1:264981/1‑151 (MQ=255)
                                                                                                         acGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTg                                >  2:422192/1‑151 (MQ=255)
                                                                                                           gCTCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGAc                              >  2:159117/1‑151 (MQ=255)
                                                                                                             tCAATCTTTGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTg                            >  1:317460/1‑151 (MQ=255)
                                                                                                                   tttGGTATTCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGAc                      >  1:60048/1‑151 (MQ=255)
                                                                                                                          ttCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCAt               >  2:184333/1‑151 (MQ=255)
                                                                                                                          ttCCCGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCAt               >  1:196201/1‑151 (MQ=255)
                                                                                                                              cGAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCGAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGaa           <  2:325253/151‑1 (MQ=255)
                                                                                                                               gAGCAACAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGAAc          <  1:372373/151‑1 (MQ=255)
                                                                                                                               gAGCAACAGTTTGGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACGCCCGTCGCCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGaa           <  1:21809/150‑1 (MQ=255)
                                                                                                                                   aaCAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGAACTaa       >  2:496237/1‑150 (MQ=255)
                                                                                                                                     cAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGAACTaaaaa    <  1:258475/151‑1 (MQ=255)
                                                                                                                                     cAGTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGAACTaaaa     >  2:257133/1‑150 (MQ=255)
                                                                                                                                       gTTTTGCTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGAACTAAAAAcc  >  2:774455/1‑151 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CTCAATGATATCTGTATAAGCTAAGGAGAGGGTTATGGGGGATATTATGCGTCCCATTCCGTTTGAGGAACTTTTGACGCGCATATTTGATGAATACCAACAACAACGCTCAATCTTTGGTATTCCCGAGCAACAGTTTTACTCACCTGTAAAAGGTAAAACTGTTAGCGTCTTCGGTGAAACCTGTGCCACTCCCGTCGGCCCTGCCGCTGGCCCGCACACGCAGCTCGCGCAAAATATTGTCACTTCCTGGCTGACTGGCGGACGCTTCATCGAACTAAAAACC  >  NC_000913/3016026‑3016311

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