Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 940,557 A→G 78.0% P236P (CCA→CCG serS → seryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb940,5570AG78.0% 60.9 / 13.7 32P236P (CCA→CCGserSseryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base G (19/5);  total (22/10)
Fisher's exact test for biased strand distribution p-value = 2.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

CGATCTGTTCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGCC  >  W3110S.gb/940494‑940619
                                                               |                                                              
cGATCTGTTCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAg                                                         <  1:1470492/71‑1 (MQ=255)
    cTGTTCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGt                                                     >  1:2018149/1‑71 (MQ=255)
       ttCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTcc                                                  <  1:2150397/71‑1 (MQ=255)
         ccATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGc                                                >  1:1628221/1‑71 (MQ=255)
         ccATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTCAGc                                                >  1:2606906/1‑71 (MQ=255)
          cATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCt                                               <  1:1843902/71‑1 (MQ=255)
              cTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGAc                                            >  1:77706/1‑70 (MQ=255)
               tCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACt                                           >  1:941722/1‑70 (MQ=255)
                    cGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTg                                     <  1:718874/71‑1 (MQ=255)
                              agaagCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGt                            >  1:830043/1‑70 (MQ=255)
                               gaagCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAAc                                        >  1:1080647/1‑57 (MQ=255)
                                 agCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGaaa                        >  1:675122/1‑71 (MQ=255)
                                      acacCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGa                          <  1:995024/64‑1 (MQ=255)
                                      acacCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGa                          <  1:863501/64‑1 (MQ=255)
                                      acacCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGa                          <  1:1924683/64‑1 (MQ=255)
                                        acCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTc                                                   >  1:2112116/1‑37 (MQ=255)
                                        acCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTc                                                   >  1:34455/1‑37 (MQ=255)
                                        acCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTc                                                   >  1:452212/1‑37 (MQ=255)
                                         ccAGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCTCGGTGaaa                        <  1:1243035/63‑1 (MQ=255)
                                           aGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGaa                         >  1:119661/1‑60 (MQ=255)
                                           aGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAatc                      >  1:2333513/1‑63 (MQ=255)
                                           aGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATg               >  1:102639/1‑70 (MQ=255)
                                           aGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATg               >  1:810721/1‑70 (MQ=255)
                                            gTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGAGCGCGGTGAAATCATCGATGaa             <  1:413484/71‑1 (MQ=255)
                                               aCTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCt                                      >  1:2092595/1‑43 (MQ=255)
                                               aCTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAgat          >  1:1476963/1‑71 (MQ=255)
                                                 tATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAAc                                        >  1:785214/1‑39 (MQ=255)
                                                     cgcTGATCCCTACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTg    >  1:179723/1‑71 (MQ=255)
                                                     cgcTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAat                       >  1:2210001/1‑52 (MQ=255)
                                                     cgcTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTg    >  1:157955/1‑71 (MQ=255)
                                                        tGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGcc  >  1:2179951/1‑70 (MQ=255)
                                                          aTCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGa                          <  1:1762252/44‑1 (MQ=255)
                                                               |                                                              
CGATCTGTTCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGCC  >  W3110S.gb/940494‑940619

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: