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

AGCGCAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCG  >  NC_000913/3301463‑3301725
                                                                                                                                        |                                                                                                                              
aGCGCAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGc                                                                                                                              <  1:494825/139‑1 (MQ=255)
                 aaaTTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGaa                                                                                                             >  2:1148805/1‑139 (MQ=255)
                              gctCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGaa                                                                                                >  2:720560/1‑139 (MQ=255)
                                         ccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAGGCGGTGAGttt                                                                                     >  1:1078339/1‑139 (MQ=255)
                                           ggAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                   <  1:1636857/139‑1 (MQ=255)
                                                       gcgcTTAAGGCGGCCATCGACAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAACCGGTGAGTTTTGTCCATCAACATc                                                                       >  1:1603297/1‑139 (MQ=255)
                                                          cTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgcc                                                                    <  1:1627966/139‑1 (MQ=255)
                                                              aGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCaa                                                                <  2:1143402/139‑1 (MQ=255)
                                                                        cGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATc                                                      <  2:642943/139‑1 (MQ=255)
                                                                        cGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATc                                                      <  2:1078339/139‑1 (MQ=255)
                                                                         gAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATcg                                                     <  1:1281702/139‑1 (MQ=255)
                                                                                     gATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACtt                                         >  2:180991/1‑139 (MQ=255)
                                                                                        gCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTg                                      <  1:1148805/139‑1 (MQ=255)
                                                                                                      gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGtttt                                                                                    >  2:1270917/1‑79 (MQ=255)
                                                                                                      gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGtttt                                                                                    <  1:1270917/79‑1 (MQ=255)
                                                                                                      gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGt                        >  2:306926/1‑139 (MQ=255)
                                                                                                      gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGt                        >  1:891159/1‑139 (MQ=255)
                                                                                                        tAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTa                      <  1:720560/139‑1 (MQ=255)
                                                                                                               tgatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGt               >  1:445601/1‑139 (MQ=255)
                                                                                                                gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGc                                                                                                     <  2:1627557/52‑1 (MQ=255)
                                                                                                                gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGc                                                                                                     >  1:1627557/1‑52 (MQ=255)
                                                                                                                gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCg                              <  2:113085/123‑1 (MQ=255)
                                                                                                                gatACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCg                              >  1:113085/1‑123 (MQ=255)
                                                                                                                  tACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTgg            >  1:674498/1‑139 (MQ=255)
                                                                                                                        tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcgc      >  2:278104/1‑139 (MQ=255)
                                                                                                                            cGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCg  >  1:219140/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                              
AGCGCAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCG  >  NC_000913/3301463‑3301725

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