Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,982,814 C→G 6.9% intergenic (‑765/+34) gadE ← / ← hdeD DNA‑binding transcriptional activator/acid‑resistance membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,982,8140CG 6.9% 177.4 / 2.7 58intergenic (‑765/+34)gadE/hdeDDNA‑binding transcriptional activator/acid‑resistance membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/31);  new base G (2/2);  total (25/33)
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

CAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGAA  >  AP009048/3982745‑3982883
                                                                     |                                                                     
cAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATc                                                                       <  1:2426533/70‑1 (MQ=255)
cAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATc                                                                       >  1:1310533/1‑70 (MQ=255)
   tGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGt                                                                   <  1:321387/71‑1 (MQ=255)
     aTAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGa                                                                 <  1:3775893/71‑1 (MQ=255)
      tAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGAc                                                                <  1:2635735/71‑1 (MQ=255)
        gTGCGTAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTa                                                              <  1:2342662/71‑1 (MQ=255)
        gTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACt                                                               <  1:4652625/70‑1 (MQ=255)
              aaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGt                                                                   >  1:2821074/1‑60 (MQ=255)
              aaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATcc                                                        <  1:900854/71‑1 (MQ=255)
              aaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATcc                                                        <  1:3765548/71‑1 (MQ=255)
              aaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATcc                                                        >  1:1291582/1‑71 (MQ=255)
                ttGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACTATg                                                            >  1:4168814/1‑65 (MQ=37)
                ttGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAg                                                                    <  1:2425425/57‑1 (MQ=255)
                 tGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCggg                                                     >  1:4418086/1‑71 (MQ=255)
                 tGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCgg                                                      <  1:4714389/70‑1 (MQ=255)
                  gATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGt                                                    >  1:2345667/1‑71 (MQ=255)
                   aTCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGt                                                    >  1:5272595/1‑70 (MQ=255)
                     cTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGa                                                  >  1:53373/1‑70 (MQ=255)
                      tACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCa                                                                      <  1:4056431/49‑1 (MQ=255)
                      tACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACt                                                               <  1:5130193/56‑1 (MQ=255)
                       acaacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACt                                                               <  1:3144672/55‑1 (MQ=25)
                       acaacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATc                                                                       <  1:1149525/47‑1 (MQ=255)
                        caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACt                                                               >  1:4096504/1‑54 (MQ=25)
                        caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTg                                                   <  1:1622108/66‑1 (MQ=255)
                        caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGAc                                                 <  1:295345/68‑1 (MQ=255)
                        caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAAc                                              >  1:5314032/1‑71 (MQ=255)
                        caacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAAc                                              >  1:1200943/1‑71 (MQ=255)
                         aacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAAc                                              <  1:3472646/70‑1 (MQ=255)
                          acacTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTa                                                              >  1:520138/1‑53 (MQ=255)
                          acacTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCg                                            <  1:3477898/71‑1 (MQ=255)
                           cacTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCgg                                           >  1:1261545/1‑71 (MQ=255)
                           cacTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCgg                                           >  1:2452314/1‑71 (MQ=255)
                                 gTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACt                                                               <  1:2953008/45‑1 (MQ=18)
                                  tAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAatt                                    <  1:1033170/71‑1 (MQ=255)
                                    gCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAatta                                   <  1:362066/70‑1 (MQ=255)
                                    gCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAatta                                   <  1:29412/70‑1 (MQ=255)
                                    gCGGAGAGAGTATTAATAGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTaa                                      >  1:1187932/1‑67 (MQ=255)
                                     cGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATg                                                            <  1:4764817/44‑1 (MQ=255)
                                     cGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAattatt                                 <  1:2907352/71‑1 (MQ=255)
                                         gagagTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCggg                                          >  1:3676437/1‑58 (MQ=255)
                                         gagagTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATtgctg                             >  1:1363028/1‑71 (MQ=255)
                                             gTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTa                         <  1:1035095/71‑1 (MQ=255)
                                              tATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACaa                                               >  1:3606690/1‑48 (MQ=255)
                                                    tCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAAc                   >  1:2303733/1‑70 (MQ=255)
                                                     cGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGaa                    <  1:807912/68‑1 (MQ=255)
                                                       gATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAAc               <  1:2917406/71‑1 (MQ=255)
                                                          cATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTgg            >  1:2062430/1‑71 (MQ=255)
                                                           aTAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGc           >  1:3895263/1‑71 (MQ=255)
                                                             aGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACTAACAAACTGGCg          >  1:3317987/1‑70 (MQ=255)
                                                             aGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTgg            >  1:2711813/1‑68 (MQ=255)
                                                              gTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGa         <  1:5096878/70‑1 (MQ=255)
                                                               tCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAattat                                  <  1:3913836/44‑1 (MQ=255)
                                                               tCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATtgc                               <  1:4675342/47‑1 (MQ=255)
                                                               tCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACt              <  1:277497/64‑1 (MQ=255)
                                                                 acaTCAAGTGACTATGATCCGGGTGACAACCGGGGTAatta                                   >  1:4628273/1‑41 (MQ=255)
                                                                 acaTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGc      <  1:4272730/70‑1 (MQ=255)
                                                                 acaTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCt     <  1:4002751/71‑1 (MQ=255)
                                                                  caTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTg    <  1:1269719/71‑1 (MQ=255)
                                                                   aTCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTg    >  1:1824853/1‑70 (MQ=255)
                                                                    tCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACa                  <  1:4070923/55‑1 (MQ=255)
                                                                    tCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGaa  >  1:2371478/1‑71 (MQ=255)
                                                                     cAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACaaa                >  1:1352515/1‑56 (MQ=255)
                                                                     cAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGaa  >  1:4022068/1‑70 (MQ=255)
                                                                     |                                                                     
CAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGAA  >  AP009048/3982745‑3982883

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: