Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,350,543 G→A 40.0% 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,5430GA40.0% 26.1 / 44.7 40V337I (GTA→ATA) nrdAribonucleoside diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (4/20);  new base A (10/6);  total (14/26)
Fisher's exact test for biased strand distribution p-value = 5.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GAAAGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTT  >  W3110S.gb/2350474‑2350601
                                                                     |                                                          
gAAAGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACgggg                                                            >  1:73387/1‑70 (MQ=255)
            gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                <  1:833218/70‑1 (MQ=255)
            gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                <  1:656277/70‑1 (MQ=255)
            gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                <  1:1331176/70‑1 (MQ=255)
            gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                <  1:1333138/70‑1 (MQ=255)
            gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                <  1:1433444/70‑1 (MQ=255)
            gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                <  1:1489227/70‑1 (MQ=255)
            gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                <  1:233641/70‑1 (MQ=255)
            gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaaca                                                <  1:1673753/70‑1 (MQ=255)
                 gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGa                                          >  1:1358884/1‑71 (MQ=255)
                 gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGa                                          >  1:1456718/1‑71 (MQ=255)
                 gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGa                                          >  1:568629/1‑71 (MQ=255)
                 gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGa                                          >  1:1890965/1‑71 (MQ=255)
                     aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGta                                      <  1:1097863/71‑1 (MQ=255)
                     aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGta                                      <  1:868401/71‑1 (MQ=255)
                     aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGta                                      <  1:104930/71‑1 (MQ=255)
                      acaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGtat                                     <  1:1153922/71‑1 (MQ=255)
                      acaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGtat                                     <  1:571271/71‑1 (MQ=255)
                                tgtgGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGta                                      <  1:663406/60‑1 (MQ=255)
                                       ggCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGtat                                     <  1:234741/54‑1 (MQ=255)
                                           aCCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGaa                <  1:283827/71‑1 (MQ=255)
                                           aCCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGaa                <  1:1294307/71‑1 (MQ=255)
                                             cgcgTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGta                                      >  1:786306/1‑47 (MQ=255)
                                             cgcgTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGta                                      >  1:1521593/1‑47 (MQ=255)
                                             cgcgTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTg                  >  1:837566/1‑67 (MQ=255)
                                             cgcgTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTg                  >  1:727481/1‑67 (MQ=255)
                                              gcgtgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:1794452/71‑1 (MQ=255)
                                              gcgtgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:745771/71‑1 (MQ=255)
                                              gcgtgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:510857/71‑1 (MQ=255)
                                               cgtgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:676946/70‑1 (MQ=255)
                                                gtgcgtCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGatat           <  1:654513/71‑1 (MQ=255)
                                                gtgcgtCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGatat           <  1:1938862/71‑1 (MQ=255)
                                                 tgcgtCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGatat           <  1:595065/70‑1 (MQ=255)
                                                    gtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCAcc       <  1:218418/71‑1 (MQ=255)
                                                    gtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCAcc       <  1:1785788/71‑1 (MQ=255)
                                                          tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt  >  1:481064/1‑70 (MQ=255)
                                                          tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt  >  1:1609918/1‑70 (MQ=255)
                                                          tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt  >  1:729006/1‑70 (MQ=255)
                                                          tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt  >  1:1715049/1‑70 (MQ=255)
                                                          tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt  >  1:1105847/1‑70 (MQ=255)
                                                          tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt  >  1:1828498/1‑70 (MQ=255)
                                                                     |                                                          
GAAAGCCTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTT  >  W3110S.gb/2350474‑2350601

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