Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,705,460 T→A 15.0% intergenic (‑152/+25) kdpF ← / ← PP_4162 potassium‑transporting ATPase subunit F/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,705,4600TA15.0% 75.4 / 5.8 33intergenic (‑152/+25)kdpF/PP_4162potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (13/14);  new base A (2/3);  total (15/18)
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.25e-01

GTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACA  >  NC_002947/4705311‑4705609
                                                                                                                                                     |                                                                                                                                                     
gtgGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCtt                                                                                                                                                       >  2:638876/1‑150 (MQ=255)
 tgGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGGTGGGGTCTGCCTGGAACGTCCTCTTCGCGGGCTTg                                                                                                                                                      >  2:472607/1‑150 (MQ=255)
   gATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCGTCTTCGCGGGCTTGcc                                                                                                                                                    >  1:686097/1‑150 (MQ=255)
         ccggggtTAGGGCTGGGGCGGTTTTGGGGAAATTGTCGGGGGGAGGGGCGCAAAGGGCCGAAGGGGGGGGGGGGGTGTTGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                              <  2:627650/145‑1 (MQ=255)
         cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCCAGCCcgctca                                                                                                                                              >  2:204517/1‑147 (MQ=255)
         cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCCAGATCGGGGGATGGGGGTTGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              >  1:531212/1‑150 (MQ=255)
          ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                             >  1:587264/1‑150 (MQ=255)
            gggtTACGGCTTGGGCCGCTTTGGGGAAAGAGTCCGGGGGAAAGGAGTAAAGGGGGGAGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCAGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGaagag                                                                                                                                           <  1:561596/148‑4 (MQ=25)
             ggtTAGGTCTTGGGTGGTTTTGGTGAAATTGTCGGGGGGAGGGGGGTAAGGGGGCGAGGGGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGaagagg                                                                                                                                          <  1:204517/149‑5 (MQ=19)
               ttACGGTTGGGGCCGCTGGGGTGAAATTGCCCGAGGGGGTGGCGTAAGGGGGCGAGACCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                                                                                                        <  2:740674/150‑1 (MQ=255)
               ttACGGCTTGGGCCGCTTTGGGGAAATTGTCCGGGGGAGAGGCGTAAAGGGGCGGGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                                                                                                        <  2:587264/150‑1 (MQ=255)
                  cgccgtGGGCCGCTTTGGGGAAATTGTCCGAAGGAGAGGCGTAACGGGGCGAGGCGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGtt                                                                                                                                     <  1:29321/145‑1 (MQ=255)
                        gggCCGCTTTGGGGAAAGTGGCCGAGGGAGAGGCGTAAAGGGGCGAGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                                                                                                               <  2:382972/150‑1 (MQ=255)
                                     gAAATTGTCCGAGGGAGAGGCGTAAAGGGGCTAGATCGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgga                                                                                                                  >  1:19705/1‑149 (MQ=255)
                                            tCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgctgct                                                                                                           >  1:486221/1‑150 (MQ=255)
                                              ggggggAGGGGGGAAAAGGGGCGAGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCtgctgc                                                                                                         <  2:686097/147‑1 (MQ=255)
                                                 gggAGGGGGGTAAGGGGGGGAGGTGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGcgcc                                                                                                      <  2:827832/150‑1 (MQ=255)
                                                       ggCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTa                                                                                                <  2:531212/150‑1 (MQ=255)
                                                                 ggCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGTGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagt                                                                                      >  1:56716/1‑150 (MQ=255)
                                                                    gagaTGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGACGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                   <  1:236039/150‑1 (MQ=255)
                                                                     agaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtccag                                                                                  >  1:520951/1‑150 (MQ=255)
                                                                      gaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagc                                                                                 >  2:529336/1‑150 (MQ=255)
                                                                                 gggggtgggTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTa                                                                      <  2:56716/150‑1 (MQ=255)
                                                                                        tgtATAAAGAAAGCGTACACATCAGCGTTGTGGTCTGCCTGTAACGCCCTCTTCGCGGGCCGGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAcc                                                               <  1:162622/149‑1 (MQ=255)
                                                                                          tataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGAATGGAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGc                                                             <  1:472607/150‑1 (MQ=255)
                                                                                            taAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCa                                                           <  1:638876/150‑1 (MQ=255)
                                                                                                                    tgtgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGc                                   <  2:486221/150‑1 (MQ=255)
                                                                                                                      tgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCa                                 >  2:363050/1‑150 (MQ=255)
                                                                                                                                  aaCGGCCTCTTCGCGGGCTTGCACGCGACAGAGCCATTACATTCTTATTGTTGTTGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAg                     <  2:520951/150‑1 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCTGCTGCTTCGCCGCGTTACGGCTCCAGTCCAGCAGCACGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                 >  2:15913/1‑150 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                 >  1:310659/1‑150 (MQ=255)
                                                                                                                                        cctttCGCGGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc               >  1:51406/4‑150 (MQ=255)
                                                                                                                                          ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAg                                                                     <  2:466601/94‑1 (MQ=255)
                                                                                                                                          ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAg                                                                     >  1:466601/1‑94 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGACGTTGCTTCGCCGCGTTACGGGTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAac   >  2:29920/1‑150 (MQ=255)
                                                                                                                                                     tGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAaca  <  2:310659/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACA  >  NC_002947/4705311‑4705609

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