Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 891,153 C→T 61.8% R8H (CGT→CAT)  grxA ← glutaredoxin 1, redox coenzyme for ribonucleotide reductase (RNR1a)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb891,1530CT61.8% 22.0 / 35.1 34R8H (CGT→CAT) grxAglutaredoxin 1, redox coenzyme for ribonucleotide reductase (RNR1a)
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (11/10);  total (17/17)
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.95e-01

AATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTA  >  W3110S.gb/891085‑891217
                                                                    |                                                                
aaTCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGa                                                                >  1:17050/1‑71 (MQ=255)
aaTCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGa                                                                >  1:1672194/1‑71 (MQ=255)
        cgTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAAt                                                        <  1:959660/71‑1 (MQ=255)
         gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAAt                                                        <  1:957667/70‑1 (MQ=255)
            cATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACg                                                    >  1:296779/1‑71 (MQ=255)
            cATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACg                                                    >  1:1338946/1‑71 (MQ=255)
            cATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACg                                                    >  1:1595852/1‑71 (MQ=255)
             aTTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACgg                                                   <  1:519130/71‑1 (MQ=255)
             aTTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACgg                                                   <  1:274370/71‑1 (MQ=255)
                gCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGtt                                                 <  1:143393/70‑1 (MQ=255)
                  tCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACg                                                    <  1:1138731/65‑1 (MQ=255)
                      tttCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGt                                                  <  1:1207713/63‑1 (MQ=255)
                       ttCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGttt                                                <  1:864932/64‑1 (MQ=255)
                       ttCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCattat                                         >  1:1133398/1‑71 (MQ=255)
                                  aTCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTg                                               <  1:1372708/54‑1 (MQ=255)
                                  aTCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTg                                               <  1:756112/54‑1 (MQ=255)
                                    cTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATaa                                                      >  1:1385838/1‑45 (MQ=255)
                                     ttttGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTAttt                                       <  1:851540/59‑1 (MQ=255)
                                     ttttGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTAttt                                       <  1:337505/59‑1 (MQ=255)
                                       ttGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGtt                                                 >  1:1357532/1‑47 (MQ=255)
                                       ttGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCa                              <  1:1226779/66‑1 (MQ=255)
                                       ttGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGAt                         >  1:1352992/1‑71 (MQ=255)
                                          cacgcacACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGAttt                       <  1:1731619/70‑1 (MQ=255)
                                                 acaGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTg                >  1:1212785/1‑70 (MQ=255)
                                                 acaGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTg                >  1:475595/1‑70 (MQ=255)
                                                 acaGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTg                >  1:1431080/1‑70 (MQ=255)
                                                  caGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGt               >  1:1519329/1‑70 (MQ=255)
                                                   aGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTaa             >  1:1585547/1‑71 (MQ=255)
                                                       aGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGAt         >  1:817831/1‑71 (MQ=255)
                                                              cccGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCt   >  1:1295059/1‑70 (MQ=255)
                                                               ccGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATa                            <  1:213493/44‑1 (MQ=255)
                                                                cGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTCTAATGATCACGCTa  >  1:1559803/1‑69 (MQ=255)
                                                                  aaCGACCAAAAATAACGGTTTGCATTATTtctctc                                  <  1:1217581/35‑1 (MQ=255)
                                                                  aaCGACCAAAAATAACGGTTTGCATTATTtctctc                                  <  1:1007016/35‑1 (MQ=255)
                                                                    |                                                                
AATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTA  >  W3110S.gb/891085‑891217

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: