Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,313,027 T→G intergenic (‑36/‑323) alsR ← / → rpiB DNA‑binding transcriptional repressor AlsR/allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,0270TG94.7% 61.3 / NA 19intergenic (‑36/‑323)alsR/rpiBDNA‑binding transcriptional repressor AlsR/allose‑6‑phosphate isomerase/ribose‑5‑phosphate isomerase B
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (12/6);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAAC  >  NC_000913/4312803‑4313240
                                                                                                                                                                                                                                |                                                                                                                                                                                                                     
catcGCTTCAGATACCGCCAGAGCTTCTACGACAGCTTTAATTGCGGGTGCTCAACCCAGGTTACCGGTTGTCAGTAACCACTCCACGATGCGTGGTTCATTTTCGGTCATGCCTGCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAATCGCTGAATCAAACTCTTACTGTCTCATCACGTTCTCATGTCCTTGAATTGAAGACGCCACTAGATCTCCACAGCAAGAACAGAAAAGAt                                                                                                                                                                                             <  2:215279/251‑1 (MQ=255)
 atcGCGTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATg                                                                                                                                                                                            >  2:66894/1‑251 (MQ=255)
       tCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGAGCACCACAGCAAGAACAGAAAAGATGCTGGc                                                                                                                                                                                       <  2:37142/250‑1 (MQ=255)
         aGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGAc                                                                                                                                                                                    >  2:38550/1‑251 (MQ=255)
             aCCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTc                                                                                                                                                                                 <  1:121439/250‑1 (MQ=255)
                                        aTTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGtttat                                                                                                                                                     >  2:211642/1‑251 (MQ=255)
                                                                                 cTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACACAAAAATAAAATTTATAAAGATTATTTTATTGTGAACAATGGCGAGGGGCTTAAAATTTTATTTct                                                                                                            >  2:13972/1‑250 (MQ=255)
                                                                                        aTGCGGGTTTCATTTTCTGTCATTCCTTCCTTCTTCATTCGCAGGGAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAATAACAGAAAAGATGCTGGCGACTTCCCAAAAAAATAACATTTATAAAGTTTATTTTATTTTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTactt                                                                                                     >  2:242948/1‑248 (MQ=255)
                                                                                         agCGGTTTTCATTATATGTCATTACTACCTGATTCATTCCCCGGTAAGGCGCTAACACTATACCTTTCGGAACGGCTGAATCAAACTCTGACAGGCTCATCACGTGCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAGGAATAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATATATGTTATTGTATTTTGAACAATGGCGAGTGGCTGAAAATTTTATGTCATTTAATTa                                                                                                    <  2:144782/250‑1 (MQ=255)
                                                                                                     tttCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACat                                                                                        >  1:277269/1‑251 (MQ=255)
                                                                                                                                  aGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGtttttt                                                              >  1:131243/1‑248 (MQ=255)
                                                                                                                                           gcTAACCCTATACCGTTCGGAAGCGCTGACTCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCAACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATg                                                  >  1:296016/1‑251 (MQ=255)
                                                                                                                                                          ttCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTc                                                                                                             <  2:230574/177‑1 (MQ=255)
                                                                                                                                                          ttCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTc                                                                                                             >  1:230574/1‑177 (MQ=255)
                                                                                                                                                           tCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCCCTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAgt                                  >  1:175701/1‑251 (MQ=255)
                                                                                                                                                                 gcgcTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTTCTTGAACTGACGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGCGTGGCTTAACATTTTCTTTCATTTAATTAGTTTCAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGAGAGTGTCttt                            >  1:339163/1‑251 (MQ=255)
                                                                                                                                                                      gAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGGGTGTCTTTAACGc                       >  2:286244/1‑251 (MQ=255)
                                                                                                                                                                                      ggCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGAt       <  1:38550/251‑1 (MQ=255)
                                                                                                                                                                                           atcaCGTTCTCTTGTCCTTGAATTGAAGACGCCACTAGATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAAGTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTGTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCGGATGTAAc  <  2:129673/251‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                     
CATCGCTTCAGATACCGCCAGAGCTTCTGCGACATCTTTAATTGCGGGTGCACAACTCAGGTTACCGGGTTTGAGTAACCACTCCACGATGCGGCTTTCATTTTCTGTCATTCCTTCCTGCTTCATTCGCAGGTAAGGCGCTAACCCTATACCGTTCGGAAGCGCTGAATCAAACTCTGACTGGCTCATCACGTTCTCTTGTCCTTGAATTGAAGACGCCACTATATCACCACAGCAAGAACAGAAAAGATGCTGGCGACTTCACAAAAAAATAAAATTTATAAAGTTTATTTTATTGTGAACAATGGCGAGTGGCTTAAAATTTTATTTCATTTAATTAGTTTAAAAACATGATGTTATTTGTTTTTCTGGTTTTTTAGATGATAGATGTTGCTGGATGGTGAGTGTCTTTAACGCTATTGCGCTGTCAGATGTAAC  >  NC_000913/4312803‑4313240

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: ‑