Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,246,753 A→G 11.5% G78G (GGT→GGC)  folE ← GTP cyclohydrolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,246,7530AG11.5% 162.2 / 12.1 61G78G (GGT→GGC) folEGTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base G (2/5);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGC  >  W3110S.gb/2246684‑2246821
                                                                     |                                                                    
ggTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGAc                                                                     >  1:132218/1‑71 (MQ=255)
     ccATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGgag                                                                <  1:2134341/71‑1 (MQ=255)
      cATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGgaga                                                               <  1:979779/71‑1 (MQ=255)
          tCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGaaaa                                                            <  1:3302885/70‑1 (MQ=255)
          tCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAAt                                                           >  1:1141654/1‑71 (MQ=255)
            aTCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAAtt                                                          <  1:1811959/70‑1 (MQ=255)
             tCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTc                                                        <  1:1327936/71‑1 (MQ=255)
                   ttCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGAcc                                                                    <  1:3100312/53‑1 (MQ=255)
                    tCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGaca                                                 >  1:4993624/1‑71 (MQ=255)
                    tCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGaca                                                 >  1:2079499/1‑71 (MQ=255)
                     cATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGa                                                   <  1:5119449/68‑1 (MQ=255)
                      aTTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACAGGAGaaaa                                                            >  1:245831/1‑58 (MQ=255)
                      aTTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACAGGAGaaaa                                                            >  1:251810/1‑58 (MQ=255)
                       ttttGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGAGCGGAGAAAATTTCa                                                       <  1:1174653/62‑1 (MQ=255)
                       ttttGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAAtctc                                                        <  1:4243412/61‑1 (MQ=255)
                       ttttGTTTTCAATGAGGGTGATTATCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGa                                                   <  1:3548929/66‑1 (MQ=255)
                         ttGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGGCCGGAGAAAATTTCATCGacatacat                                            >  1:3544617/1‑71 (MQ=255)
                         ttGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGGCCGGAGAAAATTTCATCGacatacat                                            >  1:2438093/1‑71 (MQ=255)
                          tGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACAtt                                           >  1:4480712/1‑71 (MQ=255)
                          tGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACAtt                                           >  1:3908345/1‑71 (MQ=255)
                           gTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACAtt                                           <  1:855610/70‑1 (MQ=255)
                            ttttCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACAttt                                          <  1:1513864/70‑1 (MQ=255)
                              ttCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTc                                                        <  1:1293135/54‑1 (MQ=255)
                                cAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGGCCGGAGAAAATTTCATCGACATACATTTTAgc                                      <  1:4498833/70‑1 (MQ=255)
                                cAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCAt                                                      >  1:772057/1‑54 (MQ=255)
                                cAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgc                                      <  1:3704526/70‑1 (MQ=255)
                                  aTGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgcgat                                   <  1:4313291/71‑1 (MQ=255)
                                   tGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTc                                                        <  1:2341486/49‑1 (MQ=255)
                                     aGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgcgatgc                                 >  1:3744814/1‑70 (MQ=255)
                                     aGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgcgatgc                                 <  1:1319141/70‑1 (MQ=255)
                                     aGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgcgatgc                                 >  1:2022031/1‑70 (MQ=255)
                                        gTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgc                                      <  1:2938636/62‑1 (MQ=255)
                                        gTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGcgatgcgatg                             <  1:2643829/71‑1 (MQ=255)
                                         tGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGa                                                   >  1:5042852/1‑48 (MQ=255)
                                          gATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGaca                                                 >  1:1663420/1‑49 (MQ=255)
                                           aTTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgc                                      <  1:1964677/59‑1 (MQ=255)
                                           aTTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATgcgg                          >  1:3560689/1‑71 (MQ=255)
                                            ttttCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATgcggc                         >  1:1163259/1‑71 (MQ=255)
                                             tttCGGGAAATTGGCGTAATCCAGACCGGAGAAAAttt                                                         >  1:1922521/1‑38 (MQ=255)
                                             tttCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgcga                                    >  1:2651114/1‑59 (MQ=255)
                                             ttgcGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATgcggcg                        >  1:2198050/4‑71 (MQ=255)
                                                 gggAAATTGGCGTAATCCAGGCCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCTGCGtt                      <  1:1942588/69‑1 (MQ=255)
                                                 gggAAATTGGCGTAATCCAGGCCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGtt                      <  1:4459156/69‑1 (MQ=255)
                                                 gggAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTc                    <  1:4604586/71‑1 (MQ=255)
                                                 gggAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTc                    <  1:80004/71‑1 (MQ=255)
                                                        tGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCaaa             >  1:4230604/1‑71 (MQ=255)
                                                        tGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCaa              >  1:108976/1‑70 (MQ=255)
                                                            gTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTcc                   >  1:803139/1‑61 (MQ=255)
                                                            gTAATCCAGACAGGAGAAAATTTCATCGACATACATTTTAGCGatgcgatgcgt                          >  1:4584770/1‑53 (MQ=255)
                                                             tAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAgcgat                                   <  1:2913883/44‑1 (MQ=255)
                                                               aTCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcat      <  1:981671/71‑1 (MQ=255)
                                                                tCCAGACCGGAGAAAATTTCGTCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCa               <  1:2380666/61‑1 (MQ=255)
                                                                tCCAGACCGGAGAAAATTTCATCGACATACATTTTAgc                                      <  1:3224178/38‑1 (MQ=255)
                                                                tCCAGACCGGAGAAAATTTCATCGACATACATTTTAgc                                      <  1:1014452/38‑1 (MQ=255)
                                                                 ccAGACCGGAGAAAATTTCATCGACATACATTTTAgc                                      >  1:317581/1‑37 (MQ=255)
                                                                 ccAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATG‑GGCGTTTCCATCAAACTGTCATCAg   >  1:3761179/1‑71 (MQ=255)
                                                                  cAGGCCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcatca    <  1:3970640/70‑1 (MQ=255)
                                                                  cAGGCCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGtcatca    <  1:197304/70‑1 (MQ=255)
                                                                  cAGACCGGAGAAAATTTCATCGACATACATTTTAgcgatgc                                 >  1:1370780/1‑41 (MQ=255)
                                                                  cAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCaaa             >  1:4287281/1‑61 (MQ=255)
                                                                   aGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGc  <  1:1422414/71‑1 (MQ=255)
                                                                     |                                                                    
GGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAATTTCATCGACATACATTTTAGCGATGCGATGCGGCGTTTCCATCAAACTGTCATCAGC  >  W3110S.gb/2246684‑2246821

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: ‑