Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,392,088 G→T 20.8% L655M (CTG→ATG)  yjeP ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,392,0880GT20.8% 136.2 / 16.5 48L655M (CTG→ATG) yjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base G (14/24);  new base T (6/4);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAATCATCATGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCA  >  AP009048/4392019‑4392154
                                                                     |                                                                  
cgcCAATCATCATGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGc                                                                   <  1:730761/71‑1 (MQ=255)
    aaTCATCATGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGaa                                                               >  1:2007172/1‑71 (MQ=255)
        atcatGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATccc                                                           >  1:1964165/1‑71 (MQ=255)
        atcatGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATccc                                                           >  1:1821039/1‑71 (MQ=255)
         tcatGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCg                                                          >  1:2080798/1‑71 (MQ=255)
           atGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAg                                                                    >  1:1003185/1‑59 (MQ=255)
            tGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGcc                                                       <  1:1203841/71‑1 (MQ=255)
              ttCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCtt                                                     >  1:1146950/1‑71 (MQ=255)
              ttCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCt                                                      <  1:1239297/70‑1 (MQ=255)
                ccACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCtttt                                                   >  1:546557/1‑71 (MQ=255)
                ccACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCttt                                                    <  1:1559972/70‑1 (MQ=255)
                ccACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCttt                                                    <  1:2528720/70‑1 (MQ=255)
                 cacaGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCttttt                                                  <  1:811965/71‑1 (MQ=255)
                 cacaGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCttttt                                                  <  1:2337439/71‑1 (MQ=255)
                 cacaGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCttttt                                                  <  1:1296217/71‑1 (MQ=255)
                  acaGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATc                                                             >  1:3822284/1‑59 (MQ=255)
                       atatGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCt                                            <  1:1325152/71‑1 (MQ=255)
                          tggttggtAATGTTTTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGa                                          <  1:494221/70‑1 (MQ=255)
                          tggttggtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCAGGCCTTTTTCAGGCTGa                                          <  1:2838400/70‑1 (MQ=255)
                            gttggtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACATCGGAATCCCGGCCTTTTTCAGGCTGAggg                                       >  1:1823007/1‑71 (MQ=255)
                             ttggtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGt                                      <  1:575450/71‑1 (MQ=255)
                                gtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGac                                   <  1:617353/71‑1 (MQ=255)
                                gtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGac                                   >  1:3437179/1‑71 (MQ=255)
                                gtAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGa                                    <  1:808887/70‑1 (MQ=255)
                                 tAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTg                                           <  1:720221/62‑1 (MQ=255)
                                 tAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGac                                   <  1:2129614/70‑1 (MQ=255)
                                   aTGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGc                                             >  1:3661787/1‑58 (MQ=255)
                                      ttgtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCttttt                                                  >  1:915819/1‑50 (MQ=255)
                                        gtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCgg                                                         >  1:391308/1‑41 (MQ=255)
                                        gtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGcc                            <  1:231652/70‑1 (MQ=255)
                                        gtCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGc                             >  1:4051022/1‑69 (MQ=255)
                                              gcTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCtttt                                                   <  1:2829127/41‑1 (MQ=255)
                                              gcTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCAc                     >  1:2893986/1‑71 (MQ=255)
                                                   ccTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAggg                                       <  1:545726/48‑1 (MQ=255)
                                                     tctTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGa                                          >  1:4361409/1‑43 (MQ=255)
                                                      ctTTGTTGAGATACATCGGAATCCCGGCCTTTTTCa                                                >  1:1214659/1‑36 (MQ=255)
                                                      ctTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGc                                             >  1:792910/1‑39 (MQ=255)
                                                       tttGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAggg                                       <  1:3898920/44‑1 (MQ=255)
                                                       tttGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAAt              >  1:3927784/1‑69 (MQ=255)
                                                       tttGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAAt              <  1:1488310/69‑1 (MQ=255)
                                                         tgttgAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGa            <  1:430109/69‑1 (MQ=255)
                                                            tgAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGa         >  1:1455479/1‑69 (MQ=255)
                                                             gagaTACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCaa                <  1:2485354/61‑1 (MQ=255)
                                                                aTACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaag    <  1:586780/70‑1 (MQ=255)
                                                                 tACATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaagca  <  1:1553976/71‑1 (MQ=255)
                                                                 tACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaagca  <  1:3353625/71‑1 (MQ=255)
                                                                 tACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaagca  <  1:265438/71‑1 (MQ=255)
                                                                 tACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaagca  <  1:960873/71‑1 (MQ=255)
                                                                  aCATCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATaaagca  >  1:3826369/1‑70 (MQ=255)
                                                                     |                                                                  
CGCCAATCATCATGTTCCACAGCATATGGTTGGTAATGTTGTCGCCGCTGCCCTCTTTGTTGAGATACAGCGGAATCCCGGCCTTTTTCAGGCTGAGGGTGACCACCGCCAGCGCACCGCAAATGAGGATAAAGCA  >  AP009048/4392019‑4392154

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