Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,081,459 A→G 10.8% F36L (TTT→CTT)  yeeX ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,081,4590AG10.8% 115.2 / 6.5 46F36L (TTT→CTT) yeeXconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base G (2/3);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAATGTGAGAGTCG  >  W3110S.gb/2081393‑2081526
                                                                  |                                                                   
ggATCTTTTTCTCAAGGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTcc                                                                 <  1:946575/71‑1 (MQ=255)
ggATCTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTcc                                                                 <  1:1552593/71‑1 (MQ=255)
               cGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCt                                                       >  1:113365/1‑66 (MQ=255)
                 tCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATg                                                 <  1:1365587/70‑1 (MQ=255)
                   cTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACaaa                                                                     >  1:1218387/1‑48 (MQ=255)
                   cTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAg                                              >  1:390899/1‑71 (MQ=255)
                   cTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAg                                              >  1:1371699/1‑71 (MQ=255)
                   cTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGaa                                                   >  1:1101775/1‑66 (MQ=255)
                   cTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGaa                                                   >  1:1125491/1‑66 (MQ=255)
                     gAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTa                                                             >  1:455027/1‑54 (MQ=255)
                     gAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTa                                                             >  1:614500/1‑54 (MQ=255)
                     gAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGc                                            >  1:2753309/1‑71 (MQ=255)
                     gAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGc                                            >  1:1839158/1‑71 (MQ=255)
                          tCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCa                                                                >  1:1483735/1‑46 (MQ=255)
                          tCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGaa                                               >  1:2516577/1‑63 (MQ=255)
                          tCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTgg                                        >  1:1362338/1‑70 (MQ=255)
                          tCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTgg                                        >  1:2064649/1‑70 (MQ=255)
                          tCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTgg                                        <  1:2130154/70‑1 (MQ=255)
                              gTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGtt                                   <  1:1721911/71‑1 (MQ=255)
                              gTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGtt                                   <  1:2280990/71‑1 (MQ=255)
                                  cAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGt                                    >  1:1060782/1‑66 (MQ=255)
                                  cAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGt                                    >  1:961926/1‑66 (MQ=255)
                                  cAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTcc                                <  1:375110/70‑1 (MQ=255)
                                     tttGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAt                              <  1:2536343/69‑1 (MQ=255)
                                       tGTTCTTTCGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           <  1:2597373/70‑1 (MQ=255)
                                       tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGc                                            >  1:1961768/1‑53 (MQ=255)
                                       tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGc                                            >  1:2534065/1‑53 (MQ=255)
                                       tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           >  1:2218438/1‑70 (MQ=255)
                                       tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           >  1:1290238/1‑70 (MQ=255)
                                       tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           <  1:945015/70‑1 (MQ=255)
                                       tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           >  1:354538/1‑70 (MQ=255)
                                       tGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           <  1:2143089/70‑1 (MQ=255)
                                         ttCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTtat                        <  1:475842/71‑1 (MQ=255)
                                         ttCTTACGACGGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                         <  1:2372379/70‑1 (MQ=255)
                                         ttCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                         <  1:2749419/70‑1 (MQ=255)
                                         ttCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                         <  1:1630043/70‑1 (MQ=255)
                                         ttCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                         <  1:714240/70‑1 (MQ=255)
                                              acgacgGAACAGACGAACAAGTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATAcc                     <  1:164557/69‑1 (MQ=255)
                                              acgacgGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           <  1:859970/63‑1 (MQ=255)
                                               cgacgGAACAGACGAACAAATTCCAGTACGTCCTgg                                                     <  1:114893/36‑1 (MQ=255)
                                                  cgGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGAt               <  1:2762404/71‑1 (MQ=255)
                                                         gACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTg                                         <  1:1349083/38‑1 (MQ=255)
                                                         gACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           <  1:743338/52‑1 (MQ=255)
                                                           cGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTta                         >  1:1143609/1‑52 (MQ=255)
                                                           cGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAAtgtta        >  1:2310826/1‑67 (MQ=255)
                                                             aaCAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAATGTGagag     >  1:576937/1‑70 (MQ=255)
                                                               cAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCAATTTTTATACCCCTGATAATGTGAGAGTCg  >  1:583477/1‑71 (MQ=255)
                                                                aaaTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           >  1:563537/1‑45 (MQ=255)
                                                                aaaTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           >  1:560411/1‑45 (MQ=255)
                                                                aaaTTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCAtttt                           >  1:1488328/1‑45 (MQ=255)
                                                                  |                                                                   
GGATCTTTTTCTCAACGTCCTGAATTTCACGTTGCAGTTTGTTCTTACGACGGAACAGACGAACAAATTCCAGTACGTCCTGGAATGAAGGCTTGGTAGTTTCCATTTTTATACCCCTGATAATGTGAGAGTCG  >  W3110S.gb/2081393‑2081526

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: