Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,312,151 G→C Q281E (CAA→GAA)  alsR ← DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,1510GC95.2% 63.6 / ‑6.0 21Q281E (CAA→GAA) alsRDNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (11/9);  minor base A (0/1);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCC  >  NC_000913/4311908‑4312396
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
cgctccagttaatcaaatcacgcatacgagctctcgcctacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGCGAGAAGAAATCAACAATTCAGCCGGTTTTTTCCatatta                                                                                                                                                                                                                                                >  2:167232‑M2/150‑250 (MQ=255)
                atcacgcctacgcgctctcgactacagtaagcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcaattttgacgtgaGGTAAGGGGACATTAACTCCTCGATGGCGGAATTCCCGTTTTCTTAGTTCAGCGCTCCTGGTGATAAGAAATCAACAATTGCGCCAGTTTTTTCCATATTAATATTAGCTTGTTCAAt                                                                                                                                                                                                                               <  2:349311‑M2/118‑1 (MQ=255)
                      catacgcgctctcgactaccgtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaaaaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTg                                                                                                                                                                                                                           >  2:60217‑M2/128‑249 (MQ=255)
                       atacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggacattccctgatattccgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTCTTTCAGAGCGCCTGGTGAGAAGAAATCAACAATTGAGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTTAGc                                                                                                                                                                                                                        >  1:218129‑M2/127‑251 (MQ=255)
                         tcgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCt                                                                                                                                                                                                                       <  2:190452‑M2/126‑1 (MQ=255)
                                       acagtatccatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgccgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTCTTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTACTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACa                                                                                                                                                                                                         >  1:291118‑M2/111‑250 (MQ=255)
                                                       atgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGATAAAAAAAGCGTCCAGCaa                                                                                                                                                                                        <  1:296551‑M2/157‑1 (MQ=255)
                                                                            cattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGTTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGc                                                                                                                                                                   <  2:130394‑M2/178‑1 (MQ=255)
                                                                               tccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCTGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCg                                                                                                                                                                 <  1:191848‑M2/180‑1 (MQ=255)
                                                                                                   atttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTGCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAATGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACg                                                                                                                                            <  1:218794‑M2/201‑1 (MQ=255)
                                                                                                                                     gtaagttttgaggtgaGGTAAGGGGACACTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCCAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATc                                                                                                          >  2:132046‑M2/17‑251 (MQ=255)
                                                                                                                                         gttttgacgtgaGGTAAGGGGCCAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGCCAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGatcgtcc                                                                                                      >  2:13409‑M2/13‑248 (MQ=255)
                                                                                                                                                cgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAAGATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCg                                                                                               <  1:207685‑M2/246‑1 (MQ=255)
                                                                                                                                                      gTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTGTTTCCATATTAATATTAGCTTGTTCAATGTTGGGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGGAATTTTGTTGCCGATGCATTACGACCTAATAGCGTCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCg                                                                                         <  1:140880/251‑1 (MQ=255)
                                                                                                                                                                     aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGCAGCAATCAACAATTGCGCCGGTTTTTTCCATATTATTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCt                                                                            >  2:153563/1‑250 (MQ=255)
                                                                                                                                                                                              ttcgtATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTTCGCCGGTTTTTTACATATTAAGATTAGCTTGTTCAATGTTGAGCTTGGCGACAGAGACAAAAAAAGCGTCCATCAAGGTTGTTTGTAATATTCTTGCCGAGGAATTACGACCTAATAACGGCGTTTCCGGGGCTTGTTAGCAAATAATATGATCTGCCAGTTTCGCTATCGGTTAGGGGTAGCTATTGTTTATACAAATAAGCTTTGcc                                                   <  2:218129/246‑1 (MQ=255)
                                                                                                                                                                                                   aTTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTc                                             <  1:77328/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggTGAGAAGAACTCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCAGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGttt                               >  1:215279/1‑250 (MQ=255)
                                                                                                                                                                                                                   tGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCt                              >  1:163096/1‑250 (MQ=255)
                                                                                                                                                                                                                         gaaATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGCGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCg                       >  2:41970/1‑251 (MQ=255)
                                                                                                                                                                                                                                              tttttCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTATACGGCCGCTTTTACATCACTGGTTCGcc  >  2:316105/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCC  >  NC_000913/4311908‑4312396

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