Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,350,543 G→A 80.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,5430GA80.0% 138.8 / 31.7 65V337I (GTA→ATA) nrdAribonucleoside diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (5/8);  new base A (23/29);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGC  >  W3110S.gb/2350485‑2350605
                                                          |                                                              
ggTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATca                                                       >  1:1846289/1‑68 (MQ=255)
ggTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATca                                                       >  1:3460862/1‑68 (MQ=255)
   gtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATcaacaa                                                   >  1:1490950/1‑69 (MQ=255)
   gtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATcaaca                                                    <  1:2047345/68‑1 (MQ=255)
      gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACt                                                <  1:1948629/69‑1 (MQ=255)
      gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACt                                                >  1:2197172/1‑69 (MQ=255)
      gAAAAACAACCGTGGTGTGGAAGGCAACAGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACt                                                >  1:2614821/1‑69 (MQ=255)
       aaaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACt                                                <  1:280856/68‑1 (MQ=255)
        aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATcaacaa                                                   <  1:2412709/64‑1 (MQ=255)
        aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATcaacaa                                                   <  1:2575583/64‑1 (MQ=255)
        aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTg                                               <  1:53522/68‑1 (MQ=255)
        aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTg                                               <  1:128753/68‑1 (MQ=255)
        aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTA‑‑‑GGTACAAATCAACAAACTGATg                                            <  1:1757499/68‑1 (MQ=255)
          aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                            <  1:1336480/69‑1 (MQ=255)
          aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                            >  1:1405542/1‑69 (MQ=255)
          aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                            >  1:3702010/1‑69 (MQ=255)
          aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                            <  1:2469406/69‑1 (MQ=255)
          aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            <  1:2195233/69‑1 (MQ=255)
          aacaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            <  1:1940831/69‑1 (MQ=255)
           acaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            >  1:1649181/1‑68 (MQ=255)
           acaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            >  1:3001049/1‑68 (MQ=255)
           acaacCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            >  1:46658/1‑68 (MQ=255)
                cGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATAc                                       <  1:1186768/68‑1 (MQ=255)
                    gtgtGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACt                                                <  1:2783141/55‑1 (MQ=255)
                       tgGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTc                                  <  1:2813158/66‑1 (MQ=255)
                       tgGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTc                                  <  1:2644022/66‑1 (MQ=255)
                       tgGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTc                                  <  1:1584555/66‑1 (MQ=255)
                         gAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTctgct                              >  1:3304083/1‑68 (MQ=255)
                         gAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTctgct                              >  1:3502960/1‑68 (MQ=255)
                         gAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTctgct                              >  1:2865319/1‑68 (MQ=255)
                           aGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTg                                               >  1:2115185/1‑49 (MQ=255)
                           aGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTg                                               >  1:3521505/1‑49 (MQ=255)
                           aGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTg                                               >  1:1659907/1‑49 (MQ=255)
                           aGGCAACCGCGTGCGTCATATGGACTAAGGGATACAAATCAACAAACTg                                               >  1:1158942/1‑49 (MQ=255)
                              cAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            >  1:3237329/1‑49 (MQ=255)
                              cAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            <  1:1407149/49‑1 (MQ=255)
                              cAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            >  1:1681584/1‑49 (MQ=255)
                              cAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            <  1:2719341/49‑1 (MQ=255)
                              cAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            >  1:2601416/1‑49 (MQ=255)
                              cAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                            <  1:2294950/49‑1 (MQ=255)
                               aaCCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAgg                        >  1:3563393/1‑68 (MQ=255)
                                  cgcgTGCGTCATATGGACTACGGGGTACAAATCAACAACCTGATg                                            <  1:1963/45‑1 (MQ=255)
                                  cgcgTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATg                                            <  1:3014654/45‑1 (MQ=255)
                                  cgcgTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAgg                        >  1:7444/1‑65 (MQ=255)
                                    cgtgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                  <  1:3037031/69‑1 (MQ=255)
                                    cgtgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa                  <  1:2001491/69‑1 (MQ=255)
                                       gcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGatat               <  1:2672417/69‑1 (MQ=255)
                                        cgtCATATGGACTACGGGATACAAATCAACAAACt                                                >  1:605948/1‑35 (MQ=255)
                                        cgtCATATGGACTACGGGATACAAATCAACAAACt                                                >  1:1302012/1‑35 (MQ=255)
                                        cgtCATATGGACTACGGGATACAAATCAACAAACt                                                >  1:3085332/1‑35 (MQ=255)
                                        cgtCATATGGACTACGGGATACAAATCAACAAACt                                                >  1:2791640/1‑35 (MQ=255)
                                         gtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATc              <  1:2968647/68‑1 (MQ=255)
                                          tCATATGGACTACGGGATACAAATCAACAAACTGATg                                            <  1:3139183/37‑1 (MQ=255)
                                            atatGGACTACGGGATACAAATCAACAAACTGATg                                            <  1:1589195/35‑1 (MQ=255)
                                            atatGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCAccc          <  1:1785667/69‑1 (MQ=255)
                                              atGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTg        >  1:3066306/1‑69 (MQ=255)
                                               tGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGaaa                          <  1:643750/50‑1 (MQ=255)
                                               tGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGaaa                          <  1:819722/50‑1 (MQ=255)
                                                 gACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTc     >  1:218628/1‑69 (MQ=255)
                                                 gACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCAc            <  1:1436659/62‑1 (MQ=255)
                                                 gACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt      <  1:3124213/68‑1 (MQ=255)
                                                   cTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTaaa                    <  1:1418430/52‑1 (MQ=255)
                                                   cTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGaa                    <  1:3223773/52‑1 (MQ=255)
                                                   cTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGaa                    <  1:3465990/52‑1 (MQ=255)
                                                    tACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGc  <  1:655136/69‑1 (MQ=255)
                                                          gTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGatat               >  1:3252325/1‑50 (MQ=255)
                                                          |                                                              
GGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTTCAGC  >  W3110S.gb/2350485‑2350605

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: