Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,982,812 A→T 7.2% intergenic (‑763/+36) 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,8120AT 7.2% 153.1 / 2.8 55intergenic (‑763/+36)gadE/hdeDDNA‑binding transcriptional activator/acid‑resistance membrane protein
Reads supporting (aligned to +/- strand):  ref base A (22/29);  new base T (2/2);  total (24/31)
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

CAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTG  >  AP009048/3982745‑3982881
                                                                   |                                                                     
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)
           cgcAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACAt                                                                      >  1:1269010/1‑58 (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)
                                                                   |                                                                     
CAATGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTG  >  AP009048/3982745‑3982881

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: