Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,491,067 G→A 100% V56V (GTG→GTA yjgP → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,491,0670GA96.2% 78.1 / ‑3.6 26V56V (GTG→GTAyjgPconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (6/19);  total (6/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAC  >  W3110S.gb/4491013‑4491123
                                                      |                                                        
ggTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGc                                          <  1:1875134/71‑1 (MQ=255)
 gTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCt                                         <  1:357735/71‑1 (MQ=255)
  ttGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCt                                         <  1:1602155/70‑1 (MQ=255)
              ttCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGa                                                     <  1:669061/46‑1 (MQ=255)
              ttCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGa                                                     <  1:1505797/46‑1 (MQ=255)
              ttCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGa                                                     <  1:1686319/46‑1 (MQ=255)
               tCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGa                                                     <  1:1838185/45‑1 (MQ=255)
                  ggCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTACCGGa                                                     <  1:136926/42‑1 (MQ=255)
                                tctcCCTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGc                                          <  1:1027448/39‑1 (MQ=255)
                                    ccTTCTCGGGTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGgctgct        <  1:1198209/69‑1 (MQ=255)
                                    ccTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGgctgct        <  1:1758278/69‑1 (MQ=255)
                                    ccTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGgctgct        <  1:1849142/69‑1 (MQ=255)
                                    ccTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGgctgct        <  1:1443479/69‑1 (MQ=255)
                                     cTTCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCtgat     <  1:970241/71‑1 (MQ=255)
                                      ttCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTa                                       <  1:1964865/36‑1 (MQ=255)
                                      ttCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTa                                       <  1:347714/36‑1 (MQ=255)
                                      ttCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTa                                       <  1:1328866/36‑1 (MQ=255)
                                      ttCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTa                                       <  1:390396/36‑1 (MQ=255)
                                      ttCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTa                                       <  1:533369/36‑1 (MQ=255)
                                      ttCTCGGGTTGGGCGTACCGGAAATGGCGCAGCTTa                                       <  1:646913/36‑1 (MQ=255)
                                        ctcGGGTTGGGCGTACCGGAAATGGCGCAGCTTATc                                     >  1:966870/1‑36 (MQ=255)
                                        ctcGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCtgatga   >  1:451101/1‑70 (MQ=255)
                                        ctcGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc  >  1:1614856/1‑71 (MQ=255)
                                        ctcGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc  >  1:1199922/1‑71 (MQ=255)
                                        ctcGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc  >  1:733547/1‑71 (MQ=255)
                                        ctcGGGTTGGGCGTACCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAc  >  1:758758/1‑71 (MQ=255)
                                                      |                                                        
GGTTGACGGCGATATTCCGGCGAATCTGGTGCTCTCCCTTCTCGGGTTGGGCGTGCCGGAAATGGCGCAGCTTATCCTGCCATTAAGCCTGTTCCTCGGGCTGCTGATGAC  >  W3110S.gb/4491013‑4491123

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