Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,560,351 G→A 19.1% P62L (CCC→CTC)  rbn ← tRNA processing exoribonuclease BN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,560,3510GA19.1% 153.8 / 28.4 68P62L (CCC→CTC) rbntRNA processing exoribonuclease BN
Reads supporting (aligned to +/- strand):  ref base G (34/21);  new base A (8/5);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

TAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAG  >  W3110S.gb/3560289‑3560421
                                                              |                                                                      
tAGCAGGCAGAAAGTTGTCAAAAATAAAGTGACGTAACTGAATGCTGATGTCGGAAAACATgaga                                                                      <  1:2604191/65‑1 (MQ=255)
tAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATgaga                                                                      <  1:522803/65‑1 (MQ=255)
 aGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAgcgg                                                               >  1:3264704/1‑71 (MQ=255)
  gcaggcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGaaa                                                                    >  1:2593881/1‑65 (MQ=255)
  gcaggcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGaaa                                                                    >  1:2473033/1‑65 (MQ=255)
  gcaggcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGaaa                                                                    >  1:3891415/1‑65 (MQ=255)
      gcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGaaa                                                          <  1:4056210/71‑1 (MQ=255)
        agAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAg                                                        >  1:1012641/1‑71 (MQ=255)
        agAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAg                                                        >  1:2207902/1‑71 (MQ=255)
        agAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAg                                                        >  1:3492996/1‑71 (MQ=255)
         gAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAgc                                                                 <  1:160778/61‑1 (MQ=255)
          aaaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAgcg                                                      >  1:3630349/1‑71 (MQ=255)
           aaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcg                                                      >  1:4001572/1‑70 (MQ=255)
           aaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcg                                                      >  1:1693684/1‑70 (MQ=255)
           aaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcg                                                      >  1:559211/1‑70 (MQ=255)
           aaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcg                                                      >  1:2400676/1‑70 (MQ=255)
           aaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcg                                                      >  1:2849235/1‑70 (MQ=255)
            aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGaaaa                                                         >  1:2885621/1‑66 (MQ=255)
                ggCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCaaaaa                                                >  1:3621216/1‑71 (MQ=255)
                ggCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCaaaaa                                                >  1:3046442/1‑71 (MQ=255)
                 gCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                               <  1:302412/71‑1 (MQ=255)
                 gCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAaac                                               <  1:425503/71‑1 (MQ=255)
                 gCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAaac                                               <  1:3259071/71‑1 (MQ=255)
                  cAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaaca                                              <  1:180136/71‑1 (MQ=255)
                    aaaaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGa                                                                      >  1:2028173/1‑45 (MQ=255)
                      aaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACgg                                          >  1:482372/1‑71 (MQ=255)
                      aaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACgg                                          >  1:982836/1‑71 (MQ=255)
                        tAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGata                                                                    <  1:3138353/43‑3 (MQ=255)
                        tAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGaaa                                                                    <  1:1882950/43‑1 (MQ=255)
                        tAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACTACGGCa                                        <  1:522409/71‑1 (MQ=255)
                        tAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCa                                        <  1:2379224/71‑1 (MQ=255)
                                cGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                                <  1:3850866/71‑1 (MQ=255)
                                 gTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                                <  1:254481/70‑1 (MQ=255)
                                    aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                                >  1:2497653/1‑67 (MQ=255)
                                    aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCAc                             >  1:3554827/1‑70 (MQ=255)
                                    aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACt                            <  1:3715504/71‑1 (MQ=255)
                                    aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACt                            <  1:2538621/71‑1 (MQ=255)
                                    aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACt                            <  1:2012246/71‑1 (MQ=255)
                                    aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACt                            <  1:2455042/71‑1 (MQ=255)
                                      tGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTaa                          >  1:2402478/1‑71 (MQ=255)
                                      tGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTaa                          >  1:2669681/1‑71 (MQ=255)
                                      tGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTaa                          >  1:3275583/1‑71 (MQ=255)
                                      tGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaaa                                                  <  1:47424/47‑1 (MQ=255)
                                        aaTGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAAt                         <  1:2682182/70‑1 (MQ=255)
                                        aaTGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaaaaa                                                >  1:2414019/1‑47 (MQ=255)
                                        aaTGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaaaaa                                                >  1:302856/1‑47 (MQ=255)
                                        aaTGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaaaaa                                                >  1:908235/1‑47 (MQ=255)
                                         aTGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAAt                         <  1:3349648/69‑1 (MQ=255)
                                            cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAg                                                        >  1:3483050/1‑35 (MQ=255)
                                            cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                    >  1:3052717/1‑71 (MQ=255)
                                            cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                    >  1:11051/1‑71 (MQ=255)
                                            cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                    >  1:1982288/1‑71 (MQ=255)
                                            cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                    >  1:474472/1‑71 (MQ=255)
                                            cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                    >  1:2826090/1‑71 (MQ=255)
                                            cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                    >  1:2807250/1‑71 (MQ=255)
                                               aCGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCaa                 >  1:780758/1‑71 (MQ=255)
                                               aCGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGc                   <  1:2469169/69‑1 (MQ=255)
                                                   cGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGa              >  1:3607616/1‑70 (MQ=255)
                                                    ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            >  1:919524/1‑71 (MQ=255)
                                                    ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            >  1:1125429/1‑71 (MQ=255)
                                                    ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            >  1:118355/1‑71 (MQ=255)
                                                    ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca            <  1:1494379/71‑1 (MQ=255)
                                                    ggAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCCGCACTAATGAGAGCAACGaca            <  1:3167326/71‑1 (MQ=255)
                                                      aaaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca          <  1:1856709/71‑1 (MQ=255)
                                                       aaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACAt         <  1:2837863/71‑1 (MQ=255)
                                                       aaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACAt         <  1:1013905/71‑1 (MQ=255)
                                                          cATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAgg      >  1:1761749/1‑71 (MQ=255)
                                                          cATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAgg      >  1:3502915/1‑71 (MQ=255)
                                                          cATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAgg      >  1:3156470/1‑71 (MQ=255)
                                                           aTGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGc     <  1:2436363/71‑1 (MQ=255)
                                                             gggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca          <  1:2974573/64‑1 (MQ=255)
                                                              ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCaa   <  1:1890338/70‑1 (MQ=255)
                                                              ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAg  >  1:409993/1‑71 (MQ=255)
                                                              ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAg  >  1:2039734/1‑71 (MQ=255)
                                                              ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAg  >  1:2714678/1‑71 (MQ=255)
                                                              ggAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAg  >  1:3543342/1‑71 (MQ=255)
                                                              |                                                                      
TAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCAAG  >  W3110S.gb/3560289‑3560421

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: