Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 751,737 (G)6→7 coding (433/2448 nt) ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913751,7311.G95.0% 66.7 / ‑2.7 20coding (439/2448 nt)ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (7/12);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

GCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATACACCGATCTATCT  >  NC_000913/751594‑751862
                                                                                                                                          |                                                                                                                                   
gcgcCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTagag                                                                                                                                    >  1:170046/1‑139 (MQ=255)
gcgcCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTagag                                                                                                                                    <  2:369691/139‑1 (MQ=255)
                             tACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACg                                                                                                        <  1:747592/139‑1 (MQ=255)
                                 gTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACgccagc                                                                                                   <  2:170046/139‑1 (MQ=255)
                                       aTTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGc                                                                                              <  2:255156/139‑1 (MQ=255)
                                               tAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                                                      <  1:21715/139‑1 (MQ=255)
                                               tAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                                                      <  1:154558/139‑1 (MQ=255)
                                               tAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                                                      <  1:392238/139‑1 (MQ=255)
                                               tAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgggggggagc                                                                                                                           >  1:434987/1‑97 (MQ=255)
                                               tAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCGGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgggggggagc                                                                                                                           <  2:434987/101‑5 (MQ=255)
                                                                         cGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTg                                                                                           >  2:216259/1‑108 (MQ=255)
                                                                         cGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTg                                                                                           <  1:216259/108‑1 (MQ=255)
                                                                          ggTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCtt                                                           <  1:737596/139‑1 (MQ=255)
                                                                                      cAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGa                                               >  2:651894/1‑139 (MQ=255)
                                                                                                   aTCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                                                      <  1:259358/87‑1 (MQ=255)
                                                                                                   aTCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGaa                                                                                      >  2:259358/1‑87 (MQ=255)
                                                                                                                ccGGCAACACCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGaaa                     >  1:245133/1‑139 (MQ=255)
                                                                                                                       cacCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTc                                     <  1:678355/116‑1 (MQ=255)
                                                                                                                       cacCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTc                                     >  2:678355/1‑116 (MQ=255)
                                                                                                                       cacCTTCTTTCCATGTAGAGGGGGGGAGCCCAGTTTTCGGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATac              <  2:251835/139‑1 (MQ=255)
                                                                                                                           ttcttTCCATGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATACACCg          >  1:24280/1‑139 (MQ=255)
                                                                                                                                   aTGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATACACCGatctatct  <  2:24280/139‑1 (MQ=255)
                                                                                                                                    tGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGa                                                     <  2:212777/87‑1 (MQ=255)
                                                                                                                                    tGTAGAGGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGa                                                     >  1:212777/1‑87 (MQ=255)
                                                                                                                                          |                                                                                                                                   
GCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTGCCACGCCAGCCAGGCTTGCGGAATACTAATATTCAGTTGCTGATTGGCTTGATCGAAATTGAAGAGCATTTCAGGTCGGGAACTGAAATCAATACACCGATCTATCT  >  NC_000913/751594‑751862

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