Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 313,745 C→T intergenic (‑503/‑1546) ykgR ← / → insE1 putative membrane protein YkgR/IS3 element protein InsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913313,7450CT91.7% 19.9 / NA 12intergenic (‑503/‑1546)ykgR/insE1putative membrane protein YkgR/IS3 element protein InsE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/5);  minor base G (0/1);  total (6/6)
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

TCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGGGGTTATGAACGAAATTCATAAATCTGTTTGAGTGTTGCGATGGGTAGTGCAAGTTCGAT  >  NC_000913/313552‑313929
                                                                                                                                                                                                 |                                                                                                                                                                                        
tCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGTCTGTATCGGTTCCTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAATCCAGCCAGTTCTTTTCGATGCCTGAATTTGATCCCATAGTGTATCCGGAATTGTTAATTCATGAACaa                                                                                                                                  <  1:81671/250‑1 (MQ=255)
                   ngCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCTTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATg                                                                                                                >  1:26113/2‑249 (MQ=255)
                                  ggATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGCAGGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTAGGGAGAATAATCGGTATTTAATCCATTAACGGCAACCAGCCAGTTCTTTTCGATGACTGAATTTGATCCCATAGTTTATCCGGAATAGTTAATTCATGAACAATTTGTGTTATCGTTCATGAGAAGCATAACGTaa                                                                                                 >  1:33709/1‑249 (MQ=255)
                                            ncAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCTTTTCGATGACTGAATTTGATCCCATAGTTTATcc                                                                                                                                                        >  1:159550/2‑184 (MQ=255)
                                                                   gTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCTTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAg                                                              <  1:83330/251‑1 (MQ=255)
                                                                                     anNNNNCTGNAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCTTTTCGATGCCTGAATTTGATCCCATAGTTTATcc                                                                                                                                                        <  2:159550/143‑1 (MQ=255)
                                                                                                                               gatttgatAGTTTTACTCTCGGGAGAATAATAGATATTTAAGCCATTAACGGAAACCAGCCAGGTCTTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCTTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAAGTTGTCAGGGGTTATGAACGAAATTCATAAATCTGTTTTAGTGTTGCGATTTGTAGTGCAAGGTCGan  <  1:64321/251‑2 (MQ=255)
                                                                                                                                atttgatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCTTGTCGATGCCTGAATTTGATCCCATAGTTTATTCGGAATGGTTAATTCATGAACAAGTTGTTTTATTGTTCAAGAGAAGCATAACGCAAAGGGAAAAGCGCGATTAGACGGCAGAAGTTGTCAGGGGTTATGAAc                                                   <  1:234071/201‑1 (MQ=255)
                                                                                                                                    gatAGTTTTACTCTACGGAGTATAATCGATATTTAATCCATTAACGGAAACCGGCCAGTTCTTTTCGCTGCCTGCATTTGCTCCCATCGTTTTTCCGGCATGGTTAATTCATGAACAAGGTGTCTTCTCGTTAATGAGaa                                                                                                            >  1:34886/1‑140 (MQ=255)
                                                                                                                                                                                cGGAANNCANNNAGNNCTTTTCGATGCCTGAATTTGATCNCATAGTTTATCCGGAATGGTTAATTCATGNANANGTTGNNTTANCGTTCATGAGNAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGGNGNTATGAACGAAATTCATAAATCTGTTTGAg                             >  2:188450/1‑175 (MQ=255)
                                                                                                                                                                                cGGAAACCAGCCAGTTCGTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAAGGGTTAAGTCATGAACATGTGGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGGGGTTATGAACGAAATTCATAAATCTGTTTGAg                             <  1:188450/175‑1 (MQ=255)
                                                                                                                                                                                      ncagccagTTCTTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTACTTCATGAACACGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGGGGTTATGAACGAAATTCATAAATCTGTTTGAGTGTTGCGATGGGTAGTGCacg        >  1:30316/2‑188 (MQ=255)
                                                                                                                                                                                                 |                                                                                                                                                                                        
TCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGGGGTTATGAACGAAATTCATAAATCTGTTTGAGTGTTGCGATGGGTAGTGCAAGTTCGAT  >  NC_000913/313552‑313929

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