Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 71,214 T→C 29.9% F276F (TTT→TTC)  araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904871,2140TC29.9% 107.4 / 63.8 77F276F (TTT→TTC) araCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base C (12/11);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGA  >  AP009048/71147‑71281
                                                                   |                                                                   
gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                                                                   <  1:4748361/70‑1 (MQ=255)
   ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                                                                <  1:5128710/70‑1 (MQ=255)
   ttGAAGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                                                                <  1:2964427/70‑1 (MQ=255)
          tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt                                                         >  1:317514/1‑70 (MQ=255)
                   tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa                                                  <  1:4615251/68‑1 (MQ=255)
                   tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa                                                  <  1:1400938/68‑1 (MQ=255)
                     tctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt                                              >  1:549115/1‑70 (MQ=255)
                     tctcGCGAGTATTTAAAAAATGCACAGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt                                              <  1:2729000/70‑1 (MQ=255)
                      ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg                                             <  1:583383/70‑1 (MQ=255)
                      ctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTg                                             <  1:2639173/70‑1 (MQ=255)
                          cgAGAATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                                                        >  1:5367163/1‑55 (MQ=255)
                           gAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgat                                       >  1:1397413/1‑71 (MQ=255)
                            aGTATTTAAAAAATGCACCGGGGCCAGCACGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg                                      <  1:1097584/71‑1 (MQ=255)
                              tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaag                                                     <  1:591246/54‑1 (MQ=255)
                               aTTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa                                                  <  1:5370207/56‑1 (MQ=255)
                                tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGc                                  >  1:2906428/1‑71 (MQ=255)
                                  tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgt                                                         >  1:2328841/1‑46 (MQ=255)
                                  tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa                                                  <  1:4017268/53‑1 (MQ=255)
                                  tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGcc                                 >  1:3477600/1‑70 (MQ=255)
                                   aaaaaaTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                                                                >  1:1987420/1‑38 (MQ=255)
                                   aaaaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGt                               >  1:656817/1‑71 (MQ=255)
                                    aaaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCg                                <  1:1818348/69‑1 (MQ=255)
                                     aaaaTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                              <  1:1471173/70‑1 (MQ=255)
                                     aaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                              >  1:5274842/1‑70 (MQ=255)
                                     aaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCa                             <  1:344445/71‑1 (MQ=255)
                                      aaaTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                              >  1:1431764/1‑69 (MQ=255)
                                         tGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGtt                         >  1:3787250/1‑71 (MQ=255)
                                           cACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTg                        <  1:1362352/70‑1 (MQ=255)
                                           cACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGt                       >  1:3783011/1‑71 (MQ=255)
                                            aCCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa                                                  <  1:2270166/43‑1 (MQ=255)
                                            aCCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGt                       >  1:170011/1‑70 (MQ=255)
                                            aCCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGt                       >  1:4322474/1‑70 (MQ=255)
                                              cGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCa                     >  1:2322955/1‑70 (MQ=255)
                                              cGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCa                     >  1:5293410/1‑70 (MQ=255)
                                               ggggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaaga                                                    <  1:2728928/38‑1 (MQ=255)
                                               ggggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                    <  1:4154375/70‑1 (MQ=255)
                                               ggggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATa                   <  1:2749120/71‑1 (MQ=255)
                                               ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATa                   <  1:2727368/71‑1 (MQ=255)
                                                 ggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAt                 >  1:4214309/1‑71 (MQ=255)
                                                  gCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAag                >  1:2480447‑M1/1‑69 (MQ=255)
                                                   ccAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa                                                 <  1:2634455/37‑1 (MQ=255)
                                                   ccAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGt                       <  1:1004798/63‑1 (MQ=255)
                                                   ccAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATa                   <  1:362209/67‑1 (MQ=255)
                                                   ccAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAt                 <  1:1222584/69‑1 (MQ=255)
                                                   ccAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAt                 <  1:2448243/69‑1 (MQ=255)
                                                    cAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                              <  1:513820/55‑1 (MQ=255)
                                                    cAGCCCGAGCGAGTTCCGTGCCGGTTGTTAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAag                >  1:4868877‑M1/1‑67 (MQ=255)
                                                    cAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGt                          <  1:909659/59‑1 (MQ=255)
                                                     aGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGt             <  1:802790/71‑1 (MQ=255)
                                                      gCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                              <  1:1383579/53‑1 (MQ=255)
                                                      gCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAg                           >  1:960183/1‑56 (MQ=255)
                                                       cccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGATTGATGTAGCCGTCAAGt                          >  1:2726838/1‑56 (MQ=255)
                                                       cccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAt                 >  1:4097156/1‑65 (MQ=255)
                                                       cccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtcc           >  1:4928918‑M1/1‑64 (MQ=255)
                                                        ccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAAc          >  1:4034351/1‑71 (MQ=255)
                                                        ccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGCCATAAagggtccc          >  1:1384858‑M1/1‑63 (MQ=255)
                                                         cgagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTgg              >  1:4169202/1‑66 (MQ=255)
                                                         cgagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAAc          >  1:4202623/1‑70 (MQ=255)
                                                         cgagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATaa                  >  1:839082/1‑62 (MQ=255)
                                                          gagcgagTTTCGTGCCGGTTGTGAAGAATAAGTGAATGATGTAGCCGGCAAGtt                         <  1:2198812/54‑1 (MQ=39)
                                                          gagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGaa                                           <  1:3357684/36‑1 (MQ=255)
                                                          gagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCaa                            <  1:4769908/51‑1 (MQ=255)
                                                          gagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACg         >  1:1626679/1‑70 (MQ=255)
                                                          gagcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACg         >  1:1632853/1‑70 (MQ=255)
                                                          gagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtcccta        <  1:2290621‑M1/71‑11 (MQ=255)
                                                           agcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGaa       >  1:189769/1‑71 (MQ=255)
                                                            gcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg                                      >  1:4979617/1‑39 (MQ=255)
                                                            gcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg                                      >  1:4397880/1‑39 (MQ=255)
                                                            gcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg                                      >  1:4382511/1‑39 (MQ=255)
                                                            gcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa            <  1:3403308/65‑1 (MQ=255)
                                                            gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtccctaag      >  1:143539‑M1/1‑59 (MQ=255)
                                                             cgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTc                              >  1:1389662/1‑46 (MQ=255)
                                                             cgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATaa                  <  1:4106545/58‑1 (MQ=255)
                                                               agTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtccctaagtct   <  1:3717012‑M1/71‑16 (MQ=255)
                                                                gTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGa  >  1:2596396/1‑71 (MQ=255)
                                                                gTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGa  >  1:5330138/1‑71 (MQ=255)
                                                                gTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGtt                         <  1:5084649/48‑1 (MQ=255)
                                                                 tttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAAc          <  1:3625724/62‑1 (MQ=255)
                                                                 tttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGa        <  1:2997585/64‑1 (MQ=255)
                                                                 tttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAg   <  1:3517817/69‑1 (MQ=255)
                                                                 tttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGa  >  1:3330476/1‑70 (MQ=255)
                                                                   tCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCa                     >  1:2520889/1‑49 (MQ=255)
                                                                   |                                                                   
GTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGA  >  AP009048/71147‑71281

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: ‑