Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,909,021 G→A 13.8% N107N (AAC→AAT mazG ← nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,909,0210GA13.8% 145.0 / 12.6 58N107N (AAC→AATmazGnucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (29/21);  new base A (5/3);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATCGCTAATAGCAGCGC  >  W3110S.gb/2908952‑2909088
                                                                     |                                                                   
cGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGt                                                                    >  1:222586/1‑71 (MQ=255)
      ctgcGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCgg                                                              <  1:667319/71‑1 (MQ=255)
         cgcTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAg                                                           <  1:819268/71‑1 (MQ=255)
              tCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGTCAAGCACTTCACTACTGTTTTCGGCAGAACt                                                       <  1:1601750/70‑1 (MQ=255)
                     gcgcTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCaa                                               <  1:1961861/71‑1 (MQ=255)
                        ctctTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCaaaa                                             <  1:1637604/70‑1 (MQ=255)
                           ttCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTc                                                                <  1:1595570/48‑1 (MQ=255)
                           ttCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACt                                                       <  1:2054089/57‑1 (MQ=255)
                           ttCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCaaaa                                             <  1:3831333/67‑1 (MQ=255)
                           ttCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACAt                                         >  1:1981790/1‑71 (MQ=255)
                             cGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGaa                                                         >  1:516816/1‑53 (MQ=255)
                             cGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGaa                                                         >  1:2861893/1‑53 (MQ=255)
                             cGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGc                                       >  1:3754057/1‑71 (MQ=255)
                             cGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGc                                       >  1:2471929/1‑71 (MQ=255)
                                ttttGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCaa                                               <  1:1124226/60‑1 (MQ=255)
                                  ttGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATg                                  >  1:42847/1‑71 (MQ=255)
                                   tGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCa                                                            >  1:3676412/1‑44 (MQ=255)
                                   tGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATgg                                 >  1:2501363/1‑71 (MQ=255)
                                    gATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGc                                >  1:1692569/1‑71 (MQ=255)
                                     aTTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGc                                <  1:2297952/70‑1 (MQ=255)
                                        tGCTCCCAAGGGGCAAGCACTTCACTACTATTTTCGGCAGa                                                          <  1:588922/41‑1 (MQ=38)
                                        tGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGa                              >  1:3410939/1‑69 (MQ=255)
                                        tGCTCCCAACGGGCAAGCACTTCACTACTATTTTCGGCAGa                                                          <  1:3002169/41‑1 (MQ=255)
                                        tGCTCCCAACGGGCAAGCACTTCACTACTATTTTCGGCAGa                                                          <  1:2871104/41‑1 (MQ=255)
                                        tGCTCCCAACGGGCAAGCACTTCACTACTATTTTCGGCAGAACTGTCAg                                                  >  1:601109/1‑49 (MQ=255)
                                        tGCTCCCAACGGGCAAGCACTTCACTACTATTTTCGGCAGAACTGTCAg                                                  >  1:1567707/1‑49 (MQ=255)
                                        tGCTCCCAACGGGCAAGCACTTCACTACTATTTTCGGCAGAACTGTCAg                                                  >  1:2242030/1‑49 (MQ=255)
                                        tGCTCCCAACGGGCAAGCACTTCACTACTATTTTCGGCAGAACTGTCAg                                                  >  1:1268963/1‑49 (MQ=255)
                                        tGCTCCCAACGGGCAAGCACTTCACTACTATTTTCGGCAGAACTGTCAg                                                  >  1:755666/1‑49 (MQ=255)
                                         gCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACg                            <  1:2463288/70‑1 (MQ=255)
                                          cTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGct                          >  1:3319076/1‑71 (MQ=255)
                                          cTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGc                           >  1:1957864/1‑70 (MQ=255)
                                             ccAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGctct                        >  1:770316/1‑70 (MQ=255)
                                               aaCGGGCAAGCACTTCACTACTGTTTTCGGCAGAACt                                                       <  1:2280270/37‑1 (MQ=255)
                                               aaCGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCaaaaa                                            <  1:2320820/48‑1 (MQ=255)
                                                aCGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTAt                                                     >  1:2423809/1‑38 (MQ=255)
                                                aCGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTAt                                                     >  1:382177/1‑38 (MQ=255)
                                                 cGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTAt                                                     >  1:3244492/1‑37 (MQ=255)
                                                 cGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAAttt                   <  1:31270/71‑1 (MQ=255)
                                                  gggCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCaaaaa                                            <  1:1210401/45‑1 (MQ=255)
                                                  gggCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGa                                    >  1:2809483/1‑53 (MQ=255)
                                                    gCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATc                >  1:1080475/1‑71 (MQ=255)
                                                       aGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGa                                    >  1:3764826/1‑48 (MQ=255)
                                                        gCACTTCACTACTGTTTTCGGCAGAACTATCAGCaaaaa                                            <  1:1898271/39‑1 (MQ=255)
                                                        gCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATCGCt             >  1:2034606/1‑70 (MQ=255)
                                                         cACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAACCATGCGGATGGCGACGCTCTAATTTATCGCTaa           <  1:1880184/71‑1 (MQ=255)
                                                         cACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATCGCTaa           <  1:1131853/71‑1 (MQ=255)
                                                          aCTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATCGCTaa           >  1:2099241/1‑70 (MQ=255)
                                                          aCTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATCGCTAAt          <  1:204144/71‑1 (MQ=255)
                                                            ttCACTACTGTTTTCGGCAGAACTATCAGCaaaaa                                            <  1:1508947/35‑1 (MQ=255)
                                                             tCACTACTGTTTTCGGCAGAACTATCAGCAAAAACa                                          >  1:3558978/1‑36 (MQ=255)
                                                             tCACTACTGTTTTCGGCAGAACTATCAGCAAAAACa                                          >  1:153789/1‑36 (MQ=255)
                                                             tCACTACTGTTTTCGGCAGAACTATCAGCAAAAACa                                          >  1:2011163/1‑36 (MQ=255)
                                                             tCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATCGCTAATag        >  1:1557915/1‑70 (MQ=255)
                                                              cACTACTGTTTTCGGCAGAACTATCAGCAAAAAct                                          >  1:1936061/1‑34 (MQ=255)
                                                                   ctGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAAtt                    >  1:1194072/1‑52 (MQ=255)
                                                                   ctGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATCGCTAATAGCAgcgc  >  1:530816/1‑70 (MQ=255)
                                                                    tGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAAttt                   <  1:1800668/52‑1 (MQ=255)
                                                                     |                                                                   
CGAATGCTGCGCTTTCTGCGCGCGCTCTTCGGTTTTGATTTGCTCCCAACGGGCAAGCACTTCACTACTGTTTTCGGCAGAACTATCAGCAAAAACATGCGGATGGCGACGCTCTAATTTATCGCTAATAGCAGCGC  >  W3110S.gb/2908952‑2909088

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: