Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,070,248 +GA intergenic (+48/+94) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,2461.G92.3% 76.8 / 0.5 27intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (14/10);  total (16/11)
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.94e-01
*NC_0029471,070,2462.A92.6% 76.7 / 0.2 27intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (15/10);  total (16/11)
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.10e-01

CGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGT  >  NC_002947/1070106‑1070374
                                                                                                                                               ||                                                                                                                                
cgccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑Gagat                                                                                                                             >  5:26621/1‑143 (MQ=255)
         gCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCcgcg                                                                                                                    >  1:19855/1‑148 (MQ=255)
         gCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCcgcg                                                                                                                    >  4:164623/1‑148 (MQ=255)
            tctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCgcgg                                                                                                                 >  3:6489/1‑148 (MQ=255)
                  tAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGgcggcg                                                                                                           >  3:117435/1‑148 (MQ=255)
                     tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGGCAAGGGCTGATCGGGGTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                        <  5:159774/148‑1 (MQ=39)
                             ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCGCGCCAGGCCAAGGGCTGAGCGGGTTTGGTCGGAAGGGGTAGGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                <  7:149872/148‑1 (MQ=37)
                                          ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAg                                                                                    >  2:185406/1‑147 (MQ=255)
                                          ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATAGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                   >  2:185411/1‑148 (MQ=255)
                                              gctgtgGGTGGCTTTCACGCTGAAGAAGGCGGGCCAGGCCAAGGGCTGAGCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaa                                                                               <  6:113383/148‑1 (MQ=37)
                                                ggggggtggCTGGCAAGCTGAAGACGGCGCGCCAGGCCAAGGGCGGAGCGGGGTTGAGCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCGGGGA‑‑GAGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                              <  3:57533/144‑1 (MQ=21)
                                                      tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                       >  2:115789/1‑148 (MQ=255)
                                                          tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTTAGATCGAGCGCCGCGCGGGCGGCGCGCGAGCGCACAGGCGATGAAACCCTAAAGGCATgaggct                                                                   >  7:166701/1‑143 (MQ=38)
                                                          tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                                   >  2:7498/1‑148 (MQ=255)
                                                             cAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                  <  2:12525/146‑1 (MQ=255)
                                                                      gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                       >  6:161232/1‑148 (MQ=255)
                                                                            tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                 <  4:131208/148‑1 (MQ=255)
                                                                              gCCAGGCCAAGGGCGGATCGGGTTTGTTCGTACGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                               <  3:836/148‑1 (MQ=38)
                                                                                   gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGCCGCGTTCCCCTTcc                                          >  1:13298/1‑148 (MQ=255)
                                                                                    gatctggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGAGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                           <  5:14254/141‑1 (MQ=38)
                                                                                       ggggCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                         <  8:69847/112‑1 (MQ=255)
                                                                                            tGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCACCt                                  >  2:122202/1‑147 (MQ=255)
                                                                                               tCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATAGAGCGCCGCGCGGGCGTCGCTCGATCTCACAGGCGCTGAAAGCCTAAAGGCATTCGTCTGGCGGCCCTGACGCGTTCCCCGTCCATGACCCTCAAt                              >  7:73873/1‑148 (MQ=38)
                                                                                                  ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                           <  2:142525/148‑1 (MQ=255)
                                                                                                               ggggTATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGCCGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc              >  3:64330/1‑148 (MQ=255)
                                                                                                                   tATGCCGTTAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCCCAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGCCGCGGTCCCCTTCCCTGCCCCTCAATGCTGCCTCCGCCAAgcacgc          >  3:133326/1‑148 (MQ=38)
                                                                                                                           tAGGTTTTCGGCGCC‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt  <  3:183431/148‑1 (MQ=255)
                                                                                                                                               ||                                                                                                                                
CGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGT  >  NC_002947/1070106‑1070374

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