Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,342 A→T 11.8% V713V (GTT→GTA ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,3420AT11.8% 93.5 / 4.8 34V713V (GTT→GTAybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base T (2/2);  total (18/16)
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.54e-01

TTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCACC  >  NC_000913/843196‑843488
                                                                                                                                                  |                                                                                                                                                  
tttATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAAc                                                                                                                                                  >  1:694649/1‑149 (MQ=255)
    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)
                                                                                                                                                  |                                                                                                                                                  
TTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCACC  >  NC_000913/843196‑843488

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