Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,246:1 +GA 100% intergenic (+46/+96) 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.G94.7% 54.1 / ‑2.7 19intergenic (+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/0);  new base G (6/12);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0029471,070,2462.A94.7% 52.3 / ‑2.7 19intergenic (+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/0);  new base A (6/12);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACA  >  NC_002947/1070099‑1070385
                                                                                                                                                     ||                                                                                                                                           
cTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCtgt                                                                                                                                               >  5:360606/1‑148 (MQ=255)
        gccgATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCATCGCCTGT‑‑Gagatcgag                                                                                                                                    >  5:180035/1‑142 (MQ=255)
                   tctcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGACTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcggg                                                                                                                           >  7:135866/1‑149 (MQ=255)
                       ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGACTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcggg                                                                                                                           <  8:135866/145‑1 (MQ=255)
                           aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTCAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                    <  1:115452/149‑1 (MQ=38)
                           aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                    <  4:372215/148‑1 (MQ=39)
                           aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGCATGCCTTTAGGCTTTCATCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                    <  2:304199/148‑1 (MQ=25)
                            tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                  <  6:360606/149‑1 (MQ=38)
                              cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                 <  5:287242/148‑1 (MQ=255)
                              cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCATCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                <  6:180035/149‑1 (MQ=37)
                              cTGATCCCGGCATGGATCGCTGTGCTGGGGGTGGCTTTCAAGCTGAAGACGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                <  7:185734/149‑1 (MQ=38)
                                     cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                         <  8:108842/149‑1 (MQ=255)
                                                                 tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                <  8:313891/114‑1 (MQ=37)
                                                                 tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                >  7:313891/1‑114 (MQ=37)
                                                                                 gCTCGCCAGGCCAAGGGCTGCTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                    <  1:22318/94‑1 (MQ=25)
                                                                                 gCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                    >  2:22318/1‑94 (MQ=25)
                                                                                                     aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggcg                                                                                                                         >  7:378842/1‑69 (MQ=21)
                                                                                                     aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCgcgggcg                                                                                                                         <  8:378842/69‑1 (MQ=21)
                                                                                                              gTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                                >  5:210037/1‑149 (MQ=255)
                                                                                                                                             gcgcCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAaca  >  3:58317/1‑149 (MQ=255)
                                                                                                                                                     ||                                                                                                                                           
CTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACA  >  NC_002947/1070099‑1070385

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: