Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,803,912 G→T 18.3% H125N (CAT→AAT)  thrS ← threonyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,803,9120GT18.3% 141.4 / 12.5 49H125N (CAT→AAT) thrSthreonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base T (5/4);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AGTTTCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAA  >  W3110S.gb/1803844‑1803980
                                                                    |                                                                    
aGTTTCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCa                                                                    >  1:1331867/1‑71 (MQ=255)
aGTTTCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCa                                                                    >  1:682676/1‑71 (MQ=255)
aGTTTCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCa                                                                    >  1:4180337/1‑71 (MQ=255)
aGTTTCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCa                                                                    <  1:2246958/71‑1 (MQ=255)
aGTTTCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCa                                                                    >  1:2246826/1‑71 (MQ=255)
aGTTTCACGCGCTTCGTGCCAGCTGACTTTTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCa                                                                    >  1:553741/1‑70 (MQ=255)
  tttCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATc                                                                  >  1:321243/1‑71 (MQ=255)
    tCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATcc                                                                 <  1:2141857/70‑1 (MQ=255)
      aCGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGc                                                               <  1:1571826/70‑1 (MQ=255)
        gcgcTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTtc                                                            <  1:509741/71‑1 (MQ=255)
          gcTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATCCGCTtctc                                                          <  1:144276/71‑1 (MQ=255)
                  ccAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCtt                                                  <  1:4067506/71‑1 (MQ=255)
                     gCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGa                                               <  1:639469/71‑1 (MQ=255)
                        gACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACAt                                            >  1:4046100/1‑71 (MQ=255)
                           tttCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTc                                                 <  1:4168538/63‑1 (MQ=255)
                            ttcttcTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCt                                                              <  1:2298911/49‑1 (MQ=255)
                             tcttctTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGa                                               >  1:4517498/1‑63 (MQ=255)
                              cttcttAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATc                                                                  >  1:318803/1‑43 (MQ=255)
                               ttcttAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTg                                     >  1:211989/1‑71 (MQ=255)
                               ttcttAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTg                                     >  1:989255/1‑71 (MQ=255)
                                tcttAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTtctc                                                          <  1:1464693/49‑1 (MQ=255)
                                   tAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGa                                               >  1:436137/1‑57 (MQ=255)
                                     aTGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTaa                               <  1:3553073/71‑1 (MQ=255)
                                           tcgtAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTa                                <  1:855545/64‑1 (MQ=255)
                                            cgtAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                                       >  1:747592/1‑40 (MQ=255)
                                             gtAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGtt                                 <  1:2067621/61‑1 (MQ=255)
                                             gtAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                                       >  1:4607022/1‑39 (MQ=255)
                                             gtAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTaa                               >  1:4176329/1‑63 (MQ=255)
                                             gtAGTTTTTCTCAGCAAGCTCATGCATCAGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGt                       >  1:2128649/1‑71 (MQ=255)
                                                gTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGAc                                              <  1:2766693/45‑1 (MQ=255)
                                                 tttttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACg                             <  1:2747707/61‑1 (MQ=255)
                                                  ttttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTc                                                 <  1:735694/40‑1 (MQ=255)
                                                   tttCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGc                                                    >  1:3061030/1‑36 (MQ=255)
                                                   tttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTc                      >  1:1505470/1‑66 (MQ=255)
                                                   tttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGa                  <  1:2922245/70‑1 (MQ=255)
                                                    ttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTc                                        >  1:4235890/1‑47 (MQ=255)
                                                     tctcAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATc                >  1:1756448/1‑70 (MQ=255)
                                                     tctcAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATc                <  1:3807153/70‑1 (MQ=255)
                                                         aGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAAcgt           <  1:4036902/71‑1 (MQ=255)
                                                         aGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCaa              <  1:4575031/68‑1 (MQ=255)
                                                           cAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCTGTCAAGATCAAc             <  1:2323417/67‑1 (MQ=255)
                                                            aaGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTc                                        >  1:4116110/1‑39 (MQ=255)
                                                             aGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGtt                                 <  1:2928442/45‑1 (MQ=255)
                                                              gCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGtaa      <  1:3324213/71‑1 (MQ=255)
                                                                tCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGtt                                 <  1:2744236/42‑1 (MQ=255)
                                                                tCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGtaat     >  1:3550249/1‑70 (MQ=255)
                                                                tCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGtaa      >  1:3644525/1‑69 (MQ=255)
                                                                  aTGCATCCGCTTCTAGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATaaa  <  1:3733765/71‑1 (MQ=255)
                                                                   tGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTg                                     <  1:3861454/35‑1 (MQ=255)
                                                                    |                                                                    
AGTTTCACGCGCTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATCAACGTCGTAATAAA  >  W3110S.gb/1803844‑1803980

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: