Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,092,365 A→G K13E (AAA→GAA)  hisG → ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,092,3650AG97.8% 157.4 / ‑4.0 46K13E (AAA→GAA) hisGATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (31/14);  total (32/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTT  >  W3110S.gb/2092324‑2092430
                                         |                                                                 
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaaa                                                                 >  1:3007409/1‑44 (MQ=255)
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaaa                                                                 >  1:1942325/1‑44 (MQ=255)
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGaaa                                                                 >  1:172754/1‑44 (MQ=255)
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTa                                                     >  1:2709503/1‑56 (MQ=255)
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGtgatg                                              >  1:1342463/1‑63 (MQ=255)
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCAcgc                                      >  1:1302555/1‑71 (MQ=255)
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCAc                                        >  1:1020439/1‑69 (MQ=255)
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcgc                                      >  1:1325620/1‑71 (MQ=255)
aCAAAATGACAGACAACACTCGTTTACGCATAGCTATCCAGGAATCCGGCCg                                                         >  1:2703716/1‑52 (MQ=255)
     aTGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAg                                                   >  1:2364942/1‑53 (MQ=255)
     aTGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGAc                                            >  1:2166439/1‑60 (MQ=255)
     aTGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt                                 >  1:581552/1‑71 (MQ=255)
     aTGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAAtt                                 >  1:2447026/1‑71 (MQ=255)
     aTGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAAt                                  >  1:2075026/1‑70 (MQ=255)
     aTGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAAt                                  >  1:3230494/1‑70 (MQ=255)
     aTGACAGACAACACTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAAt                                  >  1:1114469/1‑70 (MQ=255)
             caacaCTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGgcgc                         >  1:1116361/1‑71 (MQ=255)
                 acTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCa                                         >  1:3969622/1‑51 (MQ=255)
                 acTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtg                     >  1:409134/1‑71 (MQ=255)
                 acTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtg                     >  1:3929789/1‑71 (MQ=255)
                 acTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtg                     >  1:4136607/1‑71 (MQ=255)
                 acTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtg                     >  1:2415091/1‑71 (MQ=255)
                 acTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtg                     >  1:3266035/1‑71 (MQ=255)
                 acTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCtg                     >  1:2696424/1‑71 (MQ=255)
                 acTCGTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCt                      >  1:1502579/1‑70 (MQ=255)
                     gTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTgg                  <  1:1257780/70‑1 (MQ=255)
                     gTTTACGCATAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGAGCTGTgg                  <  1:3226217/70‑1 (MQ=255)
                             aTAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAAt         >  1:2225179/1‑71 (MQ=255)
                             aTAGCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAAt         >  1:2162534/1‑71 (MQ=255)
                                gCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGc                 >  1:3861479/1‑60 (MQ=255)
                                gCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGc                 >  1:295546/1‑60 (MQ=255)
                                gCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTaa      >  1:845574/1‑71 (MQ=255)
                                gCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTa       <  1:2435010/70‑1 (MQ=255)
                                gCTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTa       >  1:2196316/1‑70 (MQ=255)
                                 cTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAAt     <  1:2460956/71‑1 (MQ=255)
                                 cTATGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAAt     <  1:3748265/71‑1 (MQ=255)
                                    tGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:313036/71‑1 (MQ=255)
                                    tGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:590706/71‑1 (MQ=255)
                                    tGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:5297/71‑1 (MQ=255)
                                    tGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:489192/71‑1 (MQ=255)
                                    tGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:1137290/71‑1 (MQ=255)
                                    tGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:4041749/71‑1 (MQ=255)
                                    tGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:3992003/71‑1 (MQ=255)
                                    tGCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:3708381/71‑1 (MQ=255)
                                     gCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  >  1:2130566/1‑70 (MQ=255)
                                     gCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  >  1:1103718/1‑70 (MQ=255)
                                     gCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  >  1:3165635/1‑70 (MQ=255)
                                     gCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCtt  <  1:1008701/70‑1 (MQ=255)
                                     gCAGGAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCt   >  1:2317042/1‑69 (MQ=255)
                                         |                                                                 
ACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTT  >  W3110S.gb/2092324‑2092430

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: