Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 70,740 T→G 32.9% G118G (GGT→GGG araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904870,7400TG32.9% 20.6 / 65.2 61G118G (GGT→GGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (19/22);  new base G (11/9);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

TTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGG  >  AP009048/70674‑70809
                                                                  |                                                                     
ttACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                                   >  1:4245972/1‑71 (MQ=255)
 tACTTTCGTCCGCGCGCCTACTGGCATGATTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                   <  1:3516725/70‑1 (MQ=255)
 tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctt                                                                  >  1:3569318/1‑71 (MQ=255)
 tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctt                                                                  >  1:2963030/1‑71 (MQ=255)
 tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgtt                                                                     >  1:3964631/1‑68 (MQ=255)
 tACTTTCGTCCGCGCGCCTAATGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttctt                                                                  >  1:608339/1‑71 (MQ=255)
   cTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctttc                                                                <  1:1281467/71‑1 (MQ=255)
      tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGcc                                                             >  1:5367706/1‑71 (MQ=255)
      tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGcc                                                             >  1:617527/1‑71 (MQ=255)
        gTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCAGTCAATATTTGCCAATACGGGt                                                                       <  1:774923/59‑1 (MQ=255)
         tCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCgg                                                          <  1:2520259/71‑1 (MQ=255)
           cgcgcgCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGa                                                         >  1:2209278/1‑70 (MQ=255)
            gcgcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGAt                                                        >  1:5234993/1‑70 (MQ=255)
            gcgcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGa                                                         >  1:4547956/1‑69 (MQ=255)
            gcgcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGAt                                                        >  1:2362488/1‑70 (MQ=255)
             cgcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGa                                                      <  1:1549176/71‑1 (MQ=255)
              gcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATg                                                       >  1:3172315/1‑69 (MQ=255)
                 ccTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAgcg                                                  <  1:4246650/71‑1 (MQ=255)
                 ccTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAgcg                                                  <  1:809638/71‑1 (MQ=255)
                  cTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGAt                                                        >  1:1380674/1‑64 (MQ=255)
                    aCTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttctttc                                                                >  1:4580242/1‑54 (MQ=255)
                     cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctttc                                                                >  1:5219854/1‑53 (MQ=255)
                     cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCAc                                               <  1:3375814/70‑1 (MQ=255)
                     cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCAc                                               >  1:3601465/1‑70 (MQ=255)
                      tGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgtt                                                                     <  1:1534842/47‑1 (MQ=255)
                      tGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAg                                                    >  1:439725/1‑64 (MQ=255)
                         cATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGcc                                          >  1:802254/1‑71 (MQ=255)
                            gAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgtt                                                                     >  1:3117057/1‑41 (MQ=255)
                            gAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgtt                                                                     >  1:440355/1‑41 (MQ=255)
                             aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGaa                                                     >  1:4128756/1‑56 (MQ=255)
                              aTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAg                                            <  1:3398838/64‑1 (MQ=255)
                               tGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                    <  1:2919303/71‑1 (MQ=255)
                               tGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                    <  1:3769807/71‑1 (MQ=255)
                                ggCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGa                                                      <  1:5344294/52‑1 (MQ=255)
                                 gCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                    >  1:847941/1‑69 (MQ=255)
                                 gCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                    >  1:1488631/1‑69 (MQ=255)
                                 gCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCa                                  <  1:5254066/71‑1 (MQ=255)
                                  cTTAACTGGCCGTCAATATTTGCCAATACGGGgttctt                                                                  <  1:3956221/38‑1 (MQ=255)
                                   ttAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCg                                                           <  1:5025662/44‑1 (MQ=255)
                                   ttAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCg                                                           >  1:1111312/1‑44 (MQ=255)
                                   ttAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGc                                >  1:2771243/1‑71 (MQ=255)
                                    tAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGc                                >  1:3622466/1‑70 (MQ=255)
                                        ttgCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAtt                                     <  1:5290682/59‑1 (MQ=255)
                                        tGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAtt                                     <  1:4710078/61‑1 (MQ=255)
                                          gCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTg                          <  1:100443/70‑1 (MQ=255)
                                           ccGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                        <  1:4527320/71‑1 (MQ=255)
                                            cGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCa                                  >  1:4054469/1‑60 (MQ=255)
                                            cGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGc                                >  1:2068848/1‑62 (MQ=255)
                                            cGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                       <  1:887984/71‑1 (MQ=255)
                                            cGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                       <  1:2872370/71‑1 (MQ=255)
                                             gTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                       <  1:940554/70‑1 (MQ=255)
                                              tCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGa                                                      <  1:812007/38‑1 (MQ=255)
                                              tCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTgg                     >  1:4074349/1‑71 (MQ=255)
                                              tCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCAcc                                              <  1:1531359/46‑1 (MQ=255)
                                               cAATATTTGCCAATAAGGGTTTCTTTCGCCCGGATGAAGCGCACCAGc                                           >  1:979936/1‑48 (MQ=255)
                                                aaTATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTggg                    >  1:807900/1‑70 (MQ=255)
                                                 atatTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGc                   >  1:839084/1‑70 (MQ=255)
                                                  tatTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                         <  1:259531/47‑1 (MQ=255)
                                                    tttGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaaa                <  1:3653909/70‑1 (MQ=255)
                                                    tttGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAAt               <  1:1511980/71‑1 (MQ=255)
                                                      tGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc              <  1:5172645/70‑1 (MQ=255)
                                                      tGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCa             >  1:4167319/1‑71 (MQ=255)
                                                       gCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCa             <  1:3753525/70‑1 (MQ=255)
                                                        cctaTACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAt            <  1:290759/67‑1 (MQ=255)
                                                        ccAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAt            <  1:297053/70‑1 (MQ=255)
                                                        ccAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAtt           <  1:922093/71‑1 (MQ=255)
                                                         cAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAtt           >  1:1312382/1‑70 (MQ=255)
                                                           aTACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTaa         >  1:4581635/1‑70 (MQ=255)
                                                           aTACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAAc        >  1:791922/1‑71 (MQ=255)
                                                            tACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAt            <  1:793115/66‑1 (MQ=255)
                                                            tACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAt            <  1:2492207/66‑1 (MQ=255)
                                                            tACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTa          <  1:2290674/68‑1 (MQ=255)
                                                              cGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAAc        >  1:1474242/1‑68 (MQ=255)
                                                                 gTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                    >  1:3454044/1‑37 (MQ=255)
                                                                 gTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTgg                     >  1:5310094/1‑52 (MQ=255)
                                                                  tttctttcGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACTCCggg  >  1:173393/1‑70 (MQ=255)
                                                                  |                                                                     
TTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGG  >  AP009048/70674‑70809

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: