Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,707 C→A intergenic (+160/+148) nadR → / ← ettA DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,7070CA93.1% 74.0 / ‑0.1 29intergenic (+160/+148)nadR/ettADNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (15/12);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

TGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCCACGAAATTGAGGA‑TGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGG  >  NC_000913/4628472‑4628940
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         
tGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCAcgccgtat                                                                                                                                                                                                                           >  1:497523/1‑248 (MQ=255)
                           ccTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATg                                                                                                                                                                                                >  1:47133/1‑251 (MQ=255)
                             tgcgctgcgTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCg                                                                                                                                                                                              >  2:104251/1‑251 (MQ=255)
                                                               ggtgCAGATATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATTAGTTGAAATCCTCTCTGGAAGCGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGGGTCCCGCTACGAAATTGCGGA‑TACTAATGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtggagattca                                                                                                                                                            <  1:485711/248‑10 (MQ=255)
                                                                      agataccctcgataaaaCTGCTCAAAGGCGGGGTATAAAAGAAGTTTTTTAG‑TGAGTGGAAACCCTCTCTGGAAGAGGCTGCCGAGATTAAACAAGGTTGGAAAAAACC‑GGCACGGGTGACCCGCTACTAAATTGCGGA‑TACGAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTa                                                                                                                                                     <  2:269064/237‑1 (MQ=255)
                                                                      agaTAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtggagattcacaacgt                                                                                                                                                      >  1:94665/1‑235 (MQ=255)
                                                                        aTAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGAAAGAGCCTGTCGATATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTAtt                                                                                                                                                   <  1:294025/251‑1 (MQ=255)
                                                                                    aaaCTGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GAGGGGATGAAACCCTCTCTGGAAGAGGCTGTCGAAATTAAACAAGGTTTGAAAAACCC‑TGCAGGGGTGACCCGCTCCGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttgg                                                                                                                                       <  2:190454/251‑7 (MQ=255)
                                                                                         gCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGcc                                                                                                                                                                                                     <  2:427739/184‑1 (MQ=255)
                                                                                         gCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGcc                                                                                                                                                                                                     >  1:427739/1‑184 (MQ=255)
                                                                                                        ataAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCC‑TGCACGGGGGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtgga                                                                                                                                                                  <  2:325104/204‑4 (MQ=255)
                                                                                                        ataAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtgga                                                                                                                                                                  >  1:325104/1‑201 (MQ=255)
                                                                                                        ataAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCg                                                                                                                   >  1:387235/1‑251 (MQ=255)
                                                                                                        ataAAATAAGTTTTTT‑GAGGAGTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGc                                                                                                                    <  1:104251/250‑1 (MQ=255)
                                                                                                               aaagttttgtGATGAGTGTAAAACCTCTTGGGAAGAGGCTGTCGAGATAAAAAAAGGTTTGAAAAAACC‑TGCACGGGGGACCCGCTACGAAATTGCGGA‑TACGAGGGCGGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTcacg                                                                                                             <  2:120610/247‑4 (MQ=255)
                                                                                                                      tttGATGAGTTGAAAACCTCTCTGGATGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCTCGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTTGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                       <  1:462134/250‑1 (MQ=255)
                                                                                                                                                     cTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTcac                                                                                                              >  1:310207/1‑210 (MQ=255)
                                                                                                                                                        gcgaaaatAAAAAAGGTTTTAAAAAACCCTGCAGGGGTGCCCCGCCACGAAATTGAGGA‑TGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCtt                                                                      <  1:19146/243‑1 (MQ=255)
                                                                                                                                                           aaataaaaCAAGGTT‑GAAAAAACCCGGCAGGGGTGCCCCGCCACGAAATTAAGGA‑TGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACt                                                                 <  1:169341/246‑1 (MQ=255)
                                                                                                                                                             aTTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTg                                                               >  1:9455/1‑251 (MQ=255)
                                                                                                                                                             aTTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGCCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTCCCTCGCAATACGCTTGCACTTg                                                               >  1:471118/1‑251 (MQ=255)
                                                                                                                                                                 accAAGGTTTAAAAAACCC‑TGCCCGGGTGACCCGCAACGAAATTGCGGA‑TACTAGGGCTTTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATAc                                                           <  1:283516/249‑1 (MQ=255)
                                                                                                                                                                 aaCAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATAc                                                           >  1:82749/1‑251 (MQ=255)
                                                                                                                                                                           gAAAAAACC‑TGCACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTc                                                 >  2:441187/1‑251 (MQ=255)
                                                                                                                                                                              aaaaCC‑TGCAGGGGTGCCCCGCTAAGAATTTG‑GGACTACTAAGGCTGTGGTAGGCCGAAAAAGGCGTCCACGCCGTATCCGGCTTTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAg                                              <  2:268221/251‑1 (MQ=255)
                                                                                                                                                                                       cACGGGTGACCCGCTACGAAATTGCGGA‑TACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGACATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTg                                      >  1:189775/1‑251 (MQ=255)
                                                                                                                                                                                                                    tacTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTc          >  2:60468/3‑251 (MQ=255)
                                                                                                                                                                                                                            cTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATCCGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTc    >  2:116342/1‑249 (MQ=255)
                                                                                                                                                                                                                              gttgtAGGCGGGAAAAGGCGTTCCCGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCTTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCGTGTACTCTTCGTACTcag  <  1:127061/249‑3 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         
TGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACC‑TGCACGGGTGACCCGCCACGAAATTGAGGA‑TGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGG  >  NC_000913/4628472‑4628940

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