Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,256,639 G→C 5.5% R133R (CGC→CGG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,256,6390GC 5.5% 234.9 / 3.0 72R133R (CGC→CGGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (26/42);  new base C (2/2);  total (28/44)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

ACGTTCCAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTC  >  Exported/2256570‑2256708
                                                                     |                                                                     
aCGTTCCAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATcg                                                                       <  1:2621594/70‑1 (MQ=255)
     ccAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAg                                                                  >  1:636292/1‑70 (MQ=255)
     ccAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGc                                                                 >  1:165971/1‑71 (MQ=255)
     ccAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGc                                                                 >  1:2018698/1‑71 (MQ=255)
     ccAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGc                                                                 >  1:750534/1‑71 (MQ=255)
             gTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCatat                                                          <  1:23661/70‑1 (MQ=255)
             gTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATc                                                         <  1:1309723/71‑1 (MQ=255)
              tCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCg                                                        <  1:1901823/71‑1 (MQ=255)
               cATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCg                                                        >  1:3527816/1‑70 (MQ=255)
                aTACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCa                                                                   >  1:1123132/1‑58 (MQ=255)
                aTACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGag                                                      >  1:3196960/1‑71 (MQ=255)
                aTACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGag                                                      >  1:2528396/1‑71 (MQ=255)
                aTACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGag                                                      >  1:889063/1‑71 (MQ=255)
                 tACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagc                                                     <  1:3640092/71‑1 (MQ=255)
                 tACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagc                                                     <  1:1273968/71‑1 (MQ=255)
                 tACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagc                                                     <  1:2664361/71‑1 (MQ=255)
                 tACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagc                                                     <  1:2346255/71‑1 (MQ=255)
                   ccGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagca                                                    <  1:837355/70‑1 (MQ=255)
                   ccGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagca                                                    <  1:2416458/70‑1 (MQ=255)
                    cGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGagcagc                                                  <  1:2322040/71‑1 (MQ=255)
                     gCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAgcagca                                                 <  1:3274822/71‑1 (MQ=255)
                            aTAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATcgcg                                                                     <  1:853271/44‑1 (MQ=255)
                            aTAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATcgcg                                                                     <  1:97646/44‑1 (MQ=255)
                            aTAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATcgcg                                                                     <  1:1221695/44‑1 (MQ=255)
                            aTAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATcgcg                                                                     <  1:2076947/44‑1 (MQ=255)
                            aTAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGAt                                          <  1:2156873/71‑1 (MQ=255)
                             tAACCACATAGGATTCGAAGTACAGTACGCGTTCGATATcgcg                                                                     <  1:2342604/43‑1 (MQ=255)
                             tAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATa                                         <  1:687848/71‑1 (MQ=255)
                              aaCCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATAc                                        >  1:159683/1‑71 (MQ=255)
                              aaCCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATAc                                        >  1:3023021/1‑71 (MQ=255)
                               aCCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACg                                       >  1:2767317/1‑71 (MQ=255)
                               aCCACATAGGATTCAAAGTACAGTACGCGCTCGATATGCCGCTGCgc                                                               <  1:1204921/47‑2 (MQ=21)
                               aCCACATAGGATTCAAAGTACAGTACGCGCTCGATATGCCGCTGCgc                                                               <  1:1617076/47‑2 (MQ=21)
                                ccACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACgg                                      <  1:2382304/71‑1 (MQ=255)
                                ccACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACgg                                      <  1:1426142/71‑1 (MQ=255)
                                  acaTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGa                                    <  1:2567060/71‑1 (MQ=255)
                                   caTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGc                                                              <  1:607181/44‑1 (MQ=255)
                                   caTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGa                                    >  1:1939258/1‑70 (MQ=255)
                                    aTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATc                                                         >  1:3602075/1‑48 (MQ=255)
                                    aTAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATa                                         <  1:464415/64‑1 (MQ=255)
                                         aTTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGa                                    >  1:742142/1‑64 (MQ=255)
                                         aTTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGa                                    >  1:496180/1‑64 (MQ=255)
                                          ttCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACgg                                      <  1:1502088/61‑1 (MQ=255)
                                            cAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGa                           <  1:3724320/70‑1 (MQ=255)
                                            cAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAt                          >  1:735935/1‑71 (MQ=255)
                                             aaaGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAt                          <  1:657277/70‑1 (MQ=255)
                                             aaaGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAt                          <  1:1880169/70‑1 (MQ=255)
                                              aaGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGAc                                              <  1:741210/49‑1 (MQ=255)
                                              aaGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAttt                        <  1:2599756/71‑1 (MQ=255)
                                                    agtacGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATa                                         <  1:2295406/48‑1 (MQ=255)
                                                      tacGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGAc                                              <  1:3396431/41‑1 (MQ=255)
                                                      tacGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAAcc                <  1:3590960/71‑1 (MQ=255)
                                                      tacGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAAc                 >  1:83095/1‑70 (MQ=255)
                                                        cgcgTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGTAAc                 >  1:1220244/1‑68 (MQ=255)
                                                         gcgTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGa                   <  1:1139027/65‑1 (MQ=255)
                                                           gTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCg                                            <  1:1513901/38‑1 (MQ=255)
                                                           gTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCg                                            <  1:2409202/38‑1 (MQ=255)
                                                           gTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCg                                            <  1:543964/38‑1 (MQ=255)
                                                           gTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAttt                        <  1:138096/58‑1 (MQ=255)
                                                            ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGAt                          <  1:2313438/55‑1 (MQ=255)
                                                            ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCATGAACCAGAtgt          <  1:959721/71‑1 (MQ=255)
                                                            ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAg              <  1:1393936/67‑1 (MQ=255)
                                                            ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAg              <  1:3736660/67‑1 (MQ=255)
                                                            ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAg              <  1:304898/67‑1 (MQ=255)
                                                            ttCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAtgt          <  1:2965208/71‑1 (MQ=255)
                                                            ctCGATATGCCGCTGCGCGATATCGAACAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAtgt          >  1:3287004/2‑71 (MQ=21)
                                                            ctCGATATGCCGCTGCGCGATATCGAACAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAtgt          >  1:944790/2‑71 (MQ=21)
                                                              cGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAtgtg         >  1:390952/1‑70 (MQ=255)
                                                              cGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgc        >  1:896976/1‑71 (MQ=255)
                                                              cGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgc        >  1:1595031/1‑71 (MQ=255)
                                                              cGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgc        >  1:2152800/1‑71 (MQ=255)
                                                               gATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACgg                                      <  1:2996281/40‑1 (MQ=255)
                                                               gATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgcg       >  1:1532327/1‑71 (MQ=255)
                                                               gATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTgcg       >  1:2285704/1‑71 (MQ=255)
                                                                 tatCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGAt            >  1:2201370/1‑64 (MQ=255)
                                                                  atCGCGCAGCGGCATATCGAGCAGCAGAGCGATACGGGACGGCAGCGATTTCAGGAAcc                <  1:556636/59‑1 (MQ=255)
                                                                  atCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAAcc                <  1:337916/59‑1 (MQ=255)
                                                                  atCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAAcc                <  1:2727854/59‑1 (MQ=255)
                                                                    cgcgCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTc  >  1:2990736/1‑71 (MQ=255)
                                                                    cgcgCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTc  >  1:2186607/1‑71 (MQ=255)
                                                                     gcgcAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTc  <  1:2587135/70‑1 (MQ=255)
                                                                     |                                                                     
ACGTTCCAGGTTGGTCATACCGCCTTCGATAACCACATAGGATTCAAAGTACAGTACGCGTTCGATATCGCGCAGCGGCATATCGAGCAGCAGACCGATACGGGACGGCAGCGATTTCAGGAACCAGATGTGCGCAGTC  >  Exported/2256570‑2256708

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: