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,5990GA98.0% 162.7 / NA 51G39S (GGC→AGC) yhbUU32 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/23);  minor base C (1/0);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGTGGATATG  >  NC_000913/3301467‑3301733
                                                                                                                                    |                                                                                                                                      
cAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTa                                                                                                                                  >  2:693507/1‑139 (MQ=255)
   ccTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACt                                                                                                                               >  1:721797/1‑139 (MQ=255)
         ggAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCg                                                                                                                         <  1:1185327/139‑1 (MQ=255)
                          gctCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGaa                                                                                                        <  1:457379/139‑1 (MQ=255)
                             ctgccctgccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAg                                                                                                           >  2:95409/1‑133 (MQ=255)
                             ctgccctgccGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAg                                                                                                           <  1:95409/133‑1 (MQ=255)
                                          aaTCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTGAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCgaga                                                                                                                      >  2:1641153/1‑109 (MQ=255)
                                          aaTCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTGAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCgaga                                                                                                                      <  1:1641153/109‑1 (MQ=255)
                                                        tAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                           <  1:1635807/122‑1 (MQ=255)
                                                        tAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                           >  2:1635807/1‑122 (MQ=255)
                                                                 cATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgcc                                                                            <  2:414489/128‑1 (MQ=255)
                                                                 cATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgcc                                                                            >  1:414489/1‑128 (MQ=255)
                                                                    cGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATc                                                              <  2:721797/139‑1 (MQ=255)
                                                                      aaaaCGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATcgc                                                            >  1:621442/1‑139 (MQ=255)
                                                                                      tGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTgcg                                            <  1:693507/139‑1 (MQ=255)
                                                                                          tataTCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATc                                        >  2:534483/1‑139 (MQ=255)
                                                                                              tCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGa                                    <  2:621442/139‑1 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAgg                                                                                                          >  2:1332484/1‑65 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAgg                                                                                                          <  1:1332484/65‑1 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAg                                                                                                >  2:1122719/1‑75 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAg                                                                                                <  1:1122719/75‑1 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                           >  1:1003566/1‑80 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                           <  2:1003566/80‑1 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGt                                <  1:810982/139‑1 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGt                                >  2:323051/1‑139 (MQ=255)
                                                                                                  gCTAAAAGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGt                                >  2:576031/1‑139 (MQ=255)
                                                                                                      aaaGATGATACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACg                            <  1:534483/139‑1 (MQ=255)
                                                                                                              tACCAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTgg                    >  2:1258737/1‑139 (MQ=255)
                                                                                                                ccAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGaa                                                                                                                     <  1:1479510/40‑1 (MQ=255)
                                                                                                                ccAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGaa                                                                                                                     >  2:1479510/1‑40 (MQ=255)
                                                                                                                ccAATGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCa                  >  2:942047/1‑139 (MQ=255)
                                                                                                                   aTGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcg               <  1:576031/139‑1 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACAt                                                                                >  1:551874/1‑73 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACAt                                                                                <  2:551874/73‑1 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcg                                                                              >  2:1085579/1‑75 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcg                                                                              <  1:1085579/75‑1 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGAt                                                         <  2:1327199/96‑1 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGAt                                                         >  1:1327199/1‑96 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCAt                                         >  1:518541/1‑112 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCAt                                         <  2:518541/112‑1 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTa                              <  2:1443587/123‑1 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTa                              >  1:1443587/1‑123 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcgc              >  2:1052410/1‑139 (MQ=255)
                                                                                                                    tGCCCGTCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAgcgc              >  2:113732/1‑139 (MQ=255)
                                                                                                                         gTCACTTCTCCCGCCTTCACTCTTGCTAGAAAAAATTGCATGAGGCGGCGGGATTTGTCCATCAGCATCGCCGCAAACTTCACATCGCGGTTAACACTTTTGCGCATCCGGACAGTTACGCGCCGTGGCAGCGAGACGt         >  1:1360638/1‑139 (MQ=255)
                                                                                                                          tCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                           <  2:405213/56‑1 (MQ=255)
                                                                                                                          tCACTTCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTg                                                                                           >  1:405213/1‑56 (MQ=255)
                                                                                                                              ttCGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGTGGata    >  2:563271/1‑139 (MQ=255)
                                                                                                                                cGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccgc                                                                          <  2:1061396/67‑1 (MQ=255)
                                                                                                                                cGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATcgccgc                                                                          >  1:1061396/1‑67 (MQ=255)
                                                                                                                                cGCCAGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGTGGATATg  >  2:289202/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                      
CAACCTTCAGGAAAAATTATGGAGCTGCTCTGCCCTGCCGGAAATCTCCCGGCGCTTAAGGCGGCCATCGAAAACGGCGCAGATGCTGTTTATATCGGGCTAAAAGATGATACCAATGCCCGTCACTTCGCCGGCCTTAACTTTACCGAGAAAAAATTGCAGGAAGCGGTGAGTTTTGTCCATCAACATCGCCGCAAACTTCACATCGCGATTAACACTTTTGCGCATCCGGACGGTTACGCCCGTTGGCAGCGCGCCGTGGATATG  >  NC_000913/3301467‑3301733

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