Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,571,280 C→T E99K (GAG→AAG)  yjiR ← fused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,2800CT95.8% 74.8 / NA 24E99K (GAG→AAG) yjiRfused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/17);  minor base G (0/1);  total (6/18)
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

ATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTT  >  NC_000913/4571033‑4571502
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                              
accgggttatggcggcggtgggTGAGGATGAGGCCATCGGGGGGGCCCCCCGGGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCGCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAAAACGGTTGGCAGATTATGCT‑GCACCTCCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAg                                                                                                                                                                                                                             <  2:144274/231‑1 (MQ=255)
             ggccggtggtgttgattagtgagccggcggTAACCAGCGGGCCGTTGTCGACCAGCAGGAGGGCACTCTGCCCTCTCAATACTCGCTGTCCGGCTAACAATTAATAACCGAGAACGTTTTGCAGATTTTGCTAGC‑CCCCCCAGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAAATGTCGCTATGCGCCTTCAGCATATCCAGCAcc                                                                                                                                                                                                                <  2:7772/227‑1 (MQ=255)
                 gCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGGCGCTGTCGGGCATCACGCGGGGAAACTGCTCTCGCCAAACTCGCTGGCCGGCTAACAAGGCATAAACAAGAAAGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGGGTTTC‑GGCTGGGCGCTTCGACATCCGGGGGGCTTTTGCTTTACGGGAACA‑TGGAACAGTCGCTATGCGCCTTCAGCATCTCCAGCACCTGAt                                                                                                                                                                                                            >  1:272389/1‑251 (MQ=255)
                  ggggTGAGGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGGGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATc                                                                                                                                                                                                           <  1:19652/249‑1 (MQ=255)
                                             accaccGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTg                                                                                                                                                                                >  2:267922/1‑251 (MQ=255)
                                                      ccGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACGGGCTAACAAGTCATACCCGAGAACGGTTTGCAGATTATGCT‑GCACCCCCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGAcg                                                                                                                                                                       <  2:103299/251‑1 (MQ=255)
                                                              gAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcg                                                                                                                                                               >  1:204241/1‑251 (MQ=255)
                                                               ggCATCAGGCGGCCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCCCCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcgg                                                                                                                                                              <  2:140226/250‑1 (MQ=255)
                                                                     aGGCGGGCAATCTGATCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcggcggaa                                                                                                                                                         >  2:270476/1‑248 (MQ=255)
                                                                                          aaTACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCACCCCCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTtg                                                                                                                                   <  1:188248/251‑1 (MQ=255)
                                                                                             aCTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCACCCCCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGc                                                                                                                                <  1:50966/251‑1 (MQ=255)
                                                                                                            aaaaagccaaacccGAGACCGGTTTGCAGATTATGCG‑GCCCCCCCCGGCTTGGCCCCCGCCACAGCGG‑TGTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGa                                                                                                                 <  2:237174/239‑1 (MQ=255)
                                                                                                                aGTCATAACCGAGAACGGTTGGCAGATTATGCT‑GCACCCCCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAAcgc                                                                                                             <  1:230758/251‑1 (MQ=255)
                                                                                                                           agaACGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCt                                                                                                                                                                                                               >  2:97221/1‑142 (MQ=255)
                                                                                                                           agaACGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCt                                                                                                                                                                                                               <  1:97221/142‑1 (MQ=255)
                                                                                                                            gcaggttttGCAGATCAGGCT‑GCCCCCCCCGGCTTGGCGCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGtg                                                                                                 <  2:10976/245‑1 (MQ=255)
                                                                                                                                          tgttggccgccccccccgccgcaGCTCCCCCCTCAGCGG‑TTTCAGGCGGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGATATGCGCCTGCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATTAGCTTcctc                                                                                     <  2:204241/228‑4 (MQ=255)
                                                                                                                                                  gccccccccGGCTTGGCGCCCGCCACAGCGT‑TTTCAGGCTGGTCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCCCCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTcc                                                                               <  1:218321/242‑1 (MQ=255)
                                                                                                                                                    accaccCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTTCGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGGGTGATGAGCTTCATCGTCTcccccgt                                                                          >  1:28758/1‑246 (MQ=255)
                                                                                                                                                                aGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCct                                                               <  1:42555/250‑1 (MQ=255)
                                                                                                                                                                   ccccGCCCCAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCctgctg                                                           <  2:101447/250‑1 (MQ=255)
                                                                                                                                                                                     cAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGGCTGATACGCCTGCTGCACGGTGCTGATACTg                                           <  1:270476/250‑1 (MQ=255)
                                                                                                                                                                                      aGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGAc                                         <  1:18224/251‑1 (MQ=255)
                                                                                                                                                                                                                             gg‑AACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAAttt  <  2:137320/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                              
ATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCT‑GCACCACCCGGCTTAGCTCCCGCCACAGCGG‑TTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGG‑AACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTT  >  NC_000913/4571033‑4571502

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