Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,428,249 A→G 14.6% I246T (ATT→ACT)  yfcI ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,428,2490AG14.6% 120.7 / 13.4 48I246T (ATT→ACT) yfcIconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (4/3);  total (25/23)
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.93e-01

CTGAGCAATACGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAAC  >  W3110S.gb/2428182‑2428315
                                                                   |                                                                  
cTGAGCAATACGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATg                                                                  <  1:450128/70‑1 (MQ=255)
          cGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATg                                                                  >  1:516466/1‑60 (MQ=255)
          cGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATg                                                                  >  1:3689518/1‑60 (MQ=255)
          cGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATg                                                                  >  1:2461053/1‑60 (MQ=255)
          cGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGtt                                                        >  1:2324429/1‑70 (MQ=255)
                       cGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTATCTtgt                                          >  1:966000/1‑71 (MQ=255)
                        gtgtTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTGATCAGTTTCTCCTTTTCTtgt                                          <  1:2581332/70‑1 (MQ=255)
                        gtgtTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTtgt                                          <  1:1120117/70‑1 (MQ=255)
                        gtgtTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTtgt                                          <  1:398188/70‑1 (MQ=255)
                        gtgtTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTtgt                                          <  1:1822154/70‑1 (MQ=255)
                            ttGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCttt                                                <  1:1441029/60‑1 (MQ=255)
                             tGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTgcg                                    >  1:2801559/1‑71 (MQ=255)
                             tGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTgcg                                    >  1:3228169/1‑71 (MQ=255)
                             tGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTgcg                                    >  1:2644981/1‑71 (MQ=255)
                                   gCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTc                               >  1:1068889/1‑70 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAGTGGTCATCAGTTTCTCCTTTTCTtgt                                          <  1:3048809/56‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAGTGGTCATCAGTTTCTCCTTTTCTtgt                                          <  1:2812174/56‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           <  1:3078530/71‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           <  1:1839738/71‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           <  1:2613593/71‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           <  1:2716806/71‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           <  1:3617649/71‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           <  1:2752325/71‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           <  1:3485878/71‑1 (MQ=255)
                                      ttGCGCCTTCTTCACGTAATCAGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           <  1:1407490/71‑1 (MQ=255)
                                        gcgcCTTCTTCACGTAATCTGTGAGCAGTGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTAt                         >  1:3408021/1‑71 (MQ=255)
                                        gcgcCTTCTTCACGTAATCTGTCAGCAGTGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           >  1:3386075/1‑69 (MQ=255)
                                        gcgcCTTCTTCACGTAATCTGTCAGCAGTGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTa                          >  1:1838893/1‑70 (MQ=255)
                                        gcgcCTTCTTCACGTAATCTGTCAGCAGTGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTAt                         >  1:2691601/1‑71 (MQ=255)
                                        gcgcCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTa                          >  1:1364649/1‑70 (MQ=255)
                                        gcgcCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTAt                         >  1:387757/1‑71 (MQ=255)
                                         cgcCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTAt                         <  1:1329888/70‑1 (MQ=255)
                                         cgcCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTAt                         <  1:3418682/70‑1 (MQ=255)
                                           ccTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCttt                                                <  1:684251/45‑1 (MQ=255)
                                             ttcttcACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCtt                                                 <  1:2777001/42‑1 (MQ=255)
                                             ttcttcACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCtt                                                 <  1:639237/42‑1 (MQ=255)
                                             ttcttcACGTAATCTGTCAGCAATGGTCATCAGTCTCTCCtt                                                 <  1:3331527/42‑1 (MQ=255)
                                                             tCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGaa     >  1:2579971/1‑70 (MQ=255)
                                                               aGCAGTGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAAc  <  1:393962/71‑1 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCt                           >  1:3252265/1‑46 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCAc                    >  1:1827874/1‑53 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGaaaa   >  1:1458021/1‑70 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGaaaa   >  1:2732595/1‑70 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGaaaa   >  1:2080083/1‑70 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGaaaa   >  1:998379/1‑70 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAAc  >  1:755971/1‑71 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAAc  >  1:1319650/1‑71 (MQ=255)
                                                               aGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAAc  >  1:977562/1‑71 (MQ=255)
                                                                   |                                                                  
CTGAGCAATACGCAGGGCTTCTTCGTGTTTGCCCTGCATTGCGCCTTCTTCACGTAATCTGTCAGCAATGGTCATCAGTTTCTCCTTTTCTTGTGGTGCGCGTTCCGCTATCTCACCAATAAATGCACGAAAAC  >  W3110S.gb/2428182‑2428315

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: