Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,455:1 +G 100% intergenic (+589/+167) 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,4551.G93.3% 40.5 / ‑3.4 15intergenic (+589/+167)rtcR/glpGDNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (7/7);  total (7/8)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCA  >  NC_000913/3560255‑3560707
                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
tGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATTAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTTACCGAACTGCAAGGCGGCGCCGCACAGCCATGCTGGCGCAGGATGAGGTTTACCCCCGCAGGCCCATTGCTGTCCGGGCGACAATTCCGCGGCATGACgggtcactgccct                                                                                                                                                                                     >  2:48100/1‑264 (MQ=255)
                             ttcgcTACCGCTTCCGGTGTGGTGCCGATATCGATTAACAGCGTCGAGGCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGACGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                                         <  1:48100/272‑1 (MQ=255)
                             tgcgCGACCGCTGGCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGTGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCAGTCTGGGTGGCCTTGCGCTCGTTGCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                                                                                                                                         <  1:251443/274‑1 (MQ=255)
                              gcgcTACCGCTTCCTGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGTATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGTCTTGCGATCGTGCCACGGCGCGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGCTCTGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                                                                                                                                                <  1:190368/250‑1 (MQ=255)
                                       ctcccGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcg                                                                                                                                                                                                                                                              <  1:127952/160‑2 (MQ=255)
                                                            cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                                                                                                                     <  1:78428/231‑1 (MQ=255)
                                                                 aaCAGCGTCGAGCCATTGGGGAGTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTGTCTTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCTCACCGCCATGATGGCGGAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                                                                                                                                                <  1:154136/215‑1 (MQ=255)
                                                                     gCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCGCTTTTTCTTCGGTCTGGGTGGCCTTGAGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCGGACGGAGGAATGCTCTACCAGCTCTTCGGTACCGACATAACCCTGCTTTTTTACCGGTTCGCTGTCACCGTtggttggtggtgtt                                                                                                                >  1:61400/1‑263 (MQ=255)
                                                                        tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGAGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                                                                                                                                         >  1:262782/1‑215 (MQ=255)
                                                                                                                gcgcTTTTTTTGTCCGTGCTGGGTGGCGTTGGGATCGGGCTACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGACCCGCCATGATGGGGCAGGATCAAGTTTTGCTCCGCCAGCTCATTTATGTCGCTGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTGATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttcg                                                                         <  2:61400/271‑1 (MQ=255)
                                                                                                                  gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTGTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACTGAGAAATGCTCTACCAGCTCtt                                                                                                                                                                           <  1:2582/171‑1 (MQ=255)
                                                                                                                                               cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGGTACGTTGTGTTTGTTTCATTTATAAATTCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAAACATCCCCaaaa                                      >  2:126431/1‑274 (MQ=255)
                                                                                                                                                                         aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTTAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTCCCAGACCTTCGTTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCTTTTATAAATACCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAACGCCATCGCTAAACCCACGGCCAACCCGCCGATGTGTGc             >  1:225210/1‑274 (MQ=255)
                                                                                                                                                                             gAAGCCAGCGCCCCACCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                                                                                                                                            <  1:74262/143‑1 (MQ=255)
                                                                                                                                                                                   aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTCAACCCACGGCAAACCCGGCGACGTGTGCTCCGTTCGccc  >  1:167362/1‑274 (MQ=255)
                                                                                                                                                                                     cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACaaa                                                    >  1:51651/1‑223 (MQ=255)
                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
TGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCCGTTCGCCA  >  NC_000913/3560255‑3560707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: