Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 353,067 A→G 100% intergenic (+34/+84) hisF → / ← cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,0670AG99.6% 128.2 / ‑7.1 40intergenic (+34/+84)hisF/cbcVimidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/19);  minor base C (1/0);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352925‑353209
                                                                                                                                              |                                                                                                                                               
ggAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                           >  8:8380/1‑149 (MQ=255)
   aGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGccc                                                                                                                                        >  1:51069/1‑149 (MQ=255)
    gggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCg                                                                                                                                       >  1:204953/1‑149 (MQ=255)
           ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                                 >  7:65661/1‑148 (MQ=255)
                  ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTg                                                                                                                        >  8:62182/1‑149 (MQ=255)
                  ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTg                                                                                                                        <  7:91230/149‑1 (MQ=255)
                    tGCTGGGGGGCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCGAAGGCCTACATGGCTTCGCGGGGGATCGTCGTTCGGTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCt                                                                                                                      <  6:201053/149‑1 (MQ=38)
                      cTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTga                                                                                                                    <  7:223807/149‑1 (MQ=255)
                       tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCt                                                                                                                      >  6:52982/1‑146 (MQ=255)
                         gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCt                                                                                                                      >  8:109656/1‑144 (MQ=255)
                              cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCtg                                                                                                            >  5:230452/1‑149 (MQ=255)
                               aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCtgt                                                                                                           <  2:204953/149‑1 (MQ=255)
                               aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCtgt                                                                                                           <  2:51069/149‑1 (MQ=255)
                               aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCtgt                                                                                                           <  6:259748/149‑1 (MQ=255)
                                              ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGc                                                                                             <  7:109656/148‑1 (MQ=255)
                                                  cgTACACCGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCCCTGATCCATTCCCCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGACATCTCTGTGGGAGCGGGCGTGCCcgcg                                                                                        <  3:22327/148‑1 (MQ=38)
                                                   gTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGa                                                                                       <  8:65661/149‑1 (MQ=255)
                                                      cacaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                     <  8:125921/148‑1 (MQ=255)
                                                       acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                   >  1:171104/1‑149 (MQ=255)
                                                                gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGc                                                                          >  3:166563/1‑149 (MQ=255)
                                                                  ggCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgg                                                                        <  5:52982/149‑1 (MQ=255)
                                                                   gCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                                       <  6:230452/149‑1 (MQ=255)
                                                                     cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATACAATGACCGGGGCTCCTTTGCAGCCCATCGCCGGCACGCCCGCACCCCCAG‑GTTTGATAACATCACTGTGGGAGCCGGCGTGCCCCCCAAGGGCCGCCAAGCGGccc                                                                     >  5:241838/1‑149 (MQ=255)
                                                                      aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCGTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGcccc                                                                    <  7:62182/149‑1 (MQ=255)
                                                                                tGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTc                                                           >  6:117138/1‑148 (MQ=255)
                                                                                  gCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCaa                                                         <  7:8380/148‑1 (MQ=255)
                                                                                   cTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCCGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAAcc                                                       >  1:199499/1‑149 (MQ=255)
                                                                                       gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                   <  2:171104/149‑1 (MQ=255)
                                                                                       gcgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTgg                                                   <  2:199499/149‑1 (MQ=255)
                                                                                                    gtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                                      <  5:117138/149‑1 (MQ=255)
                                                                                                             tCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCACTGTGGGAGCGGGCGTGCCCGCGAACGGCCGCCAACCGGCCCCAAATGCATCAACCGTGGTTCTTGCCCAGCAGCGCCTgg                              >  4:154222/1‑148 (MQ=37)
                                                                                                                aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTTCTTAAAAAGTGGTTCTTGCCGAGCAGCGCGTGGTaa                           >  6:266654/1‑148 (MQ=38)
                                                                                                                aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCACCCACAG‑GTTTGCTGAGATCTATGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAATGCTTCAACCGTGGTTCATGCCGAGCAGCCCGTGGTa                            >  6:270180/1‑147 (MQ=38)
                                                                                                                    ccgggggCGGCTTTGCCGCCCTGCGCGGGCACCCCCGCTCCCACCG‑GTTGGCTCACAGCTCTGGGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTc                      <  1:128684/146‑1 (MQ=37)
                                                                                                                       ggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGt                   <  4:166563/149‑1 (MQ=255)
                                                                                                                       ggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCAACCACAC‑GGTTGCTGAGAACGCTGTGGGAGCGGGCGGGCCCCCGAAGGGCCGCGCAGCCGCCCCAAATGCTTCAACCCTGGTTCTTGCCGACCAGCGCCTGGTAAAGCTccg                    >  4:95508/1‑146 (MQ=255)
                                                                                                                            tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATCGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATccc              >  2:53848/1‑149 (MQ=255)
                                                                                                                               tttGCAGCCCTTCGCGGGCACGCCCCCACCCACAG‑GTTTGCTGAGAAATCAGTGGGAGCGGGCGTGCCCGCGAAGGGCCCCCAAGCCGCCCCAAATTCATCAAACGTGGGTATTGCCGAGCAGCGCGTGGTAAAACACGGTATccccc            >  4:64878/1‑147 (MQ=255)
                                                                                                                                        cTTCGCGGGCACGCCCGCTCCCACAGAGATCTCAGCCAAC‑CTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCGTCAACCGTGGTTGTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGGATccccc  >  7:152093/1‑148 (MQ=18)
                                                                                                                                        cTTCGCGGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTAGGTATCCCCGAGTATCCcca  >  6:86043/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAG‑GTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCA  >  NC_002947/352925‑353209

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: ‑