Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,413,273 G→T 12.7% R130R (CGC→CGA ydaO ← predicted C32 tRNA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,413,2730GT12.7% 142.5 / 5.4 47R130R (CGC→CGAydaOpredicted C32 tRNA thiolase
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base T (3/3);  total (23/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 = 9.94e-01

GTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAG  >  W3110S.gb/1413207‑1413343
                                                                  |                                                                      
gTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgac                                                                    >  1:3033450/1‑71 (MQ=255)
 tCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgacg                                                                   <  1:1798715/71‑1 (MQ=255)
       tGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgacg                                                                   <  1:1928264/65‑1 (MQ=255)
         gTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCgag                                                           <  1:1600375/71‑1 (MQ=255)
         gTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCgag                                                           <  1:945838/71‑1 (MQ=255)
            aCCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGaa                                                         >  1:70598/1‑70 (MQ=255)
              ccAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAacac                                                      <  1:4484902/71‑1 (MQ=255)
              ccAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAaca                                                       >  1:3202158/1‑70 (MQ=255)
                 aCGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGt                                                   >  1:2279405/1‑71 (MQ=255)
                  cgcgATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTg                                                  <  1:543059/71‑1 (MQ=255)
                  cgcgATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTg                                                  <  1:3066636/71‑1 (MQ=255)
                  cgcgATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTg                                                  <  1:1515250/71‑1 (MQ=255)
                       tCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAacac                                                      >  1:2358137/1‑62 (MQ=255)
                         ttCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAg                                                    >  1:1551961/1‑62 (MQ=255)
                         ttCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGt                                           >  1:4208251/1‑71 (MQ=255)
                           cgtcgCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGc                                               <  1:1340863/65‑1 (MQ=255)
                           cgtcgCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGc                                               <  1:2225862/65‑1 (MQ=255)
                               gCCCCCAGTTCCGTTGCTGTACGATAAAGGATACCGCGACGAAGGCGAGaa                                                         <  1:2475242/51‑1 (MQ=255)
                                cccccAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTg                                    <  1:3354864/71‑1 (MQ=255)
                                   ccAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGtttt                                     <  1:296429/67‑1 (MQ=255)
                                    cAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCct                                >  1:1162683/1‑71 (MQ=255)
                                     aGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCctc                               >  1:2279001/1‑71 (MQ=255)
                                     aGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCct                                >  1:1772848/1‑70 (MQ=255)
                                     aGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAACTGGTTTTGCCctc                               >  1:2272187/1‑71 (MQ=255)
                                         ccGTTGCGGTACGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTgg                                         >  1:4549506/1‑57 (MQ=255)
                                         ccGTTGCGGTACGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTgg                                         >  1:320398/1‑57 (MQ=255)
                                         ccGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTgg                            >  1:329467/1‑70 (MQ=255)
                                            ttGCGGTACGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTgg                                         <  1:647407/54‑1 (MQ=255)
                                            ttGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGt                                        <  1:3076121/55‑1 (MQ=255)
                                            ttGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGt                                        <  1:1614306/55‑1 (MQ=255)
                                            ttGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGt                                        <  1:3960098/55‑1 (MQ=255)
                                            ttGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGaa                          <  1:2386349/69‑1 (MQ=255)
                                              gCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCct                                >  1:3195907/1‑61 (MQ=255)
                                                 gTACGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGCaa                                             >  1:570091/1‑45 (MQ=255)
                                                 gTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGc                                               <  1:271462/43‑1 (MQ=255)
                                                   aCGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtctt                      >  1:953561/1‑66 (MQ=255)
                                                   aCGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTGTCtt                 >  1:707637/1‑71 (MQ=255)
                                                    cGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCctct                              >  1:2790963/1‑57 (MQ=255)
                                                    cGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTgg                            <  1:4625787/59‑1 (MQ=255)
                                                     gATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTg                                    <  1:3772953/50‑1 (MQ=255)
                                                      aTAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtctt                      >  1:236077/1‑63 (MQ=255)
                                                        aaaGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCAc             >  1:886894/1‑70 (MQ=255)
                                                         aaGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGa           <  1:2667686/71‑1 (MQ=255)
                                                         aaGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGt                                           <  1:4181378/39‑1 (MQ=255)
                                                                ccGCTACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGt     <  1:1891809/70‑1 (MQ=255)
                                                                ccGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAt                         >  1:1118840/1‑50 (MQ=255)
                                                                ccGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGt     >  1:1976609/1‑70 (MQ=255)
                                                                  gcgACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAg  >  1:2981771/1‑71 (MQ=255)
                                                                  |                                                                      
GTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAG  >  W3110S.gb/1413207‑1413343

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: