Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,939,281 C→T 45.5% E207K (GAA→AAA)  ygdE ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,939,2810CT45.5% 10.5 / 57.0 44E207K (GAA→AAA) ygdEpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base T (10/10);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

AGTTGGTAGGTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCGGCGC  >  W3110S.gb/2939214‑2939350
                                                                   |                                                                     
aGTTGGTAGGTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATc                                                                    <  1:3425928/71‑1 (MQ=255)
 gTTGGTAGGTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATc                                                                    <  1:787549/70‑1 (MQ=255)
          tCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATcc                                                          >  1:2987274/1‑71 (MQ=255)
           ccAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATcc                                                          >  1:2774758/1‑70 (MQ=255)
           ccAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCg                                                         >  1:4014226/1‑71 (MQ=255)
            cAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCg                                                         >  1:589467/1‑70 (MQ=255)
               ccgccAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAgg                                                     <  1:3795491/71‑1 (MQ=255)
                    aGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATa                                                 >  1:1726376/1‑70 (MQ=255)
                        cAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGTAATAAacac                                            <  1:3570360/71‑1 (MQ=255)
                        cAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAacac                                            <  1:1230934/71‑1 (MQ=255)
                        cAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAacac                                            <  1:2229924/71‑1 (MQ=255)
                            cgcCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGa                                        <  1:240299/71‑1 (MQ=255)
                              ccTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGa                                        <  1:2012952/69‑1 (MQ=255)
                              ccTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGa                                        <  1:2812356/69‑1 (MQ=255)
                              ccTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGa                                        <  1:3565955/69‑1 (MQ=255)
                                 aaaTCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAAcaca                                           >  1:3458733/1‑63 (MQ=255)
                                 aaaTCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGaa                                       >  1:185941/1‑67 (MQ=255)
                                 aaaTCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGaa                                       >  1:398554/1‑67 (MQ=255)
                                 aaaTCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGaa                                       >  1:3016540/1‑67 (MQ=255)
                                  aaTCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGaaa                                      >  1:3490888/1‑67 (MQ=255)
                                  aaTCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGaaa                                      >  1:1582729/1‑67 (MQ=255)
                                  aaTCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGaaa                                      >  1:1578086/1‑67 (MQ=255)
                                  aaTCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGaaa                                      >  1:1417069/1‑67 (MQ=255)
                                  aaTCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGaaa                                      >  1:312301/1‑67 (MQ=255)
                                        ccGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTcc                            <  1:831964/71‑1 (MQ=255)
                                         cGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCa                           <  1:3452182/71‑1 (MQ=255)
                                              acaTCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATa                                                 >  1:1176573/1‑44 (MQ=255)
                                              acaTCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATa                                                 >  1:3967539/1‑44 (MQ=255)
                                               caTCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAAt                                     <  1:1654624/55‑1 (MQ=255)
                                               caTCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAAt                                     <  1:826749/55‑1 (MQ=255)
                                               caTCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAAt                                     <  1:1333635/55‑1 (MQ=255)
                                               caTCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAAt                                     <  1:1446326/55‑1 (MQ=255)
                                               caTCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGaa                                                   >  1:2179020/1‑41 (MQ=255)
                                               caTCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGaa                                                   >  1:1904907/1‑41 (MQ=255)
                                               caTCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGaa                                                   >  1:118830/1‑41 (MQ=255)
                                                 tCCCGTTCGCCAGGTGTTCATCCCACTCATCCGCAGTAATAAACACATGAAATGCCTCTTCCAGTTTg                      <  1:2523123/68‑1 (MQ=255)
                                                 tCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTcc                            <  1:529397/62‑1 (MQ=255)
                                                   ccGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGaaa                                      >  1:2544650/1‑50 (MQ=255)
                                                     gTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTg                >  1:3433096/1‑70 (MQ=255)
                                                     gTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTgg               <  1:3989544/71‑1 (MQ=255)
                                                     gTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTgg               <  1:66622/71‑1 (MQ=255)
                                                      ttCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGa              <  1:3115713/71‑1 (MQ=255)
                                                        cGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAAc            >  1:2491781/1‑71 (MQ=255)
                                                            aGGCGTTTATCCCACTCATCCGCAGGAATAAacac                                            <  1:3348529/35‑1 (MQ=255)
                                                                   cATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCGgcgc  <  1:714788/70‑1 (MQ=255)
                                                                   |                                                                     
AGTTGGTAGGTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCGGCGC  >  W3110S.gb/2939214‑2939350

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: