Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,514 A→G E400E (GAA→GAG 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,5140AG96.0% 74.7 / NA 25E400E (GAA→GAGnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (13/11);  minor base T (0/1);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCT  >  NC_000913/4628290‑4628738
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                
gagcATCCGTTCGTACAGGCGTTGATTGCTGAATACCGTTGCGAGCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGGTGAAGAGGACGATTATGACAGCCGTGTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCa                                                                                                                                                                                                        <  2:138130/248‑1 (MQ=255)
                   cgttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACACAACCCCGTGGGTGGCGGATGGTTTAAGCAGCCTAGGCCGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATAGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGCTGGGGGAGCAGAGATAACCGCGCTGaaa                                                                                                                                                                                      >  1:121419/4‑250 (MQ=255)
                   cgttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAAc                                                                                                                                                                                     <  1:192368/248‑1 (MQ=255)
                    gttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTAAGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGAGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAAcg                                                                                                                                                                                    >  1:120180/3‑250 (MQ=255)
                     ttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGAGGATCGCCAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAAcgg                                                                                                                                                                                   >  1:176433/2‑249 (MQ=255)
                              gaATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGc                                                                                                                                                                          <  1:70569/251‑1 (MQ=255)
                                     gTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGtat                                                                                                                                                                   >  1:264393/1‑251 (MQ=255)
                                     gTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGtat                                                                                                                                                                   <  1:91059/251‑1 (MQ=255)
                                               ggggATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGGGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATGATGACAGCCGTTTCCTGCGCTGCTTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAAta                                                                                                                                                             <  1:252470/244‑2 (MQ=255)
                                                           ggAGAACAACACGCCGTGGGGGGCGGATTGTTTACGGAGCCTCGGCAGGTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCTTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGATATAACCGCGATGAAACGGCTCAAAGGCGAGTGATAAAATAAGTTTTTTGATGAGtt                                                                                                                                            <  2:249719/251‑1 (MQ=255)
                                                                     acGCCGAGGGTGGAGGATGGTTAACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGATGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGAGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCGGAGATAACCGCGATGAAACGGCCCAAAGGCGAGGAATAAAATAAGTTTTTTGATGAGTTGAAAACctct                                                                                                                                  >  2:191165/1‑251 (MQ=255)
                                                                          gtgggtggCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGCGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTgg                                                                                                                              >  2:21223/1‑250 (MQ=255)
                                                                                  cGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTg                                                                                                                     >  1:101854/1‑251 (MQ=255)
                                                                                                          gTTCGGTTGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGATCTGTTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGATTTGAAAACCTCTCTGGAAGAGGCTGTCGAGAGTAATGAAGGGTTGaaaa                                                                                             <  2:159911/251‑1 (MQ=255)
                                                                                                              ggTGGATCGCAAAGAGTTCCAGAACTTGCTTGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAATAAAcc                                                                                         <  2:120180/251‑1 (MQ=255)
                                                                                                              ggTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAAcc                                                                                         <  1:88501/251‑1 (MQ=255)
                                                                                                                                                             agaACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCtgttgt                                            >  1:74131/1‑249 (MQ=255)
                                                                                                                                                                 cAATATCGGATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGCGCTGGTGCGGGAGATGCTGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAACGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGCGGCTGTCGGGATTAAACAAGGTTTAAAAAACCCTGCACGGGTGACCCGCTACGAACTTGCGGATACTAAGGCTGTTGTAGTCCg                                      >  2:32893/1‑251 (MQ=255)
                                                                                                                                                                        gaattcgggCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGaaaagg                               <  1:241303/243‑4 (MQ=255)
                                                                                                                                                                                 cGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCAcgc                      >  1:30078/1‑251 (MQ=255)
                                                                                                                                                                                      tGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCGCGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCAcgccgta                  <  1:21223/250‑3 (MQ=255)
                                                                                                                                                                                                 aTTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCgg      >  1:67974/1‑251 (MQ=255)
                                                                                                                                                                                                 aTTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACATCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGCGATCCGCTACGAAATTGCGGCTACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCgg      >  2:76106/1‑251 (MQ=255)
                                                                                                                                                                                                  ttatgACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGt     >  1:236371/5‑251 (MQ=255)
                                                                                                                                                                                                      gACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGGGATGAAACGGCTTAAAGGCGAGGTAAAAAATAAGTTTTTGGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCTTTCACGCCGTATCCGGCATTCGGTGCt  <  2:101854/250‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                
GAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCT  >  NC_000913/4628290‑4628738

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