Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,702,821 A→G 29.8% intergenic (‑102/+170) rnc ← / ← lepB RNase III/leader peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,702,8210AG29.8% 62.5 / 35.5 47intergenic (‑102/+170)rnc/lepBRNase III/leader peptidase
Reads supporting (aligned to +/- strand):  ref base A (17/16);  new base G (8/6);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAAT  >  W3110S.gb/2702751‑2702887
                                                                      |                                                                  
cAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtata                                                                    >  1:2652385/1‑71 (MQ=255)
    acacGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATggg                                                                <  1:3210263/71‑1 (MQ=255)
     cacGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATggg                                                                >  1:3826318/1‑70 (MQ=255)
      acGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGa                                                              >  1:19888/1‑71 (MQ=255)
       cGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATGTGGGGAc                                                             >  1:120067/1‑71 (MQ=255)
                 cTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATGTggg                                                                <  1:1234197/58‑1 (MQ=255)
                   gATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATGTGGGGa                                                              >  1:1126656/1‑58 (MQ=255)
                       taCATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCgtgt                                             >  1:2684904/1‑71 (MQ=255)
                       taCATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCgtgt                                             >  1:401977/1‑71 (MQ=255)
                           tGGTTCTGCTAACTGCTTCGCTACAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCg                                         <  1:3007267/71‑1 (MQ=255)
                            ggTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATGTGGGGACCAACGCTGTTTCGTGTGCCGt                                        >  1:1431326/1‑71 (MQ=255)
                             gTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATGTGGGGa                                                              >  1:4087169/1‑48 (MQ=255)
                              ttCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCgtgt                                             <  1:674453/64‑1 (MQ=255)
                              ttCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTgg                                      >  1:3412103/1‑71 (MQ=255)
                              ttCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTg                                       >  1:3410470/1‑70 (MQ=255)
                                cTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtata                                                                    >  1:3870222/1‑39 (MQ=255)
                                cTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGtata                                                                    >  1:3864909/1‑39 (MQ=255)
                                  gCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGa                                                              >  1:659473/1‑43 (MQ=255)
                                   cTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATGTGGGGACCAACGCTGTTTCGTGTGCCGTGGCa                                    <  1:2180566/68‑1 (MQ=255)
                                   cTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATGTGGGGACCAACGCTGTTTCGTGTGCCGTGGCa                                    <  1:4631246/68‑1 (MQ=255)
                                   cTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAAc                                  <  1:3263725/70‑1 (MQ=255)
                                    tAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGAc                                                             >  1:2130889/1‑42 (MQ=255)
                                     aaCTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGc                                                       >  1:4234329/1‑47 (MQ=255)
                                     aaCTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGc                                                       >  1:3669142/1‑47 (MQ=255)
                                      aCTGCTTCGCTGCAGTTTCAGTTTACCGGTATGTGGGGACCAACGCTg                                                     >  1:3743457/1‑48 (MQ=255)
                                      aCTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCa                                    <  1:3344772/65‑1 (MQ=255)
                                      aCTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTgg                              <  1:19579/71‑1 (MQ=255)
                                         gCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCaa                                                          >  1:2122629/1‑40 (MQ=255)
                                          cTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGt                          <  1:3154774/71‑1 (MQ=255)
                                           ttCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCa                                    <  1:2462217/60‑1 (MQ=255)
                                           ttCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTg                         >  1:1938250/1‑71 (MQ=255)
                                             cGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACg                                                        >  1:1593209/1‑38 (MQ=255)
                                              gCTGCAGTTTCAGTTTACCGGTATGTGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGcc                       <  1:1928866/70‑1 (MQ=255)
                                              gCTGCAGTTTCAGTTTACCGGTATGTGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGcc                       <  1:4283138/70‑1 (MQ=255)
                                               cTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCaa                     <  1:1200201/71‑1 (MQ=255)
                                                  cagtttcagtttACCGGTATGTGGGGACCACCGCTg                                                     <  1:4191401/36‑1 (MQ=37)
                                                    gtttcagtttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTa                <  1:3874543/71‑1 (MQ=255)
                                                    gtttcagtttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTa                <  1:1668151/71‑1 (MQ=255)
                                                     tttcagtttACCGGTATGTGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTa                >  1:1226574/1‑70 (MQ=255)
                                                     tttcagtttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAAc                                  >  1:3770044/1‑52 (MQ=255)
                                                       tcagtttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAAc             >  1:459094/1‑71 (MQ=255)
                                                         agtttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCa                                    <  1:3892941/46‑1 (MQ=255)
                                                           tttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGt                                        <  1:3760177/40‑1 (MQ=255)
                                                           tttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGt                                        <  1:418210/40‑1 (MQ=255)
                                                           tttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTgg                              <  1:4628024/50‑1 (MQ=255)
                                                            ttACCGGTATGTGGGGACCAACGCTGTTTCGTGTGCCGTGGCAAc                                  >  1:3324606/1‑45 (MQ=255)
                                                            ttACCGGTATGTGGGGACCAACGCTGTTTCGTGTGCCGTGGCAAc                                  >  1:1991735/1‑45 (MQ=255)
                                                             tACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAAc                                  >  1:1775020/1‑44 (MQ=255)
                                                             tACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGa                             >  1:90305/1‑49 (MQ=255)
                                                                  gTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAat  <  1:2623887/71‑1 (MQ=255)
                                                                   tatgTGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAat  >  1:4618800/1‑70 (MQ=255)
                                                                      |                                                                  
CAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAAT  >  W3110S.gb/2702751‑2702887

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: