Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,350,543 G→A 35.8% V337I (GTA→ATA)  nrdA → ribonucleoside diphosphate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,350,5430GA35.8% 51.9 / 49.7 53V337I (GTA→ATA) nrdAribonucleoside diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (23/11);  new base A (11/8);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGT  >  W3110S.gb/2350476‑2350609
                                                                   |                                                                  
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:1065548/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:648336/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:618788/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:3240181/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:3099931/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:3033023/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:2853616/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:1520467/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:2427502/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:2208521/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:2113075/1‑70 (MQ=255)
aaGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:1740338/1‑70 (MQ=255)
aaGCCTGATGGTGTTGAAAAACAAACGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTa                                                                  >  1:1375309/1‑70 (MQ=255)
          gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                        <  1:99629/70‑1 (MQ=255)
          gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                        <  1:909768/70‑1 (MQ=255)
          gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                        <  1:196137/70‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                                <  1:3244698/71‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                                <  1:1876099/71‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                                <  1:2054203/71‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                                <  1:1313772/71‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                                <  1:1313785/71‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                                <  1:3025608/71‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGGACAAATCAACAAACTGATg                                                <  1:1821087/71‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                                <  1:1125077/71‑1 (MQ=255)
                 aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                                <  1:525006/71‑1 (MQ=255)
                  aaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                                <  1:510650/70‑1 (MQ=255)
                  aaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                                <  1:386340/70‑1 (MQ=255)
                  aaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                                <  1:99244/70‑1 (MQ=255)
                  aaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                                <  1:3102732/70‑1 (MQ=255)
                  aaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                                <  1:3221348/70‑1 (MQ=255)
                       acCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATAcc                                          >  1:106846/1‑71 (MQ=255)
                       acCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATAcc                                          >  1:203543/1‑71 (MQ=255)
                            ggtgTGGAAGGCAACCGCGTGAGTCATATGGACTACTGGATACAAATCAACAAACTGATg                                                <  1:694110/60‑1 (MQ=255)
                                           cgcgTGCGTCATATGGACTACGGGGTACAAATCATCAAACTGATg                                                <  1:780904/45‑1 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGt  >  1:860538/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGt  >  1:3284816/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGt  >  1:3201338/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGt  >  1:1372461/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGt  >  1:171940/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGt  >  1:2413328/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGt  >  1:1722846/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGt  >  1:1736324/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCAGTTCAGCCCGt  >  1:2781667/1‑71 (MQ=255)
                                                               cGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGAGAAGATATCACCCTGTTCAGCCCg   >  1:889038/1‑70 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTTTGCTGAAAGGTGAAGa                      >  1:2446576/1‑51 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGa                         >  1:756845/1‑48 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                      >  1:1161723/1‑51 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                      >  1:3092838/1‑51 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                      >  1:2966303/1‑51 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                      >  1:2479814/1‑51 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                      >  1:1096382/1‑51 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                      >  1:881951/1‑51 (MQ=255)
                                                               cGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                      >  1:1024731/1‑51 (MQ=255)
                                                                   |                                                                  
AAGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGCCCGT  >  W3110S.gb/2350476‑2350609

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: