Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,481 A→C E389D (GAA→GAC nadR → DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,4810AC91.3% 47.3 / ‑1.3 23E389D (GAA→GACnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (15/6);  total (16/7)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

GATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAC  >  NC_000913/4628239‑4628716
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
gacactgATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGACGAGGACGATTatg                                                                                                                                                                                                                                      >  1:87440/7‑248 (MQ=255)
                    ccAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCTGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTTCGACAGTCGTTTCctgcgctgcg                                                                                                                                                                                                                 <  1:149606/250‑1 (MQ=255)
                                                    agcATCCGTTCGTACAGGCGTTGATTGCTGCATAACGTTTCGATCTGGTGATCATGCTGAAGAACAACACGCCGTGGGTGGCGGATGGTTTACGAAGCCTCGGAAGTTCGGTGGGTCGCAAAGAGTTCCAGAACTTGATGGTGGAGATGCTCGAAGAGAACCATATCGACTTCGAGCGGGTTGAAGAGGCCGATTATAACAGCCGTTTCCTGCGCAGCGTGGAGCTGGTGCGGGGGATGATGGGGGAGcc                                                                                                                                                                                  >  1:244292/3‑249 (MQ=255)
                                                           gTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAAc                                                                                                                                                                          >  1:349963/1‑251 (MQ=255)
                                                                  caggggttgagtgatgaATCCGGTGTCGATCTGGTGGTCGTGCTGGAGAACAACACGCCTTGGGTGGCGGATGGGTTACGCAGCCTCGGCCGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATTCTGGAAGAGAACAATATCGAATTCGGGCGGGTTGAAGAGGACGATTATGACAGCCGTTGCCTTCGCTGCGGGGAGCTGTTGCGGGAGATGATGGGGGAGCATAGATAACCGcgatg                                                                                                                                                                   <  2:183060/240‑1 (MQ=255)
                                                                  cAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCTAAGAGTTCCAGACCTTGCTGGTGGAGATGCTGGAAGAGCACAATATCGACTTCGTGCGGGTTGCAGAGGACGATTATGACAGCCGTTTCCTGCGATGCGTGGCGCGGGTGCGGGGGTTGATGGGGGacctgagtcaaccgcgatg                                                                                                                                                                   >  2:210425/1‑233 (MQ=255)
                                                                   aGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGCGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCGAGAACTTGCTGGTGGCGATTCTGGAAGTGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGAAAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGGGGGAGATGGTGGGGGAGCAGAGATAAcagagatg                                                                                                                                                                   >  2:310660/1‑243 (MQ=255)
                                                                       gttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAAcg                                                                                                                                                              <  1:181037/249‑2 (MQ=255)
                                                                                                cTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATCACCGCGCTGAAACGGCTCAAAGGCGAGGAATAAAAta                                                                                                                                       >  1:351594/1‑248 (MQ=255)
                                                                                                      aTCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGtttttt                                                                                                                               >  1:359247/1‑251 (MQ=255)
                                                                                                                    ccacacGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGAGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGCTTATGAAAGCGGTTTCCTGCGCTGCGTGGGGCTGGTGCGGGAGATGATGGCGGAGCAGAGATAACCGCGATGACACGGCTCAAAGGCGTGGGAAAAAATACGTTTTTTGATGAGTTGaaacc                                                                                                                >  2:86074/3‑249 (MQ=255)
                                                                                                                     aacaCGCCGTGGGTGGCGGATGGATTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGGTTTTTGATGAGTTGAAAAcc                                                                                                               >  1:154223/1‑251 (MQ=255)
                                                                                                                           ccGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAACCCTCTCTg                                                                                                         >  1:362426/1‑251 (MQ=255)
                                                                                                                                          ggTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATaaa                                                                                                                                          <  1:119127/204‑1 (MQ=255)
                                                                                                                                          ggTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGAATaaa                                                                                                                                          >  2:119127/1‑204 (MQ=255)
                                                                                                                                                   aGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACa                                                                                  >  1:4150/1‑250 (MQ=255)
                                                                                                                                                                    ggATCGCAAAGAGTTCCAGCACATGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGCGGACGATTATGACAGCCGTTTCATGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGTGATAACCGCGTTGAAACGGCTCAAAGGCGAGGAATAAAATAAGTTTTTTGATGAGTTGAAAACCTCACTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCtcc                                                                >  2:248167/1‑249 (MQ=255)
                                                                                                                                                                        cGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATGTCGAATTCGTGCGGGTTGAAGAGGACGATGATGACAGCCGTTTACTGCGCTGCGTTGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCCCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACgg                                                            <  1:24144/251‑1 (MQ=255)
                                                                                                                                                                                    cAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGctac                                                <  1:123736/251‑4 (MQ=255)
                                                                                                                                                                                                  ggAGATGCTCGAAGAGATCAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGCACTGGTGCGGGAGATGATGTGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATg                                   >  1:142978/1‑251 (MQ=255)
                                                                                                                                                                                                                                cgtgCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGGTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTc    <  1:336477/250‑1 (MQ=255)
                                                                                                                                                                                                                                  tgCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCAc  >  1:115094/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tgCGGGTTGAAGAGGACGATTATGACAGCCGATTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGATGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGGGCCCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGTCCGGAAAAGGCGTTCAc  >  1:292751/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
GATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAC  >  NC_000913/4628239‑4628716

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: