Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,159,844 G→A P435S (CCA→TCA)  ygiQ ← Radical SAM superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,159,8440GA100.0% 145.2 / NA 43P435S (CCA→TCA) ygiQRadical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/22);  total (21/22)

GATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGGCCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTCGGTGATAGAACAGAAA  >  NC_000913/3159712‑3159974
                                                                                                                                    |                                                                                                                                  
gATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACccacc                                                                                                                              <  2:80661/139‑1 (MQ=255)
                 ggttcgtggttcgtGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGATCCACCAAGATCGGAAATCACGc                                                                                                             >  2:116380/1‑139 (MQ=255)
                     cgtggttcgtgTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCAAGGCGATTTGCAGCGCCACATATACCTGTTGTCAGTTGACGCACCAAGATCTGAAATCACGCCAGt                                                                                                         >  2:260736/1‑139 (MQ=255)
                           tcgtgTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAAcc                                                                                                   <  1:47123/139‑1 (MQ=255)
                              tgtCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTgg                                                                                                >  1:453709/1‑139 (MQ=255)
                                tCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGaa                                                                                              <  1:403637/139‑1 (MQ=255)
                                        cGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTcg                                                                                      >  1:760266/1‑139 (MQ=255)
                                                      gATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGATGACCCACCAAGATCAGAAATCACGCACGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGAt                                                                        <  1:513462/139‑1 (MQ=255)
                                                       aTAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAAcc                                                                                                   >  2:23757/1‑111 (MQ=255)
                                                       aTAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAAcc                                                                                                   <  1:23757/111‑1 (MQ=255)
                                                       aTAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGAtc                                                                       >  2:27128/1‑139 (MQ=255)
                                                         aaaCGCACGACAAACGCCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGAtctc                                                                     >  2:141737/1‑139 (MQ=255)
                                                                 gACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGa                                                             >  2:754221/1‑139 (MQ=255)
                                                                  aCAAACGGCGAGAAGTTTGTGTAGCGCGGGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACACACCAAGATCGTAAATCACGCCCGTAAAACCTGGAACGGTGTGGCGGATCGCTTCGATCTCATTAATGAt                                                            <  2:708444/139‑1 (MQ=255)
                                                                     aaCGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGa                                                         <  1:27128/139‑1 (MQ=255)
                                                                            aCAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCg                                                  >  2:280062/1‑139 (MQ=255)
                                                                               aGTTTGTTCAGCGCGTGGCGATTTGCATCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGaa                                               <  2:23183/139‑1 (MQ=255)
                                                                               aGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGaa                                               <  1:324612/139‑1 (MQ=255)
                                                                               aGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGaa                                               <  2:95266/139‑1 (MQ=255)
                                                                                   tGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGc                                           <  1:14391/139‑1 (MQ=255)
                                                                                          cgcgTGGCGATTTGTAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGaat                                    >  1:491882/1‑139 (MQ=255)
                                                                                          cgcgTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTa                                                                 >  2:645763/1‑110 (MQ=255)
                                                                                          cgcgTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTa                                                                 <  1:645763/110‑1 (MQ=255)
                                                                                          cgcgTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGaat                                    >  2:348429/1‑139 (MQ=255)
                                                                                          cgcgTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGaat                                    >  1:316491/1‑139 (MQ=255)
                                                                                                gCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATgcg                              <  1:280062/139‑1 (MQ=255)
                                                                                                gCGATTTGCAGCGCAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATgcg                              <  1:394662/139‑1 (MQ=255)
                                                                                                            gcAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGTTCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTc                  <  1:754221/139‑1 (MQ=255)
                                                                                                            gcAACATATACATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTc                  <  1:348429/139‑1 (MQ=255)
                                                                                                              aaCATATACATGTTGGCAGTTGATCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTCgg                <  1:116380/139‑1 (MQ=255)
                                                                                                                    taCATGTTGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTCGGTGATAg          <  1:141737/139‑1 (MQ=255)
                                                                                                                          ttGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGaa                                                                                              >  2:445072/1‑49 (MQ=255)
                                                                                                                          ttGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGaa                                                                                              <  1:445072/49‑1 (MQ=255)
                                                                                                                          ttGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTCGGTGATAGAACAGa    >  1:727051/1‑139 (MQ=255)
                                                                                                                           tGGCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTCGGTGATAGAACAGaa   >  2:71753/1‑139 (MQ=255)
                                                                                                                            ggCAGTTGACCCACCAAGATCGGAAATCACGCCCGTaaa                                                                                                      >  1:748302/1‑39 (MQ=255)
                                                                                                                            ggCAGTTGACCCACCAAGATCGGAAATCACGCCCGTaaa                                                                                                      <  2:748302/39‑1 (MQ=255)
                                                                                                                            ggCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACg                                                                                            >  2:670869/1‑49 (MQ=255)
                                                                                                                            ggCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACg                                                                                            <  1:670869/49‑1 (MQ=255)
                                                                                                                            ggCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTCGGTGATAGAACAGaaa  >  1:367131/1‑139 (MQ=255)
                                                                                                                            ggCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTCGGTGATAGAACAGaaa  >  1:190081/1‑139 (MQ=255)
                                                                                                                             gCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTc                  >  1:709286/1‑122 (MQ=255)
                                                                                                                             gCAGTTGACCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTc                  <  2:709286/122‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
GATAGAGGTTGATCGTTGGTTCGTGGTTCGTGTCCATGTGCGGACAAATATCCGGATAAACGCACGACAAACGGCGACAAGTTTGTTCAGCGCGTGGCGATTTGCAGCGCAACATATACATGTTGGCAGTTGGCCCACCAAGATCGGAAATCACGCCCGTAAAACCTGGAACGGTGTCGCGGATCGCTTCGATCTCATTAATGATCGAATCTTCGGAACGGCTCTGAATAATGCGCCCTTCGTGCTCGGTGATAGAACAGAAA  >  NC_000913/3159712‑3159974

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