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,5430GA96.2% 80.6 / NA 26intergenic (‑301/‑1748)ykgR/insE1putative membrane protein YkgR/IS3 element protein InsE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/12);  minor base C (1/0);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTT  >  NC_000913/313325‑313788
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
aTGTATCACAAATGTGAATAGATGTATTACATCAACTATCTTTGATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCg                                                                                                                                                                                                                       <  2:87713/251‑1 (MQ=255)
                   aGATGTCTTACATCAACTATCATTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCAGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTacac                                                                                                                                                                                                    >  2:58037/1‑251 (MQ=255)
                         aTTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAgg                                                                                                                                                                                              <  2:176502/251‑1 (MQ=255)
                               tCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGcc                                                                                                                                                                                        >  2:306473/1‑251 (MQ=255)
                                 aaCTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCtt                                                                                                                                                                                      >  2:81088/1‑251 (MQ=255)
                                     aTCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATg                                                                                                                                                                                   >  1:436950/1‑250 (MQ=255)
                                                         cATTGATATATGTCGCCTGACGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGATGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACg                                                                                                                                                              >  1:156796/1‑251 (MQ=255)
                                                          aTTGATATATGGCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTTGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGGCCTGGGAGAGAAATGGAGTTCTATTAATCCGTCAGCTATGCCGGATGCAGACTGGATGACACTGCAGGCACGCCTGAGGAGAGAACGTGCCGCAGTGACgg                                                                                                                                                             <  2:190358/251‑1 (MQ=255)
                                                                                    tCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTtgtg                                                                                                                                     >  1:179422/1‑249 (MQ=255)
                                                                                                         aaGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCAttga                                                                                                              <  1:29956/251‑1 (MQ=255)
                                                                                                          aGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCAttgat                                                                                                             >  1:26459/1‑251 (MQ=255)
                                                                                                                    aGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGttt                                                                                                   <  1:93409/251‑1 (MQ=255)
                                                                                                                                    ggTGCTGAGGAGCCTATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGaat                                                                                    >  2:299485/1‑250 (MQ=255)
                                                                                                                                    ggTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGaat                                                                                    <  2:402960/250‑1 (MQ=255)
                                                                                                                                    ggTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATAGCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGGGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGaata                                                                                   <  1:393328/251‑1 (MQ=255)
                                                                                                                                          gaggagCATCTCGCATCTAATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAAAAGTCAACTATGACGGATAAATACTGGATTACACTGCAGGCAAGCCTTATGAGAGAACGTGCCGCAGCGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGAATCGGTTACTCATTGATTTGATAGTTTTACTATCGGGAGaataat                                                                                 >  1:57566/1‑247 (MQ=255)
                                                                                                                                           aggagCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTAGTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATGATCTGAAAGAATTTGTGAGGCTGTATCTGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGata                                                                            <  1:58037/251‑1 (MQ=255)
                                                                                                                                                cATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTAc                                                                                                                     <  1:292314/205‑1 (MQ=255)
                                                                                                                                                cATATCGCATCTCATCATAATGTCGTATCTCCTAGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTa                                                                        >  1:116548/1‑250 (MQ=255)
                                                                                                                                                cATATCGCATCTCATCATAATGTCGTATCTCATGGGGTGTTATACAAGATATCGTTGTTGGAGCCCTGGGAGAGACATTGAGTTCTCTTAAACCGTCAACTATGCCGGATACATACTGGATTACAATGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGCGGCTGTATCGGTTAc                                                                                                                     >  2:292314/1‑205 (MQ=255)
                                                                                                                                                                       cGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGt                                                >  1:141241/1‑251 (MQ=255)
                                                                                                                                                                              ccTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCt                                            <  1:81088/248‑1 (MQ=255)
                                                                                                                                                                              ccTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTc                                         <  2:436950/251‑1 (MQ=255)
                                                                                                                                                                                          taCAAGCTATCGTTGTTGGTGACCTGGGAGCGCAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTAAATAATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATACTAGATATTTAATCCATTAACGGAAACCAGACAGTTCATTTTGTTGCATGAAtt                             >  2:407958/1‑251 (MQ=255)
                                                                                                                                                                                           aCAAGATATCGTTGTTGGTGACCTGGGAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAAttt                            >  1:96243/1‑251 (MQ=255)
                                                                                                                                                                                                        gttgGTGACCTGGTAGAGAAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTGATGAGAGAACGTTCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATGTGATCCCATAGttt               <  2:403120/251‑1 (MQ=255)
                                                                                                                                                                                                                     gagagaAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGtt  >  1:179023/1‑251 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
ATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAACGTCATTGATATATGTCGCCTGAAGTCAGTTCCGGGAATGAGTCTGATCTCAAGACTGGCCCAGTCCGGGCGTTGATTGGTGCTGAGGAGCATATCGCATCTCATCATAATGTCGTATCTCCTGGGGTGTTATACAAGATATCGTTGTTGGTGACCTGGGAGAGGAATTGAGTTCTATTAAACCGTCAACTATGCCGGATACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTT  >  NC_000913/313325‑313788

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