Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,705,425 G→A 12.5% intergenic (‑163/+14) kdpF ← / ← PP_4162 potassium‑transporting ATPase F/conserved membrane protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,705,4250GA12.5% 84.6 / 5.7 32intergenic (‑163/+14)kdpF/PP_4162potassium‑transporting ATPase F/conserved membrane protein of unknown function
Reads supporting (aligned to +/- strand):  ref base G (20/8);  new base A (2/2);  total (22/10)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

GCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACG  >  AE015451/4705321‑4705649
                                                                                                        |                                                                                                                                                                                                                                
gCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGCGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATtgttgt                                                                                                                                                                                                           <  1:926131/128‑1 (MQ=255)
                            ggtgggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTc                                                                                                                                         >  1:832717/1‑166 (MQ=255)
                                             cGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGcc                                                                                                                                                                                                                              >  2:945948/1‑64 (MQ=255)
                                             cGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGcc                                                                                                                                                                                                                              <  1:945953/64‑1 (MQ=255)
                                             cGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCact                                                                                                                     >  1:26934/1‑167 (MQ=255)
                                              gTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGcc                                                                                                                                     >  1:978720/1‑152 (MQ=255)
                                              gTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTAACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAc                                                                                                                                                      >  2:663492/1‑135 (MQ=255)
                                               tAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCa                                                                                                              >  2:1108515/1‑174 (MQ=255)
                                                          gTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGa                                                                                           >  2:599443/1‑182 (MQ=255)
                                                           ttGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAg                                       <  1:1191845/233‑1 (MQ=255)
                                                                 tCTGCCTGTAACGGCCTCTTCGCGGACTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAg                                                                                                                                            >  1:724893/1‑126 (MQ=255)
                                                                    gCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAac                                                                                         >  2:1191840/1‑174 (MQ=255)
                                                                      cTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTGCAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGa                                                                                                                                  >  2:815745/1‑131 (MQ=255)
                                                                         tAACGGCCCTTTCGCGGGCTTGCCCGCGAAAGGGCCGCTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGGCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAAGCCAATGAAg                                                                                                                 >  1:1230783/1‑145 (MQ=39)
                                                                           aCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTa                                                                                                                                                                                                                          >  1:649410/1‑38 (MQ=255)
                                                                           aCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCcag                                                                                 >  2:792416/1‑175 (MQ=255)
                                                                               ccTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATc                                         <  1:1108520/211‑1 (MQ=255)
                                                                               ccTCTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAg                                       <  1:488075/213‑1 (MQ=255)
                                                                                 tctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCcaca              <  1:664632/236‑1 (MQ=255)
                                                                                  ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGcc                                                                                                                        <  2:1258445/129‑1 (MQ=255)
                                                                                  ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGcc                                                                                                                        >  1:1258450/1‑129 (MQ=255)
                                                                                  ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGaa                                                                                          >  2:524132/1‑159 (MQ=255)
                                                                                  ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGa                                                                                           >  2:1171120/1‑158 (MQ=255)
                                                                                  ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAaca                                                                                        >  2:1106620/1‑161 (MQ=255)
                                                                                  ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGcc                                                                                    >  2:663084/1‑165 (MQ=255)
                                                                                        gggCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACAGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTgccgc                                                    >  1:693206/1‑191 (MQ=255)
                                                                                           caagCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGa                                                                                           >  2:488070/4‑149 (MQ=255)
                                                                                             tGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACAcc            <  1:663089/226‑1 (MQ=255)
                                                                                                ccGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGctgctg                                                                                                                                                                                                >  2:462079/1‑43 (MQ=255)
                                                                                                ccGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACAGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACg  <  2:693201/233‑1 (MQ=255)
                                                                                                  gcgAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGt                   <  2:1093439/214‑1 (MQ=255)
                                                                                                  gcgAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGt                   >  1:1093444/1‑214 (MQ=255)
                                                                                                        |                                                                                                                                                                                                                                
GCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAGCACGATCACCAGCGGCAGGTTGCCGCCACGGCTGATCAGGTACGGTTTGAGCACGTTGTCCACACCGCTGATGACG  >  AE015451/4705321‑4705649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: