Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,254,623 C→T 12.0% Q805Q (CAG→CAA rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,254,6230CT12.0% 224.8 / 17.7 75Q805Q (CAG→CAArpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (27/39);  new base T (3/6);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CACCCTCGATAACCGGAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCA  >  Exported/2254553‑2254690
                                                                      |                                                                   
caccCTCGATAACCGGAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTcc                                                                     <  1:499356/71‑1 (MQ=255)
  cccTCGATAACCGGAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTg                                                                   >  1:867049/1‑71 (MQ=255)
    cTCGATAACCGGAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTgc                                                                  >  1:2919490/1‑70 (MQ=255)
              ggAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAc                                                       <  1:95799/71‑1 (MQ=255)
              ggAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAc                                                       <  1:3140067/71‑1 (MQ=255)
              ggAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAc                                                       <  1:3059152/71‑1 (MQ=255)
               gAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAcc                                                      <  1:2241685/71‑1 (MQ=255)
                     tcatGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacga                                                >  1:933175/1‑71 (MQ=255)
                       atGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                               >  1:2681904/1‑70 (MQ=255)
                       atGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGacgacg                                              <  1:617799/71‑1 (MQ=255)
                       atGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGacgacg                                              <  1:29454/71‑1 (MQ=255)
                       atGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGacgacg                                              <  1:2885991/71‑1 (MQ=255)
                        tGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgac                                                  <  1:2115601/66‑1 (MQ=255)
                        tGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgac                                                  <  1:3021307/66‑1 (MQ=255)
                        tGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgac                                                  <  1:3220793/66‑1 (MQ=255)
                         gATACCTTCATGGGTACCACACTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTc                                                          >  1:2166741/1‑57 (MQ=255)
                         gATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGcc                                                               <  1:2224104/52‑1 (MQ=255)
                         gATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGcc                                                               <  1:1084944/52‑1 (MQ=255)
                         gATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTc                                                          >  1:3005236/1‑57 (MQ=255)
                          aTACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAg                                            >  1:1895339/1‑70 (MQ=255)
                          aTACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGt                                           >  1:1057772/1‑71 (MQ=255)
                          aTACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGt                                           >  1:368078/1‑71 (MQ=255)
                           tACCTACATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGt                                           >  1:1361052/1‑70 (MQ=255)
                                tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                             >  1:3084721/1‑63 (MQ=255)
                                tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                             >  1:3122633/1‑63 (MQ=255)
                                tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                             >  1:1817848/1‑63 (MQ=255)
                                tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGt                                      >  1:2967108/1‑70 (MQ=255)
                                tCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGt                                      >  1:1600744/1‑70 (MQ=255)
                                 cATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGt                                      <  1:3270563/69‑1 (MQ=255)
                                 cATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGAAGAGTCAGGTa                                     <  1:714221/70‑1 (MQ=255)
                                  aTGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAAc                                   <  1:1723953/71‑1 (MQ=255)
                                        aCCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGt                             >  1:1814300/1‑71 (MQ=255)
                                         ccACAATCGTCTTCGGTAACCACCAGGTTCTGCGCCACGGCAACCAGACGACGAGTCAGGTAACCGGAGt                             <  1:478107/70‑1 (MQ=255)
                                          cacaATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAAc                                                       <  1:3060457/43‑1 (MQ=255)
                                          cacaATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTc                           >  1:2811659/1‑71 (MQ=255)
                                            caATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCg                          >  1:2517912/1‑70 (MQ=255)
                                             aaTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTc                                          <  1:2761326/53‑1 (MQ=255)
                                              aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTa                                     <  1:1391619/57‑1 (MQ=255)
                                              aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTa                                     <  1:411812/57‑1 (MQ=255)
                                              aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAg                       <  1:746673/71‑1 (MQ=255)
                                              aTCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAg                       <  1:3231330/71‑1 (MQ=255)
                                               tCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGt                      >  1:1741901/1‑71 (MQ=255)
                                               tCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGt                      >  1:2539298/1‑71 (MQ=255)
                                                  tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTa                                     <  1:145764/53‑1 (MQ=255)
                                                  tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                   <  1:1208562/71‑1 (MQ=255)
                                                  tCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                   <  1:3479414/71‑1 (MQ=255)
                                                    ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                               >  1:2283414/1‑41 (MQ=255)
                                                    ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                               >  1:454034/1‑41 (MQ=255)
                                                    ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                               >  1:2151196/1‑41 (MQ=255)
                                                    ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAgacgac                                               >  1:427744/1‑41 (MQ=255)
                                                    ttCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAAc                                   >  1:281974/1‑53 (MQ=255)
                                                      cGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAg                                        >  1:1319297/1‑46 (MQ=255)
                                                        gTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTg              <  1:846624/70‑1 (MQ=255)
                                                        gTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTg              <  1:2128214/70‑1 (MQ=255)
                                                         tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                          <  1:3329271/39‑1 (MQ=16)
                                                         tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                          <  1:3434053/39‑1 (MQ=16)
                                                         tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                          <  1:801822/39‑1 (MQ=16)
                                                         tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                          <  1:691904/39‑1 (MQ=16)
                                                         tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                          <  1:196598/39‑1 (MQ=16)
                                                         tcaccaccAGATCTTGCGTCACGTCCACCAGACGACGGGTc                                          <  1:606162/39‑1 (MQ=16)
                                                         tAACCACCAGGTCCTGCGCCACGTCAACCAGAcgacga                                             <  1:428827/38‑1 (MQ=255)
                                                         tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                <  1:1200365/51‑1 (MQ=255)
                                                         tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                <  1:2749328/51‑1 (MQ=255)
                                                         tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                <  1:2208723/51‑1 (MQ=255)
                                                         tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                <  1:2330950/51‑1 (MQ=255)
                                                         tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                <  1:2366973/51‑1 (MQ=255)
                                                         tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                <  1:1947422/51‑1 (MQ=255)
                                                         tAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCg            >  1:762294/1‑71 (MQ=255)
                                                           accaccAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCgg                                <  1:660541/49‑1 (MQ=255)
                                                           accaccAGATCTTGCGTCACGTCCACCAGACGACGGGTCAg                                        >  1:1191944/1‑41 (MQ=19)
                                                           accaccAGATCTTGCGTCACGTCCACCAGACGACGGGTCAg                                        >  1:2577199/1‑41 (MQ=19)
                                                           accaccAGATCTTGCGTCACGTCCACCAGACGACGGGTCAg                                        >  1:2836454/1‑41 (MQ=19)
                                                                   gTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGtttt                   <  1:3193477/54‑1 (MQ=255)
                                                                   gTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCa  <  1:1529348/71‑1 (MQ=255)
                                                                   gTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCa  <  1:1632490/71‑1 (MQ=255)
                                                                   gTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCa  <  1:2112265/71‑1 (MQ=255)
                                                                      |                                                                   
CACCCTCGATAACCGGAGTCATCATGATACCTTCATGGGTACCACAATCGTCTTCGGTAACCACCAGGTCCTGCGCCACGTCAACCAGACGACGAGTCAGGTAACCGGAGTTCGCAGTTTTCAGTGCGGTATCCGCCA  >  Exported/2254553‑2254690

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: