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% 132.3 / NA 41G39S (GGC→AGC) yhbUU32 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/18);  total (23/18)

GCAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGT  >  NC_000913/3301466‑3301726
                                                                                                                                     |                                                                                                                               
gcAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCtt                                                                                                                            >  2:1072552/1‑139 (MQ=255)
     cTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACtt                                                                                                                       >  2:1386279/1‑139 (MQ=255)
                         ctgctCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGc                                                                                                                               >  2:235186/1‑111 (MQ=255)
                         ctgctCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAgg                                                                                                   <  1:1386279/139‑1 (MQ=255)
                          tgctCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGa                                                                                                  >  1:1068444/1‑139 (MQ=255)
                                gccctgccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGt                                                                                            <  1:1072552/139‑1 (MQ=255)
                                     gccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGTCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGtt                                                                                       <  2:1068444/139‑1 (MQ=255)
                                      ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGcc                                                                                                                              >  2:453134/1‑99 (MQ=255)
                                      ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGcc                                                                                                                              <  1:453134/99‑1 (MQ=255)
                                      ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGc                                                                                                                               <  2:281031/98‑1 (MQ=255)
                                      ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCt                                                                                                                             <  2:865862/100‑1 (MQ=255)
                                      ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCt                                                                                                                             >  1:865862/1‑100 (MQ=255)
                                                       cTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGt                                                                                            <  2:1193430/116‑1 (MQ=255)
                                                       cTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGt                                                                                            >  1:1193430/1‑116 (MQ=255)
                                                        ttAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccg                                                                    >  2:1182521/1‑139 (MQ=255)
                                                               ggcCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACtt                                                             >  2:1189507/1‑139 (MQ=255)
                                                                              cgcAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGcc                                                                                                                              <  1:1148821/59‑1 (MQ=255)
                                                                              cgcAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGcc                                                                                                                              >  2:1148821/1‑59 (MQ=255)
                                                                                      cTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGc                                                                                                      <  2:301157/75‑1 (MQ=255)
                                                                                      cTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGc                                                                                                      >  1:301157/1‑75 (MQ=255)
                                                                                                cGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGttt                                                                                      >  1:72355/1‑81 (MQ=255)
                                                                                                cGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGttt                                                                                      <  2:72355/81‑1 (MQ=255)
                                                                                                   gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCg                                                                                                                  <  1:376809/50‑1 (MQ=255)
                                                                                                   gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCg                                                                                                                  >  2:376809/1‑50 (MQ=255)
                                                                                                   gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAgg                                                                                                   >  2:982728/1‑65 (MQ=255)
                                                                                                   gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAgg                                                                                                   <  1:982728/65‑1 (MQ=255)
                                                                                                   gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGt                         >  1:358425/1‑139 (MQ=255)
                                                                                                             gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                    <  1:1215426/70‑1 (MQ=255)
                                                                                                             gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                    >  2:1215426/1‑70 (MQ=255)
                                                                                                             gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccgc                                                                   >  2:241019/1‑87 (MQ=255)
                                                                                                             gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccgc                                                                   <  1:241019/87‑1 (MQ=255)
                                                                                                                     tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcgc       >  1:619366/1‑139 (MQ=255)
                                                                                                                     tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAg          <  2:159049/136‑1 (MQ=255)
                                                                                                                     tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAg          >  1:159049/1‑136 (MQ=255)
                                                                                                                         cGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCGACATCGCCGCAAACTTca                                                           <  2:837671/83‑1 (MQ=255)
                                                                                                                         cGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTca                                                           >  1:837671/1‑83 (MQ=255)
                                                                                                                          gTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGt  >  1:1320318/1‑139 (MQ=255)
                                                                                                                            cACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcgcg      <  2:1235497/133‑1 (MQ=255)
                                                                                                                            cACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcgcg      >  1:1235497/1‑133 (MQ=255)
                                                                                                                                 cGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgc         <  2:1510829/125‑1 (MQ=255)
                                                                                                                                 cGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgc         >  1:1510829/1‑125 (MQ=255)
                                                                                                                                     |                                                                                                                               
GCAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGT  >  NC_000913/3301466‑3301726

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: ‑