Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,347 C→G 12.9% G712R (GGC→CGC)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,3470CG12.9% 84.4 / 7.5 31G712R (GGC→CGC) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base C (14/13);  new base G (2/2);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

TTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCACCAATAAT  >  NC_000913/843200‑843494
                                                                                                                                                   |                                                                                                                                                   
ttGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgccc                                                                                                                                                    >  2:706887/1‑149 (MQ=255)
   ccTGGACCTTTATGACCGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAcgcccgcc                                                                                                                                                 <  1:626716/149‑1 (MQ=255)
     tGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGGGCGTTCGCGCGCCAg                                                                                                                                               <  1:240816/149‑1 (MQ=255)
       gACCTTTAGGAACGACGAGGAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGt                                                                                                                                             <  1:405350/149‑1 (MQ=255)
       gACCTGTAGGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCCCCTGATAAGTCTGCACTGGCGGGCGTTCGCGCGCCAgtg                                                                                                                                             <  1:4914/149‑3 (MQ=255)
            ttATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGaaa                                                                                                                                        <  1:16556/149‑1 (MQ=255)
                  acgacgAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGATCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGCTAAGTCTGCACTGGCGCGCGAAAGCCCGCCAGGTCGAAATGgttt                                                                                                                                  >  2:604474/1‑149 (MQ=255)
                    gacgaGCAGCGCGGCCGTTTTGTCAATCGCCAGCGTTAAAGCCGTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCCTTATCGTCATCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtttttt                                                                                                                                <  1:666563/149‑1 (MQ=255)
                           agcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTgg                                                                                                                         <  1:764333/149‑1 (MQ=255)
                                    gTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCCGCGGTAACGGTTCCGCCGGGGTCGCGCCCGGAGCAGGCAGCCCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAgg                                                                                                                >  1:630721/1‑149 (MQ=255)
                                             aTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGt                                                                                                       >  1:127943/1‑149 (MQ=255)
                                                    gCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTcc                                                                                                >  2:724807/1‑149 (MQ=255)
                                                          aaGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGtt                                                                                          >  1:494025/1‑149 (MQ=255)
                                                           aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGttt                                                                                         <  1:270453/149‑1 (MQ=255)
                                                                 ggTGCCCCCGGCGGTAACGGTTCAGCCGGGGTCGAGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGACGGCAAAGCTTACCGTCCACTGTTTGAGTgg                                                                                   <  1:417003/149‑1 (MQ=255)
                                                                         cggcggTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCgtgg                                                                           >  1:642150/1‑149 (MQ=255)
                                                                                  cGGTTCCGCCGGGGGCGCGCCCGGGGCAGGCAGCACCGGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaaa                                                                  <  2:324602/149‑1 (MQ=255)
                                                                                  cGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCCGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACGGTTTGAGTGGCAGCGTGGTGAACGaaa                                                                  <  1:362995/149‑1 (MQ=255)
                                                                                  cGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaaa                                                                  >  2:739272/1‑149 (MQ=255)
                                                                                           cGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGGGAACGAAACACGCAgtg                                                         >  2:307280/1‑149 (MQ=255)
                                                                                               gTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTaaa                                                     >  2:767591/1‑149 (MQ=255)
                                                                                                    gcCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCgg                                                >  1:344638/1‑149 (MQ=255)
                                                                                                         gAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCATAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGAAtca                                           >  2:595729/1‑147 (MQ=255)
                                                                                                            caggcagCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCtt                                        >  2:737978/1‑149 (MQ=255)
                                                                                                            caggcagCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCtt                                        <  1:159628/149‑1 (MQ=255)
                                                                                                                   ccccTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGa                                 <  2:301488/147‑1 (MQ=255)
                                                                                                                       tGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCCAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAATGAAACACGCAGTGTAAACGCGGTATTGCTTACGCCGGCatt                             >  1:373093/1‑147 (MQ=255)
                                                                                                                         aTAAGTCTGCACTGGCGCGCGATCGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTATCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGGAAACGCGGTATTGCTTAAGCCGACAATcc                           <  1:506901/149‑1 (MQ=255)
                                                                                                                         aTAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATcc                           >  1:801942/1‑149 (MQ=255)
                                                                                                                              ttcgaccTGGCGGGCGTTCGCGCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGAGTGGTGAACGAAACACGCAGTGTAAACGCGGGATTGCTTAAGCCGACAATCCCGGcc                      >  2:726556/7‑148 (MQ=255)
                                                                                                                                       gcgggcGTACGCGCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGACGCGAACCTGTTGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGt             >  1:583437/1‑149 (MQ=255)
                                                                                                                                           gcgAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCAc         <  2:694649/149‑1 (MQ=255)
                                                                                                                                            cgAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCAcc        <  1:767591/149‑1 (MQ=255)
                                                                                                                                                  cccgcccgCTCGAAATGTTTGTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCACCaataat  <  2:668794/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCACCAATAAT  >  NC_000913/843200‑843494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: