Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,560,351 G→A 13.8% 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,3510GA13.8% 236.7 / 24.9 87P62L (CCC→CTC) rbntRNA processing exoribonuclease BN
Reads supporting (aligned to +/- strand):  ref base G (36/39);  new base A (6/6);  total (42/45)
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.96e-01

ATCGCCAGTAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCA  >  W3110S.gb/3560281‑3560419
                                                                      |                                                                    
aTCGCCAGTAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATgg                                                                      <  1:1360058/71‑1 (MQ=255)
   gCCAGTAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGaa                                                                   <  1:1203778/71‑1 (MQ=255)
    ccAGTAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGaaa                                                                  <  1:2088149/71‑1 (MQ=255)
       gTAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAgc                                                               >  1:4078748/1‑71 (MQ=255)
        tAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAgcg                                                              >  1:3024132/1‑71 (MQ=255)
         aGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAgcgg                                                             >  1:4296760/1‑71 (MQ=255)
          gcaggcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAgcggc                                                            >  1:3149211/1‑71 (MQ=255)
          gcaggcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAgcggc                                                            >  1:2359394/1‑71 (MQ=255)
           caggcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAgcggc                                                            >  1:3102227/1‑70 (MQ=255)
              gcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGaaa                                                        <  1:2891408/71‑1 (MQ=255)
              gcagAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGaa                                                         >  1:1915708/1‑70 (MQ=255)
                agAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAg                                                      <  1:415541/71‑1 (MQ=255)
                 gAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGaaaa                                                       <  1:1215783/69‑1 (MQ=255)
                 gAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGaaaa                                                       <  1:132981/69‑1 (MQ=255)
                 gAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgc                                                     <  1:46592/71‑1 (MQ=255)
                   aaGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAgcgc                                                   <  1:1842423/71‑1 (MQ=255)
                    aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcgc                                                   >  1:2148275/1‑70 (MQ=255)
                    aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcgc                                                   >  1:107694/1‑70 (MQ=255)
                    aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcgc                                                   >  1:1192809/1‑70 (MQ=255)
                    aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgcg                                                    >  1:1206680/1‑69 (MQ=255)
                    aGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCa                                                  <  1:2587143/71‑1 (MQ=255)
                     gTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCa                                                  >  1:4213495/1‑70 (MQ=255)
                       tagCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaaa                                                >  1:2031758/3‑70 (MQ=255)
                        ggCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAgc                                                     >  1:2597848/1‑64 (MQ=255)
                         gCAAAAATAAAGTGACTTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                             <  1:1841551/71‑1 (MQ=255)
                         gCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaaaaa                                              >  1:4266277/1‑70 (MQ=255)
                          cAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                             <  1:2316131/70‑1 (MQ=255)
                           aaaaaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaaca                                            >  1:3288034/1‑70 (MQ=255)
                           aaaaaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGCCGAAAAGCGCaa                                                 <  1:454232/65‑1 (MQ=255)
                           aaaaaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCaaa                                                <  1:2209085/66‑1 (MQ=255)
                            aaaaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaacaac                                          >  1:3230344/1‑71 (MQ=255)
                              aaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaacaa                                           >  1:2898763/1‑68 (MQ=255)
                              aaTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaacaa                                           >  1:2585452/1‑68 (MQ=255)
                               aTAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGc                                       <  1:993620/71‑1 (MQ=255)
                                tAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCa                                      <  1:2140627/71‑1 (MQ=255)
                                    gTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                             <  1:4149654/60‑1 (MQ=255)
                                    gTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAAc                                    >  1:199241/1‑69 (MQ=255)
                                    gTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCa                                  <  1:1897750/71‑1 (MQ=255)
                                    gTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCa                                  <  1:4522374/71‑1 (MQ=255)
                                      gACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaacaa                                           >  1:147705/1‑60 (MQ=255)
                                        cGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                              <  1:285704/71‑1 (MQ=255)
                                        cGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                              <  1:1340431/71‑1 (MQ=255)
                                         gTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                              >  1:923935/1‑70 (MQ=255)
                                         gTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                              >  1:600965/1‑70 (MQ=255)
                                          tAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCgg                              <  1:4520226/69‑1 (MQ=255)
                                           aaCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCAc                           >  1:2145732/1‑71 (MQ=255)
                                           aaCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCAc                           >  1:3590098/1‑71 (MQ=255)
                                            aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGc                             >  1:4111605/1‑68 (MQ=255)
                                            aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCAc                           >  1:2864967/1‑70 (MQ=255)
                                            aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACt                          <  1:243101/71‑1 (MQ=255)
                                            aCTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACt                          <  1:1735397/71‑1 (MQ=255)
                                            aCTGAATGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAgc                                                     >  1:756977/1‑44 (MQ=255)
                                             cTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                             <  1:2704829/51‑1 (MQ=255)
                                               gAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCa                                      <  1:1583788/56‑1 (MQ=255)
                                                aaTGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGTCACTAAt                       <  1:4468560/70‑1 (MQ=255)
                                                aaTGCTGACGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATg                      <  1:3645413/71‑1 (MQ=255)
                                                aaTGCTGACGTCGGAAAACATCGGAAAAGCTGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAAt                       >  1:876562/1‑70 (MQ=255)
                                                 aTGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCaa                                     >  1:1248821/1‑55 (MQ=255)
                                                 aTGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAAt                       >  1:686063/1‑69 (MQ=255)
                                                 aTGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAAt                       >  1:1934939/1‑69 (MQ=255)
                                                  tGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCaa                                     <  1:1019566/54‑1 (MQ=255)
                                                  tGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAAt                       >  1:647038/1‑68 (MQ=255)
                                                   ggtGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATga                     >  1:393725/3‑69 (MQ=255)
                                                    cTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGagag                  >  1:2952578/1‑71 (MQ=255)
                                                     tGACGTAGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGc                                       <  1:2114740/49‑1 (MQ=255)
                                                       aCGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGc                 <  1:4337281/69‑1 (MQ=255)
                                                       aCGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGc                 <  1:1536593/69‑1 (MQ=255)
                                                       aCGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCa                <  1:581376/70‑1 (MQ=255)
                                                       aCGTCGGAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCa                <  1:3034653/70‑1 (MQ=255)
                                                        cGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAAc              <  1:2453132/71‑1 (MQ=255)
                                                        cGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAAc              <  1:3088045/71‑1 (MQ=255)
                                                          tCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                             <  1:682442/38‑1 (MQ=255)
                                                          tCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAaac                                             <  1:4327433/38‑1 (MQ=255)
                                                            ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATga                     >  1:3925909/1‑60 (MQ=255)
                                                            ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATg                      >  1:1804562/1‑59 (MQ=255)
                                                            ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCaa               >  1:685232/1‑66 (MQ=255)
                                                            ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca          >  1:3754626/1‑71 (MQ=255)
                                                            ggAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGac           <  1:130633/70‑1 (MQ=255)
                                                            ggAAAACATGAGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGaca          <  1:2526232/71‑1 (MQ=255)
                                                              aaaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca        <  1:1891524/71‑1 (MQ=255)
                                                              aaaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGAcaca        <  1:2848794/71‑1 (MQ=255)
                                                               aaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAAc                                    <  1:3605002/42‑1 (MQ=255)
                                                               aaaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACAt       <  1:1775403/71‑1 (MQ=255)
                                                                aaCATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGc                 <  1:4050954/60‑1 (MQ=255)
                                                                  cATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTTATGAGAGCAACGACACATAgg    >  1:3676448/1‑71 (MQ=255)
                                                                  cATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAgg    >  1:3264713/1‑71 (MQ=255)
                                                                   aTGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCa                            >  1:1131085/1‑46 (MQ=255)
                                                                    tGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCa  <  1:2998344/71‑1 (MQ=255)
                                                                      |                                                                    
ATCGCCAGTAGCAGGCAGAAAGTTGGCAAAAATAAAGTGACGTAACTGAATGCTGACGTCGGAAAACATGGGAAAAGCGGCGAAAAGCGCAAAAACAACGGCAACCAGCGGCACTAATGAGAGCAACGACACATAGGCA  >  W3110S.gb/3560281‑3560419

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: