Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,449,319 A→G 35.7% G671G (GGT→GGC rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,449,3190AG35.7% 20.4 / 25.8 28G671G (GGT→GGCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (9/9);  new base G (5/5);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

GATACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTC  >  W3110S.gb/3449251‑3449385
                                                                    |                                                                  
gATACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTaccagac                                                                   <  1:3117552/70‑1 (MQ=255)
gATACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTaccagac                                                                   <  1:4029270/70‑1 (MQ=255)
gATACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTaccagac                                                                   <  1:2417672/70‑1 (MQ=255)
gATACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTaccagac                                                                   <  1:2369688/70‑1 (MQ=255)
 aTACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTaccaggcca                                                                 <  1:3509271/71‑1 (MQ=255)
 aTACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTaccaggcca                                                                 <  1:3891768/71‑1 (MQ=255)
 aTACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTaccagacca                                                                 <  1:1182213/71‑1 (MQ=255)
 aTACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTaccagacca                                                                 <  1:3616398/71‑1 (MQ=255)
    cacGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTAccagaccaga                                                               <  1:1751692/70‑1 (MQ=255)
       gATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTgg                                                           <  1:3953347/71‑1 (MQ=255)
       gATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTgg                                                           <  1:347502/71‑1 (MQ=255)
       gATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTgg                                                           <  1:3938052/71‑1 (MQ=255)
       gATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTgg                                                           <  1:1488937/71‑1 (MQ=255)
       gATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTgg                                                           <  1:2040320/71‑1 (MQ=255)
          cGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAAc                                                        >  1:1868180/1‑71 (MQ=255)
          cGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAAc                                                        >  1:3222532/1‑71 (MQ=255)
          cGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAAc                                                        >  1:3971353/1‑71 (MQ=255)
              cgccgcAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTAccagaccaga                                                               >  1:2386164/1‑60 (MQ=255)
              cgccgcAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCt                                                    >  1:3687437/1‑71 (MQ=255)
              cgccgcAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCt                                                    >  1:3674838/1‑71 (MQ=255)
                   cAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCTCCTGa                                               <  1:2585574/71‑1 (MQ=255)
                   cAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCTCCTGa                                               <  1:998081/71‑1 (MQ=255)
                   cAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCTCCTGa                                               <  1:3040611/71‑1 (MQ=255)
                   cAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCTCCTGa                                               <  1:1846461/71‑1 (MQ=255)
                            tatCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGc                                       >  1:3641740/1‑70 (MQ=255)
                            tatCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGc                                       >  1:587005/1‑70 (MQ=255)
                            tatCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCa                                      <  1:1633398/71‑1 (MQ=255)
                             atCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTAccagaccaga                                                               >  1:579943/1‑45 (MQ=255)
                             atCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCaa                                     >  1:1223363/1‑71 (MQ=255)
                                gATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCTCCTGAATTTCAGCAActt                                  >  1:3633140/1‑71 (MQ=255)
                                     cTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAActtctgct                             <  1:289192/71‑1 (MQ=255)
                                      tttGTTGTAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCTCCTGAATTTCAGCaa                                     >  1:1813902/1‑62 (MQ=255)
                                       ttgttgTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAActtctgcttc                           >  1:3637196/1‑71 (MQ=255)
                                        tgttgtAGCGTTCGCCCGCAGTTACCAGGCCAGACTGGAACTGCTCCTGAATTTCAGCAAc                                    >  1:3026628/1‑61 (MQ=255)
                                             tAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCt                      <  1:1593159/70‑1 (MQ=255)
                                                    cGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGga                  >  1:3548816/1‑67 (MQ=255)
                                                            gTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCg                    <  1:2997850/57‑1 (MQ=255)
                                                             ttACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCgtgt     <  1:2915905/71‑1 (MQ=255)
                                                              tACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAActtctgctt                            >  1:83047/1‑47 (MQ=255)
                                                               accaggccagACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGtt    <  1:2774611/70‑1 (MQ=255)
                                                               accaggccagACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGtt    <  1:1742794/70‑1 (MQ=255)
                                                               accaggacagACTGGAACTGCTCCTGAATTTCAGCAActt                                  >  1:1768190/1‑40 (MQ=38)
                                                               accagaccagACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGttt   >  1:245887/1‑71 (MQ=255)
                                                               accagaccagACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGttt   >  1:684726/1‑71 (MQ=255)
                                                               accagaccagACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGtt    <  1:2774310/70‑1 (MQ=255)
                                                               accagaccagACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGtt    <  1:574375/70‑1 (MQ=255)
                                                                ccagaccagaCTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTttc  <  1:2569899/71‑1 (MQ=255)
                                                                ccagaccagaCTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTttc  <  1:1605155/71‑1 (MQ=255)
                                                                    |                                                                  
GATACACGATCGTTCGCCGCAGCCCAGATATCGATAACTTTGTTGTAGCGTTCGCCCGCAGTTACCAGACCAGACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTC  >  W3110S.gb/3449251‑3449385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: