Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,487 C→T Y391Y (TAC→TAT 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,4870CT90.9% 50.8 / ‑1.4 22Y391Y (TAC→TATnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (14/6);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

TCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAC  >  NC_000913/4628259‑4628716
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     
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)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                     
TCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAC  >  NC_000913/4628259‑4628716

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