Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 71,175 A→G 24.6% R263R (CGA→CGG araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904871,1750AG24.6% 118.0 / 40.8 65R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base G (7/9);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATG  >  AP009048/71106‑71245
                                                                     |                                                                      
cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                       <  1:2449395/71‑1 (MQ=255)
  cTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                     >  1:3169570/1‑71 (MQ=255)
    aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                   >  1:3379054/1‑71 (MQ=255)
    aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                   >  1:2275612/1‑71 (MQ=255)
     cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                                  <  1:3341393/71‑1 (MQ=255)
     cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAttt                                                                  <  1:496871/71‑1 (MQ=255)
      ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                                 <  1:2487123/71‑1 (MQ=255)
       cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                                 >  1:1026763/1‑70 (MQ=255)
        ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                                >  1:2535314/1‑70 (MQ=255)
        ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                                >  1:4156967/1‑70 (MQ=255)
         gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:4783101/71‑1 (MQ=255)
          aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                             <  1:4630084/71‑1 (MQ=255)
          aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              >  1:859365/1‑70 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          <  1:4646135/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                           <  1:568495/70‑1 (MQ=255)
              cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                           <  1:2514942/69‑1 (MQ=255)
              cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                          >  1:1082046/1‑70 (MQ=255)
               tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:2494171/65‑1 (MQ=255)
               tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                                <  1:1382987/63‑1 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                            >  1:4287596/1‑66 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                       >  1:3811489/1‑71 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                       >  1:1942463/1‑71 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                                  >  1:2680847/1‑59 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                      >  1:1250698/1‑71 (MQ=255)
                 tCGACACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                                  >  1:4574537/1‑59 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                                     <  1:2604100/70‑1 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                                     <  1:1105142/70‑1 (MQ=255)
                    ccACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                               <  1:3777040/59‑1 (MQ=255)
                     cACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                          <  1:2817506/63‑1 (MQ=255)
                      aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                             >  1:4431477/1‑59 (MQ=255)
                        cGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                                 >  1:1534111/1‑69 (MQ=255)
                         gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                             <  1:4072564/56‑1 (MQ=255)
                          tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                                     <  1:63553/63‑1 (MQ=255)
                             gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAt                                                                    <  1:675045/45‑1 (MQ=255)
                             gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                                 <  1:107940/64‑1 (MQ=255)
                             gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                           >  1:1818695/1‑70 (MQ=255)
                              tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                             <  1:2813484/51‑1 (MQ=255)
                              tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:36291/50‑1 (MQ=255)
                              tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcga                                         <  1:5029984/71‑1 (MQ=255)
                                 cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                       >  1:978387/1‑54 (MQ=255)
                                  aaTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                        >  1:1474673/1‑52 (MQ=255)
                                  aaTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                       >  1:3884949/1‑69 (MQ=255)
                                   aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                    >  1:644836/1‑71 (MQ=255)
                                    tGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                                   >  1:4628287/1‑71 (MQ=255)
                                       tGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                                  <  1:3773550/69‑1 (MQ=255)
                                       aggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttcc                                >  1:190060/2‑69 (MQ=255)
                                         gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                               <  1:1642534/70‑1 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                            >  1:3332274/1‑70 (MQ=255)
                                                 gATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                          <  1:914709/35‑1 (MQ=255)
                                                  aTCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgt                     <  1:1110629/71‑1 (MQ=255)
                                                   tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                        >  1:4375288/1‑35 (MQ=255)
                                                   tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                                 >  1:2414356/1‑42 (MQ=255)
                                                   tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                       >  1:1668604/1‑52 (MQ=255)
                                                   tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                    <  1:4623216/71‑1 (MQ=255)
                                                         ctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa               >  1:1989279/1‑70 (MQ=255)
                                                          tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa               <  1:3652301/69‑1 (MQ=255)
                                                            tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa              <  1:471114/68‑1 (MQ=255)
                                                            tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa              <  1:4197566/68‑1 (MQ=255)
                                                                 cgcgAGTATTTAAAAAATGCACCGTGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat      >  1:2139773/1‑71 (MQ=255)
                                                                  gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa        <  1:1675092/68‑1 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGCGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat      >  1:5186752/1‑70 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGaa       <  1:4342664/69‑1 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatg     <  1:1128197/71‑1 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatg     <  1:2819502/71‑1 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat      <  1:3921137/70‑1 (MQ=255)
                                                                     aGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg  <  1:2243107/71‑1 (MQ=255)
                                                                     aGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg  <  1:4187956/71‑1 (MQ=255)
                                                                     |                                                                      
CACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATG  >  AP009048/71106‑71245

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: