Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,506,593 G→A 14.8% intergenic (+52/‑147) PP_4839 → / → cycA iron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,506,5930GA14.8% 59.8 / 4.0 27intergenic (+52/‑147)PP_4839/cycAiron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base G (16/7);  new base A (2/2);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTAT  >  NC_002947/5506469‑5506741
                                                                                                                            |                                                                                                                                                    
ggTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCCCGTTGTAGGAGCGGGTTTa                                                                                                                              >  4:314983/1‑149 (MQ=255)
                            gATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCCCACCAACGACGGACCGCGTTGTAGGATCGGGTTTACCCCCGAAGAGGCCGGTACAGGCGAAAt                                                                                                  >  8:41302/1‑149 (MQ=38)
                                  gtggtggTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGGCTTGTGTTGCCTGTGCTGGCCGCTTCGCGGGTAACCCCGCTCGTCCGAGGGACGGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGa                                                                                             <  3:314983/147‑1 (MQ=38)
                                               gccgccGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACcgcg                                                                                                                                               >  6:272364/1‑85 (MQ=255)
                                               gccgccGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACcgcg                                                                                                                                               <  5:272364/85‑1 (MQ=255)
                                                            ggCCAGCGCTTGAGCCTTGTGTTGCCGGTGCTGGCCTCTTCGCGGGTAACCCCGCTCCTACGAGGGACCGCGTTGTAGCAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGc                                                                                                <  7:281192/119‑1 (MQ=255)
                                                                       gAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTg                                                       >  8:335918/1‑149 (MQ=255)
                                                                              gtgtTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcg                                                >  7:275003/1‑149 (MQ=255)
                                                                                          cTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTtagcct                                     <  7:335918/148‑1 (MQ=255)
                                                                                              ccTCTTCGCGGGTAAACCCGCTCCTCCGAGGGCCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGcc                                                               <  4:90765/118‑1 (MQ=255)
                                                                                              ccTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGGGTTGTAGGAGCGGGTTTACCCGCGAAGCGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCTGAGTCGATCTGTCGCGGCCCGCTCTGGCACGCGCGTGTTAGCCTAGGCg                                >  1:259248/1‑149 (MQ=38)
                                                                                              ccTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGcc                                                               >  3:90765/1‑118 (MQ=255)
                                                                                               ctctTCGCGGGTAAACCCGCTCCTACGAGGGCCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTgg                                                      >  5:336524/1‑126 (MQ=255)
                                                                                                  ttCGCGGGTAAACCCGCTCCTACGAGGGGCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTa                            >  4:332329/1‑149 (MQ=255)
                                                                                                          tAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGa                     >  2:67351/1‑148 (MQ=255)
                                                                                                           aaaCCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCGTCGCCGGCTCCCTCTGAc                    >  6:244817/1‑148 (MQ=255)
                                                                                                             aCCCGCTCCTACGAGAGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGAc                    >  3:111599/1‑146 (MQ=255)
                                                                                                             aCCCGCTCCTACGAGAGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGAc                    >  3:106245/1‑146 (MQ=255)
                                                                                                               ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaag                 >  4:137901/1‑147 (MQ=255)
                                                                                                                 gCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAagaaga             >  3:101736/1‑149 (MQ=255)
                                                                                                                   tCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcgc                                               <  5:169434/113‑1 (MQ=255)
                                                                                                                   tCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcgc                                               >  6:169434/1‑113 (MQ=255)
                                                                                                                     cTACGAGGGACCGCGTTGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGc         <  7:41302/149‑1 (MQ=255)
                                                                                                                     cTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTg          >  4:182665/1‑148 (MQ=255)
                                                                                                                       aCGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGTCTACCTCTGACAAGAAGACTGCAg       >  7:205273/1‑149 (MQ=255)
                                                                                                                        cGAGAGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCa        <  4:111599/147‑1 (MQ=255)
                                                                                                                        cGAGAGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCa        <  4:106245/147‑1 (MQ=255)
                                                                                                                           gggACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGACTACCTCTGACAAGAAGACTGCAGCCTa   >  2:64759/1‑149 (MQ=255)
                                                                                                                            ggACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCtatgggtgaggtttgcgtcataggctgcagtcttcttgtcagaggtagccggctaggctaacac  >  2:195769‑M1/1‑85 (MQ=255)
                                                                                                                            |                                                                                                                                                    
GGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTAT  >  NC_002947/5506469‑5506741

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: