Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,413,523 G→A 37.8% S47F (TCC→TTC)  ydaO ← predicted C32 tRNA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,413,5230GA37.8% 38.2 / 46.9 45S47F (TCC→TTC) ydaOpredicted C32 tRNA thiolase
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base A (10/7);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTAC  >  W3110S.gb/1413459‑1413592
                                                                |                                                                     
gattgGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAgg                                                                  >  1:4239867/1‑70 (MQ=255)
gattgGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAgg                                                                  >  1:517119/1‑70 (MQ=255)
  ttgGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCaa                                                               <  1:892039/71‑1 (MQ=255)
  ttgGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCaa                                                               <  1:885723/71‑1 (MQ=255)
        cgcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCat                                                          <  1:4881568/70‑1 (MQ=255)
         gcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCatga                                                        >  1:988988/1‑71 (MQ=255)
         gcTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCatga                                                        >  1:4845393/1‑71 (MQ=255)
                     aTTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCct                                             <  1:2700512/70‑1 (MQ=255)
                       tGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCcttct                                          >  1:1827795/1‑71 (MQ=255)
                       tGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCcttct                                          >  1:1953303/1‑71 (MQ=255)
                           cAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCaa                                      <  1:4016047/71‑1 (MQ=255)
                             gAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCg                                                    <  1:371340/55‑1 (MQ=255)
                             gAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCg                                                    <  1:2857228/55‑1 (MQ=255)
                              aaTCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCa                                   <  1:4184269/71‑1 (MQ=255)
                              aaTCTCCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCa                                   >  1:4757464/1‑71 (MQ=255)
                                tctcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAgg                                                                  >  1:982040/1‑38 (MQ=255)
                                tctcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCata                                 >  1:736274/1‑71 (MQ=255)
                                tctcCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGTTCGCCTTCTTCAATCata                                 >  1:1915337/1‑71 (MQ=255)
                                  tcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCg                                                    <  1:2672865/50‑1 (MQ=255)
                                  tcCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCg                                                    <  1:1718322/50‑1 (MQ=255)
                                  tcCAGCATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATAtt                               >  1:3989526/1‑71 (MQ=255)
                                       cATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGt                          <  1:4136052/71‑1 (MQ=255)
                                       cATGGTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCg                                                    <  1:1304680/45‑1 (MQ=255)
                                           gTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGc                      <  1:2649123/71‑1 (MQ=255)
                                           gTATAGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGc                      <  1:4119267/71‑1 (MQ=255)
                                            tataGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGc                      >  1:1121193/1‑70 (MQ=255)
                                            tataGCTGTCTTTACCCCCGAAGAGGCAAACCATGATGc                                                     <  1:3562744/39‑1 (MQ=255)
                                              taGCTGTCTTTACCCCCGGAGAGGCAAACCATAATGCg                                                    <  1:89605/38‑1 (MQ=255)
                                                gCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAAt                                     >  1:3363542/1‑51 (MQ=255)
                                                gCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCa                     >  1:1117111/1‑67 (MQ=255)
                                                gCTGTCTTTACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCa                     >  1:3338377/1‑67 (MQ=255)
                                                  tGTCTTTACCCCAGGAGAGGCAAACCATGATGCGAt                                                  >  1:892057/1‑36 (MQ=255)
                                                   gTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCat                                  <  1:4838586/51‑1 (MQ=255)
                                                      tttACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGc           >  1:311132/1‑71 (MQ=255)
                                                      tttACCCCCGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGc           >  1:1329488/1‑71 (MQ=255)
                                                       ttACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGcc          >  1:4412424/1‑71 (MQ=255)
                                                       ttACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGcc          >  1:1846085/1‑71 (MQ=255)
                                                          cccccGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGaa                            >  1:4355774/1‑50 (MQ=255)
                                                           ccccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAg                       <  1:3550049/54‑1 (MQ=255)
                                                            cccGAAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGt     <  1:2161468/71‑1 (MQ=255)
                                                             ccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGc                      <  1:4526284/53‑1 (MQ=255)
                                                             ccGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCa                     >  1:1251722/1‑54 (MQ=255)
                                                              cGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGaa                            >  1:2609099/1‑46 (MQ=255)
                                                              cGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGaa                            >  1:366397/1‑46 (MQ=255)
                                                              cGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTa   >  1:2896542/1‑71 (MQ=255)
                                                               ggAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCaa                    >  1:3510449/1‑53 (MQ=255)
                                                                gagaGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTac  >  1:1956180/1‑70 (MQ=255)
                                                                |                                                                     
GATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATGATGCGATCGCCTTCTTCAATCATATTGAAGTCAGCAATGGCTTCGCCCACGTTAC  >  W3110S.gb/1413459‑1413592

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: