Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,012,950 A→G 37.8% V255A (GTG→GCG)  yqeC ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,012,9500AG37.8% 25.7 / 38.1 37V255A (GTG→GCG) yqeCconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/10);  new base G (6/8);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTT  >  W3110S.gb/3012883‑3013019
                                                                   |                                                                     
tgAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTcgca                                                                     >  1:549927/1‑70 (MQ=255)
  aaaaTATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAAt                                                                  >  1:1476598/1‑71 (MQ=255)
  aaaaTATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAAt                                                                  >  1:5134036/1‑71 (MQ=255)
  aaaaTATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAAt                                                                  >  1:3684434/1‑71 (MQ=255)
           aTACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTgcgcgc                                                          <  1:1694704/70‑1 (MQ=255)
            tACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCGCaa                                                                    <  1:3639635/59‑1 (MQ=255)
                    tCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCGCAAATCTGCGCGCGAtt                                                      <  1:4324999/65‑1 (MQ=255)
                      aaTGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCGCAAATCTGCGCGCGATTGCAGGa                                                <  1:3493989/69‑1 (MQ=255)
                      aaTGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCGCAAATCTGCGCGCGATTGCAGGATg                                              <  1:3765859/71‑1 (MQ=255)
                          gATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCGCaa                                                                    >  1:1071950/1‑45 (MQ=255)
                            tatCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGTATGTTCTTg                                        <  1:3005187/71‑1 (MQ=255)
                                cTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCGCAAATCTGCGCGCGATTGCAGGATGTTCTTGtata                                    >  1:2512609/1‑71 (MQ=255)
                                 tttCAGTAACCCGGAATACCCGGGCCGCTAATTCGCAAATCTGCGCGCGATTGCAGGATGTTCTTGTatat                                   >  1:2469412/1‑71 (MQ=255)
                                 tttCAGTAACCCGGAATACCCGGGCCGCTAATTCACaa                                                                    <  1:4648050/38‑1 (MQ=255)
                                 tttCAGTAACCCGGAATACCCGGGCGCTAATTCACaa                                                                    <  1:4207870/37‑1 (MQ=38)
                                     aGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGa                                                        >  1:2382749/1‑46 (MQ=255)
                                        aaCCCGGAATACCCGGGCCGCTAATTCACAAATCTgcgcgc                                                          <  1:1260501/41‑1 (MQ=255)
                                        aaCCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccag                            >  1:2307608/1‑71 (MQ=255)
                                         aCCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                           <  1:3922342/71‑1 (MQ=255)
                                          cccGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGa                                                        >  1:3341587/1‑41 (MQ=255)
                                          cccGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGa                                                        >  1:412102/1‑41 (MQ=255)
                                           ccGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagcca                         <  1:2908215/71‑1 (MQ=255)
                                            cGGAATACCCGGGCCGCTAATTCGCAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagccag                        <  1:2736118/71‑1 (MQ=255)
                                             ggAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagc                           >  1:4495733/1‑67 (MQ=255)
                                                 tACCCGGGCCGCTAATTCGCAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGc                   >  1:4120450/1‑71 (MQ=255)
                                                 tACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAtt                     <  1:4838671/69‑1 (MQ=255)
                                                            cTAATTCACAAATCTGCGCGCGATTGCAGGATGtt                                            <  1:577830/35‑1 (MQ=255)
                                                            cTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGtat                                     <  1:2157056/42‑1 (MQ=255)
                                                            cTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGcc                  >  1:2254084/1‑61 (MQ=255)
                                                             tAATTCGCAAATCTGCGCGCGATTGCAGGATGTTCTTg                                        <  1:4211348/38‑1 (MQ=255)
                                                              aaTTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGc                   >  1:240629/1‑58 (MQ=255)
                                                               aTTCGCAAATCTGCGCGCGATTGCAGGATGTTCTTGta                                      >  1:3863061/1‑38 (MQ=255)
                                                               aTTCGCAAATCTGCGCGCGATTGCAGGATGTTCTTGta                                      >  1:2803741/1‑38 (MQ=255)
                                                                ttCGCAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAtt                     <  1:2224639/54‑1 (MQ=255)
                                                                ttCGCAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGAtt                     <  1:1174630/54‑1 (MQ=255)
                                                                 tCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCccagccag                        >  1:4816494/1‑50 (MQ=255)
                                                                 tCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGa                       >  1:2072326/1‑51 (MQ=255)
                                                                 tCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGtgttgt   >  1:2686722/1‑71 (MQ=255)
                                                                  cacaAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCATCCAGATTGCctc                <  1:3713822/57‑1 (MQ=255)
                                                                  cacaAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTgttgtt  >  1:3844102/1‑71 (MQ=255)
                                                                   |                                                                     
TGAAAATATTCATACCAACCTCAATGGATATCCTTTCAGTAACCCGGAATACCCGGGCCGCTAATTCACAAATCTGCGCGCGATTGCAGGATGTTCTTGTATATCGCCCAGCCAGATTGCCTCTACGTCGTGTTGTT  >  W3110S.gb/3012883‑3013019

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: