Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,377 A→G V21V (GTA→GTG)  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,627,3770AG96.5% 87.2 / NA 29V21V (GTA→GTG) nadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/17);  minor base T (1/0);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGAC  >  NC_000913/4627133‑4627617
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
ccGGTGCCCAGCGTTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTAGGCGTTAAAAAAATCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATAAGGAGTGAGATAGGTCGTCATTTGATTACCTGAACACTGCCATCAAGCAAATGGACTGCAAGCTACAGCAGGTGGCTGAt                                                                                                                                                                                                                                            >  2:207595/1‑251 (MQ=255)
    ggCCCGGCGGTCAGGAACGAATTTCAGCAGCGGCGAAAAACGGTGTTCGCCGGGGGATGGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTGTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAAGACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCa                                                                                                                                                                                                                                        <  2:188408/250‑1 (MQ=255)
      cccAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGc                                                                                                                                                                                                                                      >  1:350511/1‑251 (MQ=255)
         aGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGt                                                                                                                                                                                                                                   <  1:371700/251‑1 (MQ=255)
                        atttcagAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGTGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTAt                                                                                                                                                                                                                    <  1:207164/245‑1 (MQ=255)
                                 gcggcgAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAg                                                                                                                                                                                                           <  1:242210/251‑1 (MQ=255)
                                      gAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACt                                                                                                                                                                                                      <  1:246524/251‑1 (MQ=255)
                                       aaaCACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTg                                                                                                                                                                                                     >  1:252518/1‑251 (MQ=255)
                                                                       ggcTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACGCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGTGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCaaaa                                                                                                                                                                     <  1:193848/251‑1 (MQ=255)
                                                                                       aaaGCGCCGGAAGAGATGCAGATTTTTGACGTTAAAAAACTCTACGATGAGCTTAGCGTGTTCGAAGACTTATAATGAGGGATACGGAGGGAGATATGTTGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCAGAACGCTACAGCAGGAGGCTGATGCCAGCGGGATGACCAAAGGGTATTAAAGCCAGTGACTGAATGCAAAAATACAATGCCCCAGGGCGCAAAAGCTGGAGGCGTTgaac                                                                                                                                                     >  2:193118/1‑248 (MQ=255)
                                                                                           cgcCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGGGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTTGCTGATGCCAGCGGTATGACCACAGGGTGTTGACGCAAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGAGCAAAAGCTGGAGGCGTTGCACCGtg                                                                                                                                                 >  2:149187/1‑250 (MQ=255)
                                                                                               ggAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGtttttt                                                                                                                                             >  1:100762/1‑251 (MQ=255)
                                                                                                gAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTg                                                                                                                                            <  1:17380/251‑1 (MQ=255)
                                                                                                   gggATGCAGATTTTTGGCGTTAAAAAACTCGCCGACGCGCTAAGAGTGTTCGACGACTTATAATGGGGAATACGGAGGGAGATATTTCGTCATTTGATTACCTGAAAACTGCCATCACGCAACAGGGATGCACGCTAAAGCAGGTGGCTGATGCCAGCGGTATGCCCAAAGGGTAGTTACGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGATGCACCGGTTTTTgtgg                                                                                                                                         >  2:311302/1‑248 (MQ=255)
                                                                                                    tgATGAAGATTTTTGGCGTTAAACAACTCTCAGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGATGGAGATATGTCGTCATTTGATTACCTGAAAGCTGCCATCAAGAAACAGGGCTGCATGCTACAGCAGTTGGCTGATGCCAGCGGTATGCCCACAGGGTATTTAAGCCAGTTACTGTATGCCTAAATCAAAAGCCCCATCGCGAACAAGCTGGAGGCGTGGCACCGTTTTTAGGGGc                                                                                                                                        >  2:314910/2‑251 (MQ=255)
                                                                                                           gATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGTGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACTGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTGTGGGGCTTGAGtt                                                                                                                                 <  1:67945/251‑1 (MQ=255)
                                                                                                           gATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGtt                                                                                                                                 <  1:116860/251‑1 (MQ=255)
                                                                                                                  ggCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTAt                                                                                                                                                                                                                                 <  1:351541/148‑1 (MQ=255)
                                                                                                                  ggCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTAt                                                                                                                                                                                                                                 >  2:351541/1‑148 (MQ=255)
                                                                                                                                                     cgacgtCTTATAATGAGGAATAAGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCATCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGt                                                                                        <  2:149417/250‑1 (MQ=255)
                                                                                                                                                     cgacgaCTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGtt                                                                                       <  1:380341/251‑1 (MQ=255)
                                                                                                                                                     cgacgaCTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGtt                                                                                       <  1:51748/251‑1 (MQ=255)
                                                                                                                                                                                      ttttcgtcATTGGATTACTCGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGGCCAAAGGGTATTTAAGCCAGTTATCGAATGCCAAAATCAAAAGCCTCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTa                                                       <  2:379366/247‑1 (MQ=255)
                                                                                                                                                                                       tGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGATTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGAAAGTTCTACCCGCTGCATACCGGACATATCTCCCTTATc                                                     >  2:233267/1‑251 (MQ=255)
                                                                                                                                                                                                  gATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTg                                          <  1:233267/251‑1 (MQ=255)
                                                                                                                                                                                                  gATTACCTGAAAACTGCCATCAAGCAACAGGGATGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTg                                          <  2:100762/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGAATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAACGCTGGAGGCGGTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGGGTGGTATTCGTTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCAAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATgg          >  2:208543/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGtt       <  1:321908/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGAc  >  1:166284/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
CCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGAC  >  NC_000913/4627133‑4627617

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: ‑