Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,198,610 A→G 58.3% F338L (TTC→CTC)  ygiF ← predicted adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,198,6100AG58.3% 18.6 / 40.8 36F338L (TTC→CTC) ygiFpredicted adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (6/9);  new base G (6/15);  total (12/24)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCATGCTTTGCTTACCAAC  >  W3110S.gb/3198545‑3198674
                                                                 |                                                                
cTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGa                                                               <  1:91570/69‑1 (MQ=255)
      gAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCt                                                         <  1:111858/69‑1 (MQ=255)
      gAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCt                                                         <  1:24161/69‑1 (MQ=255)
      gAAAACGCTTTTCAGTTCAGCGGCATGACGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCt                                                         <  1:1180225/69‑1 (MQ=255)
      gAAAACGCTTTTCAGTTCAGCGGCATGACGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCt                                                         <  1:3670665/69‑1 (MQ=255)
      gAAAACGCTTTTCAGTTCAGCGGCATGACGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCt                                                         <  1:2226602/69‑1 (MQ=255)
        aaaCGCTTTTCAGTTCAGCGGCATAGCGGGAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGa                                                       >  1:1641904/1‑69 (MQ=255)
        aaaCGCTTTTCAGTTCAGAGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGa                                                       >  1:2569044/1‑69 (MQ=255)
         aaCGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGAt                                                      <  1:1976260/69‑1 (MQ=255)
           cGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGAttt                                                    <  1:1856962/69‑1 (MQ=255)
                     ttcagCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTGGCTGATTTTGCCCTgcgc                                          <  1:2902672/69‑1 (MQ=255)
                     ttcagCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTgcgc                                          <  1:1512472/69‑1 (MQ=255)
                     ttcagCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTgcgc                                          <  1:1126021/69‑1 (MQ=255)
                                ggCGGGAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTa                               <  1:1004344/69‑1 (MQ=255)
                                 gCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTa                               >  1:1969140/1‑68 (MQ=255)
                                   gggAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaaa                             <  1:3561228/68‑1 (MQ=255)
                                     gAAAGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAAt                          <  1:2860428/69‑1 (MQ=255)
                                      aaaGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTg                                      >  1:668525/1‑56 (MQ=255)
                                      aaaGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATg                         >  1:121888/1‑69 (MQ=255)
                                      aaaGATGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATg                         >  1:1363262/1‑69 (MQ=255)
                                      aaaGAAGGATATCGGCAAAGCGGTTGAGGGAGTCGCTGAatt                                                    >  1:2045447/1‑42 (MQ=25)
                                        aGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTg                                      <  1:1224372/54‑1 (MQ=255)
                                           tGGATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGc                    <  1:2708831/69‑1 (MQ=255)
                                           tGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGc                    >  1:3065036/1‑69 (MQ=255)
                                           tGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGc                    >  1:495294/1‑69 (MQ=255)
                                             gATATCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCa                  <  1:2000529/69‑1 (MQ=255)
                                               tatCGGCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCAtg                <  1:3327285/69‑1 (MQ=255)
                                                atCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCAtgc               >  1:764497/1‑69 (MQ=255)
                                                   ggCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaaa                             <  1:3036775/52‑1 (MQ=255)
                                                   ggCAAAGCGTTTGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaaa                             <  1:864978/52‑1 (MQ=255)
                                                        aGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATATAAAAATGGCTGcc                   >  1:33497/1‑57 (MQ=255)
                                                        aGCGTTTGAAGGAGTCGCTGATTTTGACCTGCGCTTTGGCATc                                 >  1:1745737/1‑43 (MQ=255)
                                                             ttGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCATGCTTTGCTTaccaac  <  1:1218253/69‑1 (MQ=255)
                                                             ttGAGGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCATGCTTTGCTTaccaac  <  1:2133514/69‑1 (MQ=255)
                                                              tGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaaa                             <  1:2655147/41‑1 (MQ=255)
                                                              tGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTaaa                             <  1:99366/41‑1 (MQ=255)
                                                                 |                                                                
CTGGCAGAAAACGCTTTTCAGTTCAGCGGCATGGCGGGAAAGATGGATATCGGCAAAGCGTTTGAAGGAGTCGCTGATTTTGCCCTGCGCTTTGGCATCTAAAAATGGCTGCCATGCTTTGCTTACCAAC  >  W3110S.gb/3198545‑3198674

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: