Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,362 G→A intergenic (+496/+260) rtcR → / ← glpG DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,3620GA91.3% 54.6 / ‑6.1 23intergenic (+496/+260)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/8);  minor base C (0/2);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

AGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTC  >  NC_000913/3560126‑3560603
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
agaTAAAAGCGAGCGCCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCtttt                                                                                                                                                                                                                                      <  2:47292/251‑1 (MQ=255)
                 gCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTAGCCACCGGAGACAATGATGCCAAAATCTTCTTTTACAATCAACGTGTTAACCACGTTG‑AGCTTGTTGGTGACAATGCGCAAATTACTGTGATTGAGCAGAGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑CGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGc                                                                                                                                                                                                                     >  1:219413/1‑251 (MQ=255)
                        cAATGATCCCGCCATCGCAGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCgg                                                                                                                                                                                                              >  2:182426/1‑250 (MQ=255)
                          aTGATCCCGCCATCGCGGCTGCGTAATTCGCCCCCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATc                                                                                                                                                                                                            <  2:151220/251‑1 (MQ=255)
                           tGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCg                                                                                                                                                                                                           >  1:55177/1‑251 (MQ=255)
                                         cgGCTTCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCAGCAGCGCTTTAGCAAC‑TTGTAGTTTGTTGTTGACAATGCGCAAATTGCTGTGATTGAGCAGTTCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGT‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGGTCGTGCCACTGCGTGtt                                                                                                                                                                                             <  2:120500/251‑1 (MQ=255)
                                           gCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCCGGGCGTGCGCTACCGCTTCCGGCGTGGTGCAGATATCGATAAACAGCG‑TAG‑AGCCATTGGGG‑ATTTG‑CTCCGACACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCACCGGCGTGTTa                                                                                                                                                                                            >  2:176835/1‑250 (MQ=255)
                                                                      gatgCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCACTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGc                                                                                                                                                                >  1:80474/1‑251 (MQ=255)
                                                                                 cttcttTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                                                     >  1:231790/1‑249 (MQ=255)
                                                                                        tACCATCAACGTGTTAGCAACGTGG‑AGATTGTTGGTGACAATGCGCAAATTGCTGGGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                               <  1:105538/250‑1 (MQ=255)
                                                                                          ccATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                                                                                                            >  1:8961/1‑251 (MQ=255)
                                                                                           catcacagggttaaaaccGTTG‑GGATTTTTGGTAACAATGAGAAAATTGCGGTGATGGAGCAGGGCGGCCCCTCCCGCTTCCGGCGGGGGGCCGATATCGAAAAACGGCG‑GCG‑‑GCCCTTGGGGGATTGC‑CTCCGCCACTTTGCGGCCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGt                                                                                                                                           <  2:75836/234‑1 (MQ=255)
                                                                                                                          ggTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCgcggc                                                                                                             >  1:68273/1‑251 (MQ=255)
                                                                                                                                       aaaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                >  1:101408/1‑251 (MQ=255)
                                                                                                                                        aaTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                               >  1:20136/1‑251 (MQ=255)
                                                                                                                                                   ttGAGCAGGGCGTGCGCTCCCGCTTCCGTCGTGGTGCCCATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCATCGCCGCACCGCCATGATCGCGCAGGATCAGGCTTTGCGCCGCCAGCTCATTGAGGTCgcgggtaatcacctgcggggtgccgggga                                                                                     >  2:6209/1‑225 (MQ=255)
                                                                                                                                                            gcgtccgCTCCCGCTTCCGGCGGGGTGCCGATTTCTA‑AAACCGGGCTCG‑CGCCATGGGGG‑TTTTC‑CTCCGCCACTTGGCGGACTATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCCCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                           <  2:151143/251‑1 (MQ=255)
                                                                                                                                                             gtggcgttACCGCTTCCGGCGTGGTGCTGATTTCGAAAACCTGC‑‑TCGGACCTATTGGGG‑‑TTTGTCTCCGCCACTTTGCGGACTATGCGCTCTGTTGCTCCGGTCTGGGGGGCTTTGCGATCGTGCAACGGCGTTTTAACAGACCTGGAAGCCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGTCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                           <  2:151130/243‑1 (MQ=255)
                                                                                                                                                                                cGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                       >  1:235030/1‑251 (MQ=255)
                                                                                                                                                                                       cctatatAGATAAACAGGG‑TCG‑GGCAATTGGGT‑ATTTC‑CCCCGCCACTTTGCGGACGATGCGCTCGTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGGGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                                <  2:235030/248‑1 (MQ=255)
                                                                                                                                                                                                             g‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                            >  2:165923/1‑250 (MQ=255)
                                                                                                                                                                                                                                ccTCCGCCACTTTGCGGCCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt           <  1:229122/250‑1 (MQ=255)
                                                                                                                                                                                                                                          ttggCGGACGATGCCTTCTTTTTCTTGGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGAAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc  <  2:180736/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
AGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTG‑AGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCG‑TCG‑AGCCATTGGGG‑ATTTG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTC  >  NC_000913/3560126‑3560603

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑