Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,160,278 G→A S171S (AGC→AGT sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,160,2780GA100.0% 83.9 / NA 25S171S (AGC→AGTsthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/12);  total (13/12)

ATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCG  >  NC_000913/4160143‑4160402
                                                                                                                                       |                                                                                                                            
aTCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTg                                                                                                                           >  1:354313/1‑139 (MQ=255)
          gCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGt                                                                                                                 <  1:124113/139‑1 (MQ=255)
           ccAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTc                                                                                                                <  2:221658/139‑1 (MQ=255)
               cagGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTg                                                                                                            <  2:82291/139‑1 (MQ=255)
                                                    tCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTg                                                                       <  2:354313/139‑1 (MQ=255)
                                                       aTACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTc                                                                                                                >  1:254248/1‑95 (MQ=255)
                                                       aTACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTc                                                                                                                <  2:254248/95‑1 (MQ=255)
                                                       aTACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgat                                                                    >  2:111010/1‑139 (MQ=255)
                                                       aTACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgat                                                                    >  2:371099/1‑139 (MQ=255)
                                                       aTACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgat                                                                    >  1:201030/1‑139 (MQ=255)
                                                                                     ccGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTGGCTGTCGTAAATGCGTGGATGGGTGAAATCAAc                                                                               >  1:245719/1‑98 (MQ=255)
                                                                                     ccGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAAc                                                                               <  2:245719/98‑1 (MQ=255)
                                                                                        aTCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGc                                                                                                              <  1:242168/64‑1 (MQ=255)
                                                                                        aTCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGc                                                                                                              >  2:242168/1‑64 (MQ=255)
                                                                                        aTCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGga                                                                     <  1:449429/105‑1 (MQ=255)
                                                                                        aTCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGga                                                                     >  2:449429/1‑105 (MQ=255)
                                                                                        aTCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAAtt                                   >  2:248039/1‑139 (MQ=255)
                                                                                                                 cATGGCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCt          <  1:248039/139‑1 (MQ=255)
                                                                                                                    ggCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGga                                                                     <  1:394214/77‑1 (MQ=255)
                                                                                                                    ggCGCGGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGga                                                                     >  2:394214/1‑77 (MQ=255)
                                                                                                                      cgcgGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATc                                                                                  <  2:163383/62‑1 (MQ=255)
                                                                                                                      cgcgGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATc                                                                                  >  1:163383/1‑62 (MQ=255)
                                                                                                                        cgGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGga                                                                     >  2:21081/1‑73 (MQ=255)
                                                                                                                        cgGTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGga                                                                     <  1:21081/73‑1 (MQ=255)
                                                                                                                         ggTTCGTGGTGCATACTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCg  >  2:193036/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                            
ATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCCG  >  NC_000913/4160143‑4160402

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