Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,197,153 C→T 27.3% G117D (GGC→GAC)  mrp ← antiporter inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,197,1530CT27.3% 52.6 / 20.0 33G117D (GGC→GAC) mrpantiporter inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (15/9);  new base T (3/6);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGATTCTTCACGCGTTTCAG  >  W3110S.gb/2197096‑2197223
                                                         |                                                                      
ccTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAAc                                                            >  1:2303928/1‑70 (MQ=255)
 cTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGTCTTTACCGGAACTc                                                          <  1:541921/71‑1 (MQ=255)
 cTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCgg                                                               <  1:921762/66‑1 (MQ=255)
 cTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACt                                                           <  1:2672502/70‑1 (MQ=255)
        cgccagcgccagTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGTCTTTACCGGAACTCACCGCGa                                                   <  1:462430/71‑1 (MQ=255)
          ccagcgccagTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGAt                                                  <  1:516879/70‑1 (MQ=255)
             gcgccagTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCgg                                                               <  1:2567018/54‑1 (MQ=255)
                      ccAGGTTTACCGCCGTGGACGATTTACCCACGCCGTCTTTACCGGAACTCACCGCGAt                                                  <  1:4060897/58‑1 (MQ=255)
                         ggTTTACCGCCGTGGACGATTTACCCCCGCCGCCTTTACCGGa                                                              >  1:2098488/1‑43 (MQ=255)
                          gTTTACCGCCGTGGACGATTTACCCACGCCGTCTTTACCGGAACTCACcgc                                                     <  1:3264548/51‑1 (MQ=255)
                           tttACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACt                                                           >  1:2651742/1‑44 (MQ=255)
                           tttACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACt                                                           >  1:1692802/1‑44 (MQ=255)
                            ttACCGCCGTGGACGATTTACCCACGCCGTCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGtt                                >  1:2487028/1‑70 (MQ=255)
                            ttACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATaa                                                >  1:1308658/1‑54 (MQ=255)
                            ttACCGCCGTGGACGATTTACCCACGCCGCCGTTACCGGaa                                                             >  1:2175693/1‑41 (MQ=255)
                                 gccgTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTcc                          <  1:3754979/71‑1 (MQ=255)
                                  ccgTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACAcc                                   <  1:1651491/61‑1 (MQ=255)
                                    gTGGACGATTTACCCACGCCGTCTTTACCGGAACTCa                                                         >  1:1570606/1‑37 (MQ=255)
                                     ttgACGATTTACCCACGCCGTCTTTACCGGAACTCACCGCGATaa                                                <  1:2290855/43‑1 (MQ=38)
                                       gACGATTTACCCACGCCGCCTTTACCGGAACTCAc                                                        >  1:3122601/1‑35 (MQ=255)
                                       gACGATTTACCCACGCCGCCTTTACCGGAACTCACcgc                                                     >  1:3571085/1‑38 (MQ=255)
                                           aTTTACCCACGCCGTCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGAttc                <  1:1877578/71‑1 (MQ=255)
                                             ttACCCACGCCGTCTTTACCGGAACTCACCGCGATAATGttc                                           >  1:3876491/1‑42 (MQ=255)
                                             ttACCCACGCCGCCTTTACCGGAACTCACCGCGATaa                                                >  1:2663029/1‑37 (MQ=255)
                                             ttACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttct                                          >  1:143322/1‑43 (MQ=255)
                                             ttACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTaa                              >  1:805782/1‑55 (MQ=255)
                                              tACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGttctt                                         >  1:1455529/1‑43 (MQ=255)
                                              tACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGAttctt              <  1:2938486/70‑1 (MQ=255)
                                              tACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGAttctt              >  1:3943588/1‑70 (MQ=255)
                                              tACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGAttctt              <  1:802110/70‑1 (MQ=255)
                                              tACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGAttctt              >  1:1659222/1‑70 (MQ=255)
                                                cccACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTa                               <  1:941938/51‑1 (MQ=255)
                                                   aCGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGATTCTTCAcgcg        >  1:3539617/1‑71 (MQ=255)
                                                      ccgtcTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGATTCTTCACGCGttt     >  1:1225466/1‑71 (MQ=255)
                                                      ccgccTTTACCGGAACTCACCGCGATAATGTTCTTcaca                                     >  1:2606812/1‑39 (MQ=255)
                                                      ccgccTTTACCGGAACTCACCGCGATAATGTTCTTcaca                                     >  1:1387442/1‑39 (MQ=255)
                                                      ccgccTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTaa                              <  1:1672200/46‑1 (MQ=255)
                                                        gccTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGATTCTTCACGCGttt     >  1:2420938/1‑69 (MQ=255)
                                                         ccTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGATTCTTCACGCGTTTCag  <  1:2624236/71‑1 (MQ=255)
                                                         |                                                                      
CCTTCTGCCGCCAGCGCCAGTGCCAGGTTTACCGCCGTGGACGATTTACCCACGCCGCCTTTACCGGAACTCACCGCGATAATGTTCTTCACACCGTTAATTCCCGGCTGATTCTTCACGCGTTTCAG  >  W3110S.gb/2197096‑2197223

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