Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,473,895 A→G 71.4% E162G (GAG→GGG)  trmA → tRNA (uracil‑5‑)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,473,8950AG71.4% 62.7 / 29.2 42E162G (GAG→GGG) trmAtRNA (uracil‑5‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base G (17/13);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CGCGATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAG  >  W3110S.gb/3473828‑3473963
                                                                   |                                                                    
cgcgATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCt                                                                   >  1:2709679/1‑71 (MQ=255)
       cACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAg                                                            >  1:553833/1‑71 (MQ=255)
       cACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAg                                                            <  1:273863/71‑1 (MQ=255)
       cACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAg                                                            <  1:1954498/71‑1 (MQ=255)
          tGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATc                                                              <  1:824313/66‑1 (MQ=255)
            cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACTAAAACCTAAATCGAGCTGGATCAGGATTa                                                       <  1:2683055/71‑1 (MQ=255)
            cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGAtt                                                        >  1:2746583/1‑70 (MQ=255)
            cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTa                                                       >  1:1220863/1‑71 (MQ=255)
            cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTa                                                       <  1:261753/71‑1 (MQ=255)
            cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTa                                                       >  1:2596303/1‑71 (MQ=255)
            cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTa                                                       <  1:1614877/71‑1 (MQ=255)
                    gAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATg                                               <  1:78377/71‑1 (MQ=255)
                       tCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAAc                                            >  1:1585510/1‑71 (MQ=255)
                        cTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGa                                                          >  1:1237045/1‑56 (MQ=255)
                        cTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGa                                                          >  1:1971792/1‑56 (MQ=255)
                         tGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGt                                          >  1:96775/1‑71 (MQ=255)
                              gtgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGaa                                             >  1:166411/1‑63 (MQ=255)
                               tgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTgg                                                                 <  1:57469/42‑1 (MQ=255)
                               tgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGa                                                 >  1:617624/1‑58 (MQ=255)
                               tgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACg                                           <  1:72337/64‑1 (MQ=255)
                               tgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACg                                           <  1:1399398/64‑1 (MQ=255)
                                   tCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGAt                                                         <  1:2432326/46‑1 (MQ=255)
                                    cTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGt                                          >  1:1799806/1‑60 (MQ=255)
                                     tGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCg                                >  1:569648/1‑69 (MQ=255)
                                     tGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCg                                >  1:566358/1‑69 (MQ=255)
                                     tGATTGGACGGGCAACGAAAACCAAAATAGGTCTGGAta                                                              >  1:1698628/1‑38 (MQ=21)
                                      aaTTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGaa                                             >  1:623519/2‑55 (MQ=255)
                                       aTTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAgg                            >  1:1066419/1‑71 (MQ=255)
                                         tGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTa                                                       >  1:333794/1‑42 (MQ=255)
                                          ggTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa                          <  1:2453458/70‑1 (MQ=255)
                                              gggCAACGAAAACCAAAATCGGGCTGGATCATGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga                      <  1:2018838/70‑1 (MQ=255)
                                              gggCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga                      <  1:1016961/70‑1 (MQ=255)
                                                   aCGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa                          <  1:529890/61‑1 (MQ=255)
                                                     gAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATc               >  1:1811202/1‑70 (MQ=255)
                                                     gAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATc               >  1:2515356/1‑70 (MQ=255)
                                                     gAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCt              <  1:2449333/71‑1 (MQ=255)
                                                         aCCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTAc            <  1:1204936/69‑1 (MQ=255)
                                                            aaaaTCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga                      <  1:758166/56‑1 (MQ=255)
                                                             aaaTCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga                      <  1:519022/55‑1 (MQ=255)
                                                             aaaTCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga                      <  1:809293/55‑1 (MQ=255)
                                                              aaTCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga                 <  1:46842/59‑1 (MQ=255)
                                                                  gAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAg  >  1:2791733/1‑70 (MQ=255)
                                                                   |                                                                    
CGCGATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAG  >  W3110S.gb/3473828‑3473963

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: