Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,301,599 G→A G39S (GGC→AGC)  yhbU → U32 peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,301,5990GA100.0% 94.1 / NA 31G39S (GGC→AGC) yhbUU32 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/14);  total (17/14)

AGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGTGGATATG  >  NC_000913/3301475‑3301733
                                                                                                                            |                                                                                                                                      
aGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTAcc                                                                                                                          <  2:850432/139‑1 (MQ=255)
                             ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGaa                                                                                                        <  2:422474/128‑1 (MQ=255)
                             ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGaa                                                                                                        >  1:422474/1‑128 (MQ=255)
                             ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGttt                                                                                             >  2:24999/1‑139 (MQ=255)
                                              cTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgcc                                                                            >  2:170479/1‑139 (MQ=255)
                                              cTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACGTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgcc                                                                            <  2:854879/139‑1 (MQ=255)
                                                tAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccgc                                                                          >  1:683368/1‑139 (MQ=255)
                                                  aGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCaa                                                                        <  1:170479/139‑1 (MQ=255)
                                                            cGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATc                                                              <  1:24999/139‑1 (MQ=255)
                                                                     cgcAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGAt                                                         >  1:774163/1‑135 (MQ=255)
                                                                     cgcAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGAt                                                         <  2:774163/135‑1 (MQ=255)
                                                                                          gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACt                                                                                                                               >  1:63333/1‑44 (MQ=255)
                                                                                          gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACt                                                                                                                               <  2:63333/44‑1 (MQ=255)
                                                                                          gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGt                                >  2:846752/1‑139 (MQ=255)
                                                                                                    gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGtt                      >  1:377750/1‑139 (MQ=255)
                                                                                                        ccAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGttt                                                                                             <  2:902179/64‑1 (MQ=255)
                                                                                                        ccAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGttt                                                                                             >  1:902179/1‑64 (MQ=255)
                                                                                                           aTGCCCGTCACTTCGCCAGCCTTACCTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTTTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcg               <  2:278248/139‑1 (MQ=255)
                                                                                                           aTGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCgg                                                                                                    >  2:981690/1‑54 (MQ=255)
                                                                                                           aTGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCgg                                                                                                    <  1:981690/54‑1 (MQ=255)
                                                                                                            tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGttt                                                                                             >  2:813826/1‑60 (MQ=255)
                                                                                                            tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGttt                                                                                             <  1:813826/60‑1 (MQ=255)
                                                                                                             gCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcgcg             >  1:941809/1‑139 (MQ=255)
                                                                                                                cGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGc                                                                                                             <  2:665922/40‑1 (MQ=255)
                                                                                                                cGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGc                                                                                                             >  1:665922/1‑40 (MQ=255)
                                                                                                                cGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCg          >  2:63550/1‑139 (MQ=255)
                                                                                                                   cACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGcc           >  1:968953/1‑135 (MQ=255)
                                                                                                                   cACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGcc           <  2:968953/135‑1 (MQ=255)
                                                                                                                      ttCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATcg                                                             >  2:238221/1‑82 (MQ=255)
                                                                                                                      ttCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATcg                                                             <  1:238221/82‑1 (MQ=255)
                                                                                                                        cGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGTGGATATg  >  1:7862/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                      
AGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGTGGATATG  >  NC_000913/3301475‑3301733

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