Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,923,946 (G)6→7 coding (646/873 nt) yebZ ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,923,9401.G94.0% 173.0 / 3.6 50coding (652/873 nt)yebZputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (24/23);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CCAGGGTGCCGCACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAACGGCCTTTGATGAGTT  >  NC_000913/1923806‑1924074
                                                                                                                                       |                                                                                                                                      
ccAGGGTGCCGCACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTggggg                                                                                                                                   <  2:1027945/140‑1 (MQ=255)
  aGGGTGCCGCACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTggggggg                                                                                                                                 <  1:247375/140‑2 (MQ=255)
       gCCGCACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATc                                                                                                                             <  2:355345/140‑1 (MQ=255)
        ccGCACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATcc                                                                                                                            >  1:260653/1‑140 (MQ=255)
         cGCACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCg                                                                                                                           >  2:161233/1‑140 (MQ=255)
           cACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATc                                                                                                                             >  1:447514/1‑136 (MQ=255)
           cACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATc                                                                                                                             <  2:447514/136‑1 (MQ=255)
                tCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATc                                                                                                                             >  1:951861/1‑131 (MQ=255)
                tCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATc                                                                                                                             <  2:951861/131‑1 (MQ=255)
                   tCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATcaaca                                                                                                                 >  1:206515/1‑140 (MQ=255)
                   tCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATcaaca                                                                                                                 <  1:287931/140‑1 (MQ=255)
                   tCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATcaaca                                                                                                                 <  1:330072/140‑1 (MQ=255)
                     aTCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAAt                                                                                                               <  2:612019/140‑1 (MQ=255)
                          gCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCAtt                                                                                                          >  1:40790/1‑140 (MQ=255)
                           ccTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTa                                                                                                         >  1:379556/1‑140 (MQ=255)
                           ccTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTa                                                                                                         >  1:583705/1‑140 (MQ=255)
                              cATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGc                                                                                                      <  2:1047672/140‑1 (MQ=255)
                                      gTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCt                                                                                              <  2:1053319/140‑1 (MQ=255)
                                           aGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCa                                                                                         <  2:710858/140‑1 (MQ=255)
                                                 tAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGc                                                                                   >  2:165756/1‑140 (MQ=255)
                                                 tAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGc                                                                                   >  1:1007189/1‑140 (MQ=255)
                                                 tAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCAGTAAGGCGGAAGCGGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCg                                                                                 >  2:997346/1‑140 (MQ=255)
                                                                  aaTGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGc                                                                  <  2:312326/140‑1 (MQ=255)
                                                                    tGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGccc                                                                >  2:12058/1‑140 (MQ=255)
                                                                    tGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCg                                                               >  1:173597/1‑140 (MQ=255)
                                                                     gACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCg                                                               <  1:273902/140‑1 (MQ=255)
                                                                          ccaTAGTCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAg                                                         >  1:353712/1‑140 (MQ=255)
                                                                           caTAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAg                                                         >  1:292807/1‑140 (MQ=255)
                                                                                  atcaCCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAAc                                                  >  2:1001051/1‑140 (MQ=255)
                                                                                  atcaCCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAAc                                                  >  2:1051564/1‑140 (MQ=255)
                                                                                         ggATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAg                                           <  2:260653/140‑1 (MQ=255)
                                                                                          gATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAgc                                          >  1:754180/1‑140 (MQ=255)
                                                                                          gATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAgc                                          <  2:40790/140‑1 (MQ=255)
                                                                                            tGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAgcgc                                        <  2:292807/140‑1 (MQ=255)
                                                                                              gCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAg                                           >  1:179200/1‑135 (MQ=255)
                                                                                              gCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAg                                           <  2:179200/135‑1 (MQ=255)
                                                                                                    aGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAg                                >  1:83460/1‑140 (MQ=255)
                                                                                                         tAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTgg                           >  1:260945/1‑140 (MQ=255)
                                                                                                         tAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTgg                           >  1:702793/1‑140 (MQ=255)
                                                                                                         tAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTgg                           >  2:487782/1‑140 (MQ=255)
                                                                                                              gTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGgtg                         >  1:389063/1‑137 (MQ=255)
                                                                                                              gTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGgtg                         <  2:389063/137‑1 (MQ=255)
                                                                                                                 ggCCCCAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCa                   >  2:973295/1‑140 (MQ=255)
                                                                                                                   ccccAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGaa                                   <  2:104304/122‑1 (MQ=255)
                                                                                                                   ccccAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGaa                                   >  1:104304/1‑122 (MQ=255)
                                                                                                                   ccccAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAAc                 <  2:206515/140‑1 (MQ=255)
                                                                                                                     ccAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAACgg               >  2:793246/1‑140 (MQ=255)
                                                                                                                      cAGTAAGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCACCAGCGCCTGAATAGCCTGGTGTCGCCAACGGc              <  1:12058/140‑1 (MQ=255)
                                                                                                                           aGTCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAACGGCCTTtg         <  1:973295/140‑1 (MQ=255)
                                                                                                                             tCGTAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAACGGCCTTtgat       <  1:161233/140‑1 (MQ=255)
                                                                                                                               gTAAGCGTGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAACGGCCTTTGATGAg    <  2:353712/140‑1 (MQ=255)
                                                                                                                                tAAGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAACGGCCTTTGATGAg    <  1:1001051/140‑1 (MQ=255)
                                                                                                                                  aGCGTGGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAACGGCCTTTGATGAGtt  <  2:960656/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
CCAGGGTGCCGCACGATCTTCATCCTGCCTCATGCGCGGTACAAGAACATAACGATTAGCCAGCGCAATGACCACCATAATCATCACCAGGATGGCTTTGAGCAATAAGAGTTGGCCCCAGTAAGTCGTAAGCGTGGGGGGAAATCCGGTAATCAACAATGCATTAAGCACGCCGCTTGCCAGTACGCCGATCACCGCAAAATGCCCGCACCAGGAAAAACGCATCAGCGCCTGAATAGCCTGGTGTCGCCAACGGCCTTTGATGAGTT  >  NC_000913/1923806‑1924074

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