Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 313,543 G→A intergenic (‑301/‑1748) 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,5430GA94.1% 50.6 / ‑5.3 18intergenic (‑301/‑1748)ykgR/insE1putative membrane protein YkgR/IS3 element protein InsE
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (6/10);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCG  >  NC_000913/313310‑313780
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             
aaCCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTaaa                                                                                                                                                                                                                              <  2:293448/251‑1 (MQ=255)
            gACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATg                                                                                                                                                                                                                  <  2:99768/251‑1 (MQ=255)
                   aTCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCt                                                                                                                                                                                                                                    >  2:225830/1‑226 (MQ=255)
                   aTCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCt                                                                                                                                                                                                                                    <  1:225830/226‑1 (MQ=255)
                                 tAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTac                                                                                                                                                                                              >  1:50120/1‑250 (MQ=255)
                                  aGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAAGGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATAGCGCATCGCATCATATTGTCGTATCGCCTGGGGTGGTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTGCTAGGAAATCGTCAACTATGCCGGATACATACTGGATTaca                                                                                                                                                                                             <  1:315854/250‑1 (MQ=255)
                                                    aTCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTTAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCCGGCACGCCTTAt                                                                                                                                                                            <  2:70549/249‑1 (MQ=255)
                                                                                            aGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAAt                                                                                                                                  <  2:50120/251‑1 (MQ=255)
                                                                                                                  gATCTCAAGACTTGCCCAGTCCGGGCGTTGATGGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGTAAGGGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTAc                                                                                                             <  2:62381/250‑1 (MQ=255)
                                                                                                                  gATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACTGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTAc                                                                                                             <  1:267818/250‑1 (MQ=255)
                                                                                                                                    ggccGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGGCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGATAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGTTTATTCATTGATTTGATAGtttt                                                                                          <  2:154081/249‑1 (MQ=255)
                                                                                                                                           cGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCg                                                                                   >  2:291747/1‑251 (MQ=255)
                                                                                                                                            gTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCgg                                                                                  >  1:291633/1‑251 (MQ=255)
                                                                                                                                                           ggagCATATTGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATAGCGGTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCATTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTCTTACTCTCGGGAGAATAATAGatat                                                                   <  2:12945/251‑1 (MQ=255)
                                                                                                                                                                                   tGTCGTATGTCCTGGGGTGTTATAAAAGATATCGTTGTTGGTGACCTTGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGAAGCAGAGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATAGGTTACTCATTGATTTGCTCGGTTTACTCTCGGGAGAATAATAGATATTTAATCTCTTAAGGGAAAccagcc                                           >  2:293747/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttgGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGAc                                                                                                                                                       >  2:108551/1‑106 (MQ=255)
                                                                                                                                                                                                                        ttgGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGAc                                                                                                                                                       <  1:108551/106‑1 (MQ=255)
                                                                                                                                                                                                                            tGACCTGGGAAAGGAATTGAGATCTATTACACCGTCAACTATGCCGGATAGATACTGGATCACACTGCAGGCACGCCTTATGCGAGAACGTGCCGCAGTGGCGGGATAATTATCTGAAAGCATTTGTGCGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGCATCATAGATATTTAATCCATTAACGGAAACCAGCACGTTCCTTTCGATGCCTGATTTTGATCCCATAGTTGATCCg  >  2:96656/1‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             
AACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCG  >  NC_000913/313310‑313780

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