Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,982,816 A→T 6.8% intergenic (‑767/+32) 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,8160AT 6.8% 163.3 / 4.1 59intergenic (‑767/+32)gadE/hdeDDNA‑binding transcriptional activator/acid‑resistance membrane protein
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base T (2/2);  total (27/32)
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.61e-01

TGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGAACA  >  AP009048/3982748‑3982885
                                                                    |                                                                     
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:3765548/71‑1 (MQ=255)
           aaTTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATcc                                                          <  1:900854/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)
                   tACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACt                                                                 <  1:5130193/56‑1 (MQ=255)
                    acaacaCTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACt                                                                 <  1:3144672/55‑1 (MQ=25)
                     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:2452314/1‑71 (MQ=255)
                        cacTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCgg                                             >  1:1261545/1‑71 (MQ=255)
                              gTAGCGGAGAGAGTATTAATCGGATCATAGTCACTTGATGTGACt                                                                 <  1:2953008/45‑1 (MQ=18)
                               tAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAatt                                      <  1:1033170/71‑1 (MQ=255)
                                 gCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAatta                                     <  1:29412/70‑1 (MQ=255)
                                 gCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAatta                                     <  1:362066/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)
                                                                   aaGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGaa    >  1:2331372/1‑69 (MQ=255)
                                                                    aGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGAACa  <  1:4051242/70‑1 (MQ=255)
                                                                    aGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGAACa  <  1:1936196/70‑1 (MQ=255)
                                                                    |                                                                     
TGATAGTGCGCAATTGATCTACAACACTGCGTAGCGGAGAGAGTATTAATCGGATCATAGTCACATCAAGTGACTATGATCCGGGTGACAACCGGGGTAATTATTGCTGCTTAACGAACAAACTGGCGAAGCTGAACA  >  AP009048/3982748‑3982885

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: