Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,350,543 G→A 100% 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,5430GA94.5% 155.5 / 1.1 54V337I (GTA→ATA) nrdAribonucleoside diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (22/29);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTT  >  W3110S.gb/2350480‑2350601
                                                               |                                                          
ctgctgGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAAt                                                     >  1:2606086/1‑71 (MQ=255)
   ctgGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATcaa                                                  >  1:714757/1‑71 (MQ=255)
   ctgGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATcaa                                                  >  1:1799119/1‑71 (MQ=255)
      gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATcaaca                                                <  1:1228456/70‑1 (MQ=255)
      gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATcaaca                                                <  1:2776189/70‑1 (MQ=255)
      gtgtTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATcaaca                                                <  1:1563385/70‑1 (MQ=255)
         ttGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACt                                            <  1:3124911/71‑1 (MQ=255)
         ttGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACt                                            <  1:2010727/71‑1 (MQ=255)
          tGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACt                                            <  1:930778/70‑1 (MQ=255)
           gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTg                                           <  1:257834/70‑1 (MQ=255)
           gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTg                                           <  1:2246380/70‑1 (MQ=255)
           gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTg                                           <  1:867758/70‑1 (MQ=255)
           gAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGa                                          >  1:1206173/1‑71 (MQ=255)
             aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGTGGTACAAATCAACAAACTGATg                                        <  1:3140658/71‑1 (MQ=255)
             aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCGACAAACTGATg                                        <  1:137297/71‑1 (MQ=255)
             aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                        <  1:471140/71‑1 (MQ=255)
             aaaaCAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATg                                        <  1:669573/71‑1 (MQ=255)
                              gAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGtata                                    >  1:545695/1‑58 (MQ=255)
                              gAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGtata                                    >  1:2690049/1‑58 (MQ=255)
                              gAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGtata                                    >  1:2561669/1‑58 (MQ=255)
                              gAAGGCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGaa                       <  1:435589/71‑1 (MQ=255)
                                 ggCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATAcc                                  >  1:1472337/1‑57 (MQ=255)
                                 ggCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATAcc                                  >  1:522788/1‑57 (MQ=255)
                                 ggCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATAcc                                  >  1:557000/1‑57 (MQ=255)
                                 ggCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATAcc                                  >  1:1578953/1‑57 (MQ=255)
                                 ggCAACCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATAcc                                  >  1:885447/1‑57 (MQ=255)
                                     aCCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGta                                      <  1:1024917/49‑1 (MQ=255)
                                     aCCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGta                                      <  1:279153/49‑1 (MQ=255)
                                     aCCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGta                                      <  1:3027423/49‑1 (MQ=255)
                                     aCCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGta                                      <  1:401397/49‑1 (MQ=255)
                                     aCCGCGTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGaa                <  1:2529875/71‑1 (MQ=255)
                                       cgcgTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGta                                      >  1:2697798/1‑47 (MQ=255)
                                       cgcgTGCGTCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa              <  1:3142270/71‑1 (MQ=255)
                                        gcgtgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAg               >  1:2354011/1‑69 (MQ=255)
                                        gcgtgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa              >  1:2349893/1‑70 (MQ=255)
                                           tgcgtCATATGGACTACGGGATACAAATCAACAAACTGATGta                                      <  1:178167/43‑1 (MQ=255)
                                                  tatGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGa                        <  1:2873001/50‑1 (MQ=255)
                                                  tatGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGa                        <  1:2011980/50‑1 (MQ=255)
                                                  tatGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGa                        <  1:1374530/50‑1 (MQ=255)
                                                  tatGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGa                        <  1:2964166/50‑1 (MQ=255)
                                                   atGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt  >  1:33148/1‑71 (MQ=255)
                                                   atGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGt   >  1:1277027/1‑70 (MQ=255)
                                                    tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:1476081/59‑1 (MQ=255)
                                                    tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:2802578/59‑1 (MQ=255)
                                                    tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:1524787/59‑1 (MQ=255)
                                                    tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:3466368/59‑1 (MQ=255)
                                                    tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGat             <  1:3461858/59‑1 (MQ=255)
                                                    tGGACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGtt  >  1:3440977/1‑70 (MQ=255)
                                                      gACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa              >  1:3117738/1‑56 (MQ=255)
                                                      gACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa              >  1:2793256/1‑56 (MQ=255)
                                                      gACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa              >  1:59043/1‑56 (MQ=255)
                                                      gACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa              >  1:2561739/1‑56 (MQ=255)
                                                      gACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa              >  1:2158345/1‑56 (MQ=255)
                                                      gACTACGGGATACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGa              >  1:1445202/1‑56 (MQ=255)
                                                               |                                                          
CTGCTGGTGTTGAAAAACAACCGTGGTGTGGAAGGCAACCGCGTGCGTCATATGGACTACGGGGTACAAATCAACAAACTGATGTATACCCGTCTGCTGAAAGGTGAAGATATCACCCTGTT  >  W3110S.gb/2350480‑2350601

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