Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,171,879 C→T *111* (TAG→TAA ygiZ ← inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,171,8790CT100.0% 84.5 / NA 29*111* (TAG→TAAygiZinner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/15);  total (14/15)

TCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGAT  >  NC_000913/3171748‑3171998
                                                                                                                                   |                                                                                                                       
tCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                  <  1:649514/139‑1 (MQ=255)
tCCAGAGAATCGCCAAATTGCATGGCATGAATGGTGAATTGGGGAATGCGGAACAAGGTGGAATTGAGGCTAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGc                                                                                                                  <  2:748770/139‑1 (MQ=255)
 ccAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTTTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                 <  2:657147/139‑1 (MQ=255)
 ccAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                 <  1:342958/139‑1 (MQ=255)
     cGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAat                                                                                                             >  1:7189/1‑139 (MQ=255)
       aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAatat                                                                                                           >  2:145755/1‑139 (MQ=255)
       aaTCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAatat                                                                                                           >  1:702963/1‑139 (MQ=255)
                         catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                                                         >  1:870385/1‑139 (MQ=255)
                         catgAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAg                                                                                         >  1:837990/1‑139 (MQ=255)
                              aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTTCCCaaaa                                                                                            >  2:827663/1‑131 (MQ=255)
                              aTGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCaaaa                                                                                            <  1:827663/131‑1 (MQ=255)
                                          ggAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCa                                                                        <  2:705747/139‑1 (MQ=255)
                                                       aGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                           <  2:734573/139‑1 (MQ=255)
                                                        ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTaa                                                                                                                      >  1:18604/1‑79 (MQ=255)
                                                        ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTaa                                                                                                                      <  2:18604/79‑1 (MQ=255)
                                                        ggTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagc                                                          >  2:57829/1‑139 (MQ=255)
                                                           ggATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgagcgag                                                       <  2:702963/139‑1 (MQ=255)
                                                             aTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                          <  2:375646/118‑1 (MQ=255)
                                                             aTTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                          >  1:375646/1‑118 (MQ=255)
                                                                   ggCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAACCCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACata                                               <  1:783837/139‑1 (MQ=255)
                                                                       aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                 >  2:63381/1‑69 (MQ=255)
                                                                       aaGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCa                                                                                                                 <  1:63381/69‑1 (MQ=255)
                                                                                 ggTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGa                                 >  1:828898/1‑139 (MQ=255)
                                                                                         ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                           >  2:234681/1‑105 (MQ=255)
                                                                                         ctctGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACcgag                                                           <  1:234681/105‑1 (MQ=255)
                                                                                            tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                          >  2:489644/1‑87 (MQ=255)
                                                                                            tGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAAc                                                                          <  1:489644/87‑1 (MQ=255)
                                                                                                        ccGTTAAAACGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGattt          >  1:690377/1‑139 (MQ=255)
                                                                                                                aCGACTCTCATAGATTTTATTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTtaagat  <  1:57829/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                       
TCCAGCGAATCGCCAAATTGCATGGCATGAATGTTGAATTTGGGAATGCGGAACAAGGTGGATTTGAGGCGAAGGTAAGCTGGTAACTTCTCTGTCTATGAGAGCCGTTAAAACGACTCTCATAGATTTTACTAATAGCAAAATATAAACCGTCCCCAAAAAAGCCACCAACCACAAACCAGACAAATACCGAGCGAGAATACATATACCCAAAACTCGAACATTTCCCGCATAAAGAGTTTTCTTAAGAT  >  NC_000913/3171748‑3171998

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