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.G80.8% 69.3 / 13.0 26intergenic (+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/4);  new base G (9/12);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01
*NC_0029471,070,2462.A80.8% 67.3 / 13.2 26intergenic (+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/4);  new base A (9/12);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA  >  NC_002947/1070103‑1070379
                                                                                                                                                    ||                                                                                                                                     
cgccgcCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Gag                                                                                                                                    <  1:115009/147‑2 (MQ=255)
cgccgcCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGAGCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Ga                                                                                                                                     <  6:84363/146‑1 (MQ=255)
    gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑Gagatcga                                                                                                                               >  7:8779/1‑142 (MQ=255)
         ttCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGcgcc                                                                                                                            >  2:130864/1‑148 (MQ=255)
              aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCcgcgcg                                                                                                                       >  5:134144/1‑148 (MQ=255)
                    gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCgcgggcg                                                                                                                   <  5:69261/146‑1 (MQ=255)
                    gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCgcggg                                                                                                                     <  5:20207/144‑1 (MQ=255)
                         gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCg                                                                                                            <  6:62512/148‑1 (MQ=255)
                           tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                          <  2:79298/148‑1 (MQ=38)
                                ccGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtctc                                                                                                       >  5:79132/1‑146 (MQ=255)
                                        gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGAt                                                                                              >  8:71786/1‑147 (MQ=255)
                                              tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                        >  6:72502/1‑147 (MQ=255)
                                                         tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGGTTGTTCGTAAGGGGTATGCCATTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                            <  6:15858/148‑1 (MQ=38)
                                                             tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑A‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                                        >  5:57491/1‑148 (MQ=255)
                                                                  aGCTGAAGAAGGCTCGCCAGTCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGAGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcgg                                                                   <  6:108490/148‑1 (MQ=38)
                                                                  aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGGAAGCCTAAAGGCGTGCGTCTggcgg                                                                   >  4:111591/1‑148 (MQ=38)
                                                                        agaagGCTCGCCAGGCCAAGGGCTGATCGGGTGTGTTCGTACGGGGTATGCGGTTAGGTTTTCGGCGCCGC‑‑TGA‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGTGGCCCt                                                               <  5:93566/146‑1 (MQ=37)
                                                                                   cAGGCCAAGGGCTGATGGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGGGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                                  <  2:28461/148‑1 (MQ=38)
                                                                                          aGGGCTGATCGGGTTTGTTCGTAAGTGGTATGCCGTTAGGTTTTCGGCGC‑‑CGGTGA‑‑GAGCGAGCGCCGCGCGGGCGGCGCGCGAGCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                           <  7:29037/148‑1 (MQ=38)
                                                                                                  tCGGGTTTGTGCGTCAGGGGTAGGCCGGTAGGGTTTCGGCGC‑‑C‑‑TGTGAGAGCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAt                                   <  3:71032/148‑1 (MQ=37)
                                                                                                     ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGACCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                <  6:23280/148‑1 (MQ=255)
                                                                                                                aaGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgc                                 >  1:96745/1‑136 (MQ=255)
                                                                                                                 aGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                <  2:96745/136‑1 (MQ=255)
                                                                                                                           cGTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTAACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgt          >  5:124492/1‑148 (MQ=255)
                                                                                                                            gTTAGGTTTTCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg         <  5:115553/148‑1 (MQ=255)
                                                                                                                                   tttCGGCGC‑‑C‑‑TGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa  <  6:134144/148‑1 (MQ=255)
                                                                                                                                                    ||                                                                                                                                     
CGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGC‑‑C‑‑TGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACA  >  NC_002947/1070103‑1070379

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