Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,691,177 G→A 14.3% A1008V (GCG→GTG)  purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,691,1770GA14.3% 71.0 / 3.8 28A1008V (GCG→GTG) purLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base G (14/9);  new base A (2/2);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGTTTCTGC  >  W3110S.gb/2691121‑2691244
                                                        |                                                                   
cGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTc                                                           >  1:545294/1‑67 (MQ=255)
cGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCt                                                         >  1:1039280/1‑69 (MQ=255)
cGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCt                                                         >  1:291661/1‑69 (MQ=255)
cGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCt                                                         >  1:1371887/1‑69 (MQ=255)
cGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCt                                                         >  1:2564232/1‑69 (MQ=255)
cGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTg                                                        <  1:2507447/70‑1 (MQ=255)
cGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTg                                                        <  1:744636/70‑1 (MQ=255)
      tatCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCt                                                                  <  1:1037078/54‑1 (MQ=255)
      tatCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTg                                                        <  1:1217272/64‑1 (MQ=255)
        tCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGAt                                                      >  1:232016/1‑64 (MQ=255)
        tCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGAt                                                      >  1:2393698/1‑64 (MQ=255)
          gAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGCTCGGcaca                                               >  1:2308740/1‑69 (MQ=255)
          gAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGcaca                                               >  1:1247210/1‑69 (MQ=255)
          gAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGcac                                                >  1:1300725/1‑68 (MQ=255)
          gAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGcac                                                >  1:334640/1‑68 (MQ=255)
          gAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACt                                             >  1:550683/1‑71 (MQ=255)
          gAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACt                                             >  1:409088/1‑71 (MQ=255)
                   ttttACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCgg                                                   <  1:875954/56‑1 (MQ=255)
                          ttCAGGCCCGGATCGGCGTCGTTAGATTTCACCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCa                              <  1:2091247/70‑1 (MQ=255)
                          ttCAGGCCCGGATCGGCGTCGTTAGATTTCACCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCa                              <  1:1616025/70‑1 (MQ=255)
                            cAGGCCCGGATCGGCGTCGTTAGATTTCACCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAgg                            >  1:987958/1‑70 (MQ=255)
                            cAGGACCGGATCGGCGTCGTTAGATTTCACCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAg                             >  1:708285/1‑69 (MQ=255)
                                 ccGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGc                      <  1:1573312/71‑1 (MQ=255)
                                     aTCGGCGTCGTTAGATTTCTCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCt                  <  1:629949/71‑1 (MQ=255)
                                           gtcgtTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAgg            <  1:2129193/71‑1 (MQ=255)
                                               ttAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACg                                <  1:541746/47‑1 (MQ=255)
                                                 aGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGtt       <  1:2163538/70‑1 (MQ=255)
                                                     ttCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGTTTCTGc  >  1:83234/1‑71 (MQ=255)
                                                        |                                                                   
CGTTGATATCGAACGACAGTTTTACATTCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGTTTCTGC  >  W3110S.gb/2691121‑2691244

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: