Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 515,934 C→T 82.0% A43V (GCG→GTG)  ybbM → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb515,9340CT82.0% 111.2 / 18.6 50A43V (GCG→GTG) ybbMpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (27/14);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCT  >  W3110S.gb/515867‑516000
                                                                   |                                                                  
gCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGAt                                                                 <  1:851049/71‑1 (MQ=255)
     cTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAAtt                                                             >  1:2818992/1‑70 (MQ=255)
     cTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAAtt                                                             >  1:1276036/1‑70 (MQ=255)
     cTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAAtt                                                             >  1:2035987/1‑70 (MQ=255)
     cTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAAt                                                              >  1:3552070/1‑69 (MQ=255)
      ttAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAAtt                                                             >  1:1996516/1‑69 (MQ=255)
      ttAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAAtt                                                             >  1:2568521/1‑69 (MQ=255)
             gCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGatt                                                    <  1:2011793/71‑1 (MQ=255)
             gCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGatt                                                    <  1:2777215/71‑1 (MQ=255)
             gCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGatt                                                    <  1:2246246/71‑1 (MQ=255)
             gCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGatt                                                    <  1:1054657/71‑1 (MQ=255)
                   aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                              >  1:956864/1‑71 (MQ=255)
                   aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTc                                              >  1:2009641/1‑71 (MQ=255)
                   aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTc                                              >  1:2018235/1‑71 (MQ=255)
                   aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTc                                              >  1:594624/1‑71 (MQ=255)
                   aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTc                                              >  1:3427165/1‑71 (MQ=255)
                     gAAAAACTGGCGCTGGAGAAAGATATTCTCTTGAGCGTCGGGCGAGTGATAATTCAGCTGattatt                                                 <  1:1675236/66‑1 (MQ=255)
                     gAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGattatt                                                 <  1:194625/66‑1 (MQ=255)
                     gAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGattatt                                                 <  1:1773226/66‑1 (MQ=255)
                     gAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGattatt                                                 <  1:544190/66‑1 (MQ=255)
                      aaaaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGattatt                                                 <  1:2610911/65‑1 (MQ=255)
                        aaaCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTa                                         >  1:1357559/1‑71 (MQ=255)
                             ggCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGat                                                     >  1:1005638/1‑54 (MQ=255)
                                 cTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGa                                                      <  1:929919/49‑1 (MQ=255)
                                   ggAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGatt                                                    >  1:2943779/1‑49 (MQ=255)
                                   ggAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGat                                                     >  1:2831390/1‑48 (MQ=255)
                                   ggAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAg                               >  1:299073/1‑70 (MQ=255)
                                   ggAGAAAGATATTCTCTGGAGCGTCGCGCGAGTGATAATTCAGCTGatt                                                    >  1:53609/1‑49 (MQ=255)
                                     agaAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGta                             >  1:177839/1‑70 (MQ=255)
                                     agaAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGta                             >  1:314989/1‑70 (MQ=255)
                                      gaAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGtata                           <  1:2603071/71‑1 (MQ=255)
                                        aaGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAtt                         <  1:2658815/71‑1 (MQ=255)
                                        aaGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAtt                         <  1:860644/71‑1 (MQ=255)
                                        aaGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAtt                         <  1:1936979/71‑1 (MQ=255)
                                         aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTAt                                        >  1:1064632/1‑55 (MQ=255)
                                         aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                        >  1:995672/1‑71 (MQ=255)
                                         aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAttt                        >  1:3443216/1‑71 (MQ=255)
                                         aGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAtt                         >  1:2541516/1‑70 (MQ=255)
                                          gATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAtttt                       <  1:145323/71‑1 (MQ=255)
                                              ttCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGaa                                >  1:758598/1‑58 (MQ=255)
                                              ttCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGaa                                >  1:3144933/1‑58 (MQ=255)
                                              ttCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGc                   >  1:1780015/1‑71 (MQ=255)
                                              ttCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTAtg                                       >  1:3498935/1‑51 (MQ=255)
                                               tctcTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCg                  >  1:1345194/1‑71 (MQ=255)
                                                ctctGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTc                      >  1:2551128/1‑66 (MQ=255)
                                                   tGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg                <  1:868707/69‑1 (MQ=255)
                                                   tGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg                <  1:1101706/69‑1 (MQ=255)
                                                     gAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat             <  1:2935727/70‑1 (MQ=255)
                                                      aGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgatga           >  1:1474517/1‑71 (MQ=255)
                                                        cGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGatgatg         >  1:681802/1‑71 (MQ=255)
                                                               cgagcgaTAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCt  >  1:1391382/1‑71 (MQ=255)
                                                                   |                                                                  
GCAATCTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCCAGCCT  >  W3110S.gb/515867‑516000

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: