Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,449,061 C→T 5.3% T757T (ACG→ACA 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,0610CT 5.3% 430.2 / 8.8 131T757T (ACG→ACArpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (71/52);  new base T (2/5);  total (73/58)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 133 total aligned reads displayed.

ACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCG  >  W3110S.gb/3448992‑3449131
                                                                     |                                                                      
 cGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGt                                                                      >  1:6594588/1‑71 (MQ=255)
   aGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                    <  1:4489252/71‑1 (MQ=255)
   aGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                    <  1:211584/71‑1 (MQ=255)
    gCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTc                                                                   <  1:1559027/71‑1 (MQ=255)
       ccGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgat                                                                >  1:5696079/1‑71 (MQ=255)
         gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                              <  1:5618265/71‑1 (MQ=255)
         gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                              <  1:4805610/71‑1 (MQ=255)
           gggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGatgatg                                                            >  1:3400570/1‑71 (MQ=255)
            ggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGatgatg                                                            <  1:3059763/70‑1 (MQ=255)
              ttgagaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                    <  1:3354661/58‑1 (MQ=255)
                gagaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatg                                                               >  1:5402488/1‑63 (MQ=255)
                 agaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCa                                                      >  1:2817089/1‑71 (MQ=255)
                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCg                                                                       <  1:855146/53‑1 (MQ=255)
                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                     <  1:5055905/71‑1 (MQ=255)
                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                     <  1:3285446/71‑1 (MQ=255)
                     gaggTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCg                                                  >  1:3405396/4‑71 (MQ=255)
                     aacgTACTACAGTACTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCg                                                  >  1:4858885/4‑70 (MQ=255)
                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCg                                                  <  1:4834278/70‑1 (MQ=255)
                        gTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCt                                               >  1:2819198/1‑71 (MQ=255)
                        gTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCt                                               <  1:6761310/71‑1 (MQ=255)
                         tACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                    <  1:7781136/49‑1 (MQ=255)
                         tACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                    <  1:1727769/49‑1 (MQ=255)
                         tACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                    <  1:5851906/49‑1 (MQ=255)
                         tACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                    <  1:4901691/49‑1 (MQ=255)
                           cTGGAGTACGTTCGGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATcc                                                   <  1:4065550/64‑1 (MQ=255)
                           cTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                    <  1:170935/47‑1 (MQ=255)
                           cTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                  >  1:5703055/1‑49 (MQ=255)
                           cTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATcc                                                   <  1:7416417/64‑1 (MQ=255)
                              gAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGa                                                          <  1:5553744/54‑1 (MQ=255)
                                gTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCAt                                        >  1:7883882/1‑70 (MQ=255)
                                 tACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGc                                                                        <  1:2734696/37‑1 (MQ=255)
                                 tACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGTGTTTCGATGATGGAGCCATCCGGCTTCGCCAt                                        >  1:3988563/1‑69 (MQ=255)
                                 tACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTc                                             >  1:369240/1‑64 (MQ=255)
                                  aCGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCAt                                        >  1:1700237/1‑68 (MQ=255)
                                   cGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGTGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGa                                    >  1:7932693/1‑71 (MQ=255)
                                   cGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCg                                                                       >  1:2892166/1‑36 (MQ=255)
                                   cGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCg                                                                       >  1:6548703/1‑36 (MQ=255)
                                   cGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGtt                                                                     >  1:3161774/1‑38 (MQ=255)
                                   cGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGtt                                                                     >  1:5380157/1‑38 (MQ=255)
                                    gTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCAt                                        <  1:4116256/66‑1 (MQ=255)
                                    gTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGa                                    >  1:5278114/1‑70 (MQ=255)
                                     ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGc                                                        <  1:411106/49‑1 (MQ=255)
                                      tCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCa                                 >  1:5403461/1‑71 (MQ=255)
                                      tCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCa                                 >  1:2338884/1‑71 (MQ=255)
                                        aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                     <  1:6832132/49‑1 (MQ=255)
                                        aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                     <  1:7196915/49‑1 (MQ=255)
                                        aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATcc                                                   >  1:4764928/1‑51 (MQ=255)
                                        aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGa                                    >  1:2419108/1‑66 (MQ=255)
                                         gACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATcc                                                   >  1:2851704/1‑50 (MQ=255)
                                         gACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATcc                                                   >  1:7330932/1‑50 (MQ=255)
                                         gACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCtt                                              >  1:4311700/1‑55 (MQ=255)
                                         gACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCt                                               >  1:1631124/1‑54 (MQ=255)
                                         gACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                                >  1:321545/1‑69 (MQ=255)
                                         gACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                                >  1:1590097/1‑69 (MQ=255)
                                           ccTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCAt                            <  1:1168661/71‑1 (MQ=255)
                                            cTTCACGGAAGTTCGCGGTGATTGGAGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACg                               <  1:5288200/67‑1 (MQ=255)
                                              tCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc                          >  1:6729138/1‑70 (MQ=255)
                                              tCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc                          >  1:180041/1‑70 (MQ=255)
                                               cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGTGCCATCCGGCTTCGCCATCa                                      >  1:5533281/1‑57 (MQ=255)
                                               cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCAtaa                                      >  1:6946680/1‑55 (MQ=255)
                                               cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGa                                    >  1:1501321/1‑59 (MQ=255)
                                               cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc                          >  1:4804969/1‑69 (MQ=255)
                                                aCGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAg                                                         >  1:2069629/1‑37 (MQ=255)
                                                aCGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAg                       >  1:318445/1‑71 (MQ=255)
                                                aCGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAg                       >  1:5846720/1‑71 (MQ=255)
                                                  ggAAGTTCGCGGTGATTGGTGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCa                     <  1:4510510/71‑1 (MQ=255)
                                                  ggAAGTTCGCGGTGATTGGTGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCa                     <  1:4509588/71‑1 (MQ=255)
                                                  ggAAGTTCGCGGTGATTGGTGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCa                     <  1:6199303/71‑1 (MQ=255)
                                                  ggAAGTTCGCGGTGATTGGTGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCa                     <  1:3413741/71‑1 (MQ=255)
                                                  ggAAGTTCGCGGTGATTGGTGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCa                     <  1:6716236/71‑1 (MQ=255)
                                                  ggAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                      >  1:7446772/1‑70 (MQ=255)
                                                   gAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc                          >  1:3106343/1‑65 (MQ=255)
                                                   gAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCaa                    <  1:3895425/71‑1 (MQ=255)
                                                    aaGTTCGCGGTGATTTGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGa                                    >  1:7307000/1‑54 (MQ=255)
                                                     aGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                                                 >  1:5827682/1‑40 (MQ=255)
                                                     aGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGa                                    >  1:2690202/1‑53 (MQ=255)
                                                      gTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt                 <  1:820396/71‑1 (MQ=255)
                                                       ttCGCGGTTATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTg                <  1:6211507/71‑1 (MQ=255)
                                                       ttCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTg                <  1:5113429/71‑1 (MQ=255)
                                                        tCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGa               >  1:2345178/1‑71 (MQ=255)
                                                            ggTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCAt                            <  1:7729207/54‑1 (MQ=255)
                                                            ggTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGaa           >  1:1923901/1‑71 (MQ=255)
                                                            ggTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGAATACCAGCAAGCTGACGa            >  1:3294459/1‑70 (MQ=255)
                                                             gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt          <  1:6571312/71‑1 (MQ=255)
                                                               gATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa                             >  1:2629437/1‑50 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCAt                            <  1:2356105/50‑1 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAATCTg                >  1:3048797/1‑62 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTg                >  1:2954820/1‑62 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTg                >  1:3944516/1‑62 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCtg       <  1:3044854/71‑1 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCtg       >  1:738776/1‑71 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCtg       >  1:5667763/1‑71 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCtg       <  1:5144955/71‑1 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt        <  1:2534179/70‑1 (MQ=255)
                                                                aTTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt        <  1:5184609/70‑1 (MQ=255)
                                                                 ttGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCa                                 >  1:2430578/1‑44 (MQ=255)
                                                                 ttGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCa                                 >  1:6190551/1‑44 (MQ=255)
                                                                 ttGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCa                                 >  1:4194275/1‑44 (MQ=255)
                                                                 ttGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCtgt      >  1:6941356/1‑71 (MQ=255)
                                                                  tGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                      <  1:2107099/54‑1 (MQ=255)
                                                                     |                                                                      
ACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCG  >  W3110S.gb/3448992‑3449131

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: