Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 693,030 C→A 8.0% V117L (GTA→TTA)  ybeZ ← heat shock protein, putative NTPase; PhoH‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913693,0300CA 8.0% 173.4 / 5.5 50V117L (GTA→TTA) ybeZheat shock protein, putative NTPase; PhoH‑like protein
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base A (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATCTGTTCCGGTTCGATATCCTGAATCTGACCGCGCATCGGGGCA  >  NC_000913/692896‑693169
                                                                                                                                      |                                                                                                                                           
cTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGtt                                                                                                                                <  2:321061/148‑1 (MQ=255)
    aGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGtt                                                                                                                                     >  1:653075/1‑139 (MQ=255)
    aGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGtt                                                                                                                                     <  2:653075/139‑1 (MQ=255)
     gggCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGtt                                                                                                                                     >  4:843359/1‑138 (MQ=255)
     gggCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGtt                                                                                                                                     >  4:431634/1‑138 (MQ=255)
     gggCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGtt                                                                                                                                     <  3:843359/138‑1 (MQ=255)
     gggCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGtt                                                                                                                                     <  3:431634/138‑1 (MQ=255)
                   gCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGACCAGCGCCGAAGGTAATGTCATAGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACTCCGCGTTTGCTTTTGATATAGAGCggaagag                                                                                                             >  3:705117/1‑142 (MQ=255)
                        cACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTGATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGtt                                                                                                                                     >  4:87296/1‑119 (MQ=255)
                        cACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGtt                                                                                                                                     <  3:87296/119‑1 (MQ=255)
                                   aggtTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGt                                                                                                                                      >  1:572647/1‑107 (MQ=255)
                                   aggtTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGt                                                                                                                                      <  2:572647/107‑1 (MQ=255)
                                      ttttACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGttt                                                                                                                                    >  2:105748/1‑106 (MQ=255)
                                      ttttACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGttt                                                                                                                                    <  1:105748/106‑1 (MQ=255)
                                      ttttACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGttt                                                                                                                               >  3:196206/1‑111 (MQ=255)
                                      ttttACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGttt                                                                                                                               <  4:196206/111‑1 (MQ=255)
                                        ttACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCcg                                                                                        >  1:718629/1‑148 (MQ=255)
                                                      cgGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGt                                                                          >  4:812226/1‑148 (MQ=255)
                                                             aCGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCcgcg                                                                   >  3:249624/1‑148 (MQ=255)
                                                                  gAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCcgcg                                                                   >  2:991152/1‑143 (MQ=255)
                                                                  gAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCcgcg                                                                   <  1:991152/143‑1 (MQ=255)
                                                                  gAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGcttct                                                              <  3:812226/148‑1 (MQ=255)
                                                                       tAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGat                                                                                                                           <  4:535405/82‑1 (MQ=255)
                                                                       tAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGat                                                                                                                           >  3:535405/1‑82 (MQ=255)
                                                                       tAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTaa                                                          >  4:628144/1‑147 (MQ=255)
                                                                       tAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTaa                                                          <  3:628144/147‑1 (MQ=255)
                                                                           gTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCa                                                     <  4:249624/148‑1 (MQ=255)
                                                                             cATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTg                                                                                                                             <  2:806697/74‑1 (MQ=255)
                                                                             cATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTg                                                                                                                             >  1:806697/1‑74 (MQ=255)
                                                                                ggTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTg                                                <  2:718629/148‑1 (MQ=255)
                                                                                   cGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCaa                                                    >  1:298799/1‑141 (MQ=255)
                                                                                   cGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCaa                                                    <  2:298799/141‑1 (MQ=255)
                                                                                          atTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGa                                                                                                                      <  3:181561/68‑1 (MQ=255)
                                                                                          atTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGa                                                                                                                      >  4:181561/1‑68 (MQ=255)
                                                                                            tGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGa                                              <  3:1111519/138‑1 (MQ=255)
                                                                                            tGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGa                                              >  4:1111519/1‑138 (MQ=255)
                                                                                                 aTGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTa                                                           >  3:585433/1‑120 (MQ=255)
                                                                                                 aTGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTa                                                           <  4:585433/120‑1 (MQ=255)
                                                                                                     aCTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAAt                                                         <  1:476994/118‑1 (MQ=255)
                                                                                                     aCTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAAt                                                         >  2:476994/1‑118 (MQ=255)
                                                                                                                        acgcgGCGTAATTAAGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATCTGTTCCGGTTCGATATCCTGAATCTGACCGCGCATc        >  3:810500/2‑148 (MQ=255)
                                                                                                                         cgcgGCGTAATTAAGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAAt                                                         <  1:567288/98‑1 (MQ=255)
                                                                                                                         cgcgGCGTAATTAAGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAAt                                                         >  2:567288/1‑98 (MQ=255)
                                                                                                                              cgtAATTAAGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATCTGTTCCGGTTCGATATCCTGAATCTGACCGCGCATCGGGGCa  <  4:810500/146‑1 (MQ=255)
                                                                                                                                tAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATCTGTTCCGGTTCGATATCCTGAATCTGACCGCGCATCg       <  4:476482/141‑1 (MQ=255)
                                                                                                                                tAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATCTGTTCCGGTTCGATATCCTGAATCTGACCGCGCATCg       >  3:476482/1‑141 (MQ=255)
                                                                                                                                   ttACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCcg                                                                     <  2:300198/76‑1 (MQ=255)
                                                                                                                                   ttACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCcg                                                                     >  1:300198/1‑76 (MQ=255)
                                                                                                                                   ttACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGcttct                                                              >  1:691271/1‑83 (MQ=255)
                                                                                                                                   ttACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGcttct                                                              <  2:691271/83‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                           
CTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGCGTTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGGCACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATCTGTTCCGGTTCGATATCCTGAATCTGACCGCGCATCGGGGCA  >  NC_000913/692896‑693169

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