Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,254,746 A→C 12.2% R764R (CGT→CGG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,254,7460AC12.2% 188.6 / 24.8 82R764R (CGT→CGGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (27/45);  new base C (6/4);  total (33/49)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

TATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAC  >  Exported/2254682‑2254813
                                                                |                                                                   
tATCCGCCAGACCTTTACGAGCACCGTGGGTGGAAATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAg                                                               >  1:4768096/1‑71 (MQ=255)
tATCCGCCAGACCTTTACGAGCACCGTCGGTGGAAATGAAGTACTGGAGTACGTTCAGACCTTCACGGaa                                                                >  1:3912804/1‑70 (MQ=255)
tATCCGCCAGACCTTTACGAGCACCGAGGGTGGAAATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAg                                                               >  1:4922957/1‑71 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTGCGTTCAGACCTTCACGGAAGtt                                                             >  1:4246128/1‑69 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGaa                                                                >  1:1829973/1‑66 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGtt                                                             >  1:18555/1‑69 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGtt                                                             >  1:3861999/1‑69 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                            <  1:712658/70‑1 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                            <  1:1835123/70‑1 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                            <  1:17428/70‑1 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                            <  1:1205087/70‑1 (MQ=255)
    cGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                            <  1:1048281/70‑1 (MQ=255)
     gCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTc                                                            <  1:3495751/69‑1 (MQ=255)
         gACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCa                                                                     >  1:4848455/1‑56 (MQ=255)
         gACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGt                                                       >  1:3964406/1‑70 (MQ=255)
         gACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGt                                                       >  1:1590413/1‑70 (MQ=255)
            cTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                              <  1:4914160/60‑1 (MQ=255)
            cTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                              <  1:3699838/60‑1 (MQ=255)
            cTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                              <  1:482234/60‑1 (MQ=255)
            cTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                              <  1:1133119/60‑1 (MQ=255)
                 cGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGt                                              <  1:3602964/71‑1 (MQ=255)
                     cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                          <  1:3680419/71‑1 (MQ=255)
                     cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                          <  1:3674579/71‑1 (MQ=255)
                     cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                          <  1:1036823/71‑1 (MQ=255)
                     cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                          <  1:866263/71‑1 (MQ=255)
                     cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                          <  1:2490041/71‑1 (MQ=255)
                       ccGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgat                                        <  1:3583096/71‑1 (MQ=255)
                         gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                      >  1:1140902/1‑71 (MQ=255)
                         gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                      >  1:2753853/1‑71 (MQ=255)
                         gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                      >  1:4459948/1‑71 (MQ=255)
                         gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                      >  1:2141417/1‑71 (MQ=255)
                         gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                      <  1:2328737/71‑1 (MQ=255)
                          tgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatgat                                     >  1:4466266/1‑71 (MQ=255)
                          tgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatgat                                     >  1:3934677/1‑71 (MQ=255)
                            ggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGatgatg                                    <  1:250970/70‑1 (MQ=255)
                              tggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGt                                                              <  1:5082925/42‑1 (MQ=255)
                                gagaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGc                                >  1:1973328/1‑70 (MQ=255)
                                 agaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCa                              <  1:3638102/71‑1 (MQ=255)
                                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTcg                                                           <  1:4113205/41‑1 (MQ=255)
                                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTcg                                                           <  1:1291735/41‑1 (MQ=255)
                                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTcg                                                           <  1:504150/41‑1 (MQ=255)
                                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTcg                                                           <  1:3058502/41‑1 (MQ=255)
                                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                             <  1:3370327/71‑1 (MQ=255)
                                  gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                             <  1:1619030/71‑1 (MQ=255)
                                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                            <  1:4338288/52‑1 (MQ=255)
                                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                            <  1:4158060/52‑1 (MQ=255)
                                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                            <  1:740292/52‑1 (MQ=255)
                                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                            <  1:2658214/52‑1 (MQ=255)
                                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                            <  1:2083262/52‑1 (MQ=255)
                                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                            <  1:5078207/52‑1 (MQ=255)
                                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                            <  1:593299/52‑1 (MQ=255)
                                      aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                         >  1:3747987/1‑71 (MQ=255)
                                       aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                        <  1:2367583/71‑1 (MQ=255)
                                       aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                        <  1:3297121/71‑1 (MQ=255)
                                             ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCg                    >  1:3662170/1‑69 (MQ=255)
                                             ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                   >  1:4010258/1‑70 (MQ=255)
                                             ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                   >  1:3499390/1‑70 (MQ=255)
                                             ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                   >  1:2802122/1‑70 (MQ=255)
                                                 tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                         <  1:3641668/60‑1 (MQ=255)
                                                 tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                         <  1:3306006/60‑1 (MQ=255)
                                                 tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                         <  1:4690708/60‑1 (MQ=255)
                                                  aCGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                         <  1:1160143/59‑1 (MQ=255)
                                                   cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                       >  1:2760828/1‑60 (MQ=255)
                                                   cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                       >  1:294702/1‑60 (MQ=255)
                                                   cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                       >  1:4047195/1‑60 (MQ=255)
                                                   cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                       >  1:1581500/1‑60 (MQ=255)
                                                   cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                       >  1:1737997/1‑60 (MQ=255)
                                                   agTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                       >  1:461485/2‑60 (MQ=255)
                                                     ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAcc          <  1:1159017/71‑1 (MQ=255)
                                                     ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc           >  1:1283729/1‑70 (MQ=255)
                                                     ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc           >  1:4213329/1‑70 (MQ=255)
                                                       cAGAGCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc        >  1:1040516/1‑71 (MQ=255)
                                                       cAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc        >  1:3658936/1‑71 (MQ=255)
                                                        aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc        <  1:4016088/70‑1 (MQ=255)
                                                        aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc        <  1:3303827/70‑1 (MQ=255)
                                                        aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACg       <  1:3205649/71‑1 (MQ=255)
                                                        aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATACGGCTTCGCCATCAGACCAc        <  1:3090693/70‑1 (MQ=255)
                                                            cTTCACGGAAGTTCGCGGTGATTGGCGTTTCGAAGATGGAGCCATCCgg                         <  1:2967508/49‑1 (MQ=255)
                                                             ttCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                         <  1:4973791/48‑1 (MQ=255)
                                                              tCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa     <  1:886499/67‑1 (MQ=255)
                                                               cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa     <  1:4818177/66‑1 (MQ=255)
                                                               cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc  >  1:3823471/1‑69 (MQ=255)
                                                               cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc  >  1:1726566/1‑69 (MQ=255)
                                                                |                                                                   
TATCCGCCAGACCTTTACGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAC  >  Exported/2254682‑2254813

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: