Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 891,153 C→T 39.0% 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,1530CT39.0% 31.4 / 42.5 41R8H (CGT→CAT) grxAglutaredoxin 1, redox coenzyme for ribonucleotide reductase (RNR1a)
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base T (9/7);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATG  >  W3110S.gb/891090‑891222
                                                               |                                                                     
tcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCaaa                                                                <  1:24328/71‑1 (MQ=255)
      tCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAAc                                                          >  1:4638300/1‑71 (MQ=255)
       cATTGCTCAATTTCTCAGCCAGATCTTTTG‑‑CGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGt                                                       >  1:4126924/1‑71 (MQ=255)
           gCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGttt                                                     >  1:4986907/1‑71 (MQ=255)
            cTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTg                                                    <  1:2463214/71‑1 (MQ=255)
            cTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATa                                                            >  1:2345411/1‑63 (MQ=255)
                   tctcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCattatt                                             >  1:1060522/1‑71 (MQ=255)
                   tctcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCattatt                                             >  1:4188044/1‑71 (MQ=255)
                   tctcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCattatt                                             >  1:2893939/1‑71 (MQ=255)
                   tctcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCattat                                              <  1:623430/70‑1 (MQ=255)
                   tctcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCattat                                              <  1:3995511/70‑1 (MQ=255)
                       aGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTtctc                                         <  1:1229288/71‑1 (MQ=255)
                       aGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTAttt                                            <  1:2159699/68‑1 (MQ=255)
                       aGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTAttt                                            <  1:4233395/68‑1 (MQ=255)
                          cAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTctcc                                      >  1:1451639/1‑71 (MQ=255)
                             aTCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTtctctc                                       <  1:1147775/67‑1 (MQ=255)
                                    gcacgcacACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGAtt                             <  1:1420551/70‑1 (MQ=255)
                                    gcacgcacACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGAtt                             <  1:254015/70‑1 (MQ=255)
                                      acgcacACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTAt                          <  1:798683/71‑1 (MQ=255)
                                      acgcacACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTAt                          <  1:1357226/71‑1 (MQ=255)
                                      acgcacACAGTAAGGGCAACCCGAACAACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTAt                          <  1:535537/71‑1 (MQ=255)
                                        gcacACAGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGc                        >  1:2525965/1‑71 (MQ=255)
                                        gcacACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGc                        >  1:3604759/1‑71 (MQ=255)
                                        gcacACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGc                        >  1:2171876/1‑71 (MQ=255)
                                            acaGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTg                     >  1:3164023/1‑70 (MQ=255)
                                            acaGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGt                    <  1:4173868/71‑1 (MQ=255)
                                              aGTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTaa                  >  1:4211103/1‑71 (MQ=255)
                                               gTAAGGGCAACCCGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAAt                 <  1:2341891/71‑1 (MQ=255)
                                               gTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAAt                 <  1:2239555/71‑1 (MQ=255)
                                                 aaGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGc                        >  1:3515662/1‑62 (MQ=255)
                                                  aGGGCAACCCGAACGACCAAAAATAACGGTTTGCattat                                              >  1:1747393/1‑39 (MQ=255)
                                                      cAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTa                   <  1:3770950/62‑1 (MQ=255)
                                                         cccGAATGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCt        >  1:321035/1‑70 (MQ=255)
                                                           cGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCAt                                  >  1:583133/1‑42 (MQ=255)
                                                            gAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTaaaa    >  1:2251341/1‑71 (MQ=255)
                                                            gAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTaaaa    >  1:4852200/1‑71 (MQ=255)
                                                            gAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTaa      >  1:1605793/1‑69 (MQ=255)
                                                             aaCGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAAt   <  1:725667/71‑1 (MQ=255)
                                                              aCGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAAt   <  1:1271888/70‑1 (MQ=255)
                                                              aCGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATg  <  1:3882374/71‑1 (MQ=255)
                                                              aCGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATg  <  1:1959798/71‑1 (MQ=255)
                                                               cGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATg  >  1:1020956/1‑70 (MQ=255)
                                                               |                                                                     
TCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATG  >  W3110S.gb/891090‑891222

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: