Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,805 G→T intergenic (+258/+50) 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,8050GT94.4% 48.5 / ‑5.2 18intergenic (+258/+50)nadR/ettADNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (12/5);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAG‑TTTTACAACGTATTGAAATTGCT‑GGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACC  >  NC_000913/4628564‑4629038
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         
cAAAGGAGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAg                                                                                                                                                                                                                                  <  1:310660/250‑1 (MQ=255)
         ggTATAAAATAAGTTTTTT‑GAGGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTAAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAg                                                                                                                                                                                                                        <  1:210425/251‑1 (MQ=255)
                 aataagttttttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGCGACCCGCTACGAAATTGCGGCTACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGACCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATACGGCGTTc                                                                                                                                                                                                                  >  1:156731/7‑250 (MQ=255)
                                          ccTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCAtt                                                                                                                                                                                        >  1:351445/1‑251 (MQ=255)
                                                   aaGAGGCGGTCGAAATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTAAGAAATGGCGGAAACTAAGGTTTTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCaa                                                                                                                                                                               <  1:346790/251‑1 (MQ=255)
                                                                                 aaaacCCTGCCCGGGTGACCCGCTACGAAATTGCGTATACAAAGGCGGTTGTAGGCCGGAAAAGGCTTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCa                                                                                                                                                 <  1:248167/251‑1 (MQ=255)
                                                                                   aaaCCTGCCCGAGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGACTGATGCGTCGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACACCGTATTGCAATTTCTTGGATTT‑TAGGCCGGATTAGGCGTTCACGCCGACTCCGGCATTTTACGCATTACTTCGCAATACGCTTGTGCTTTCTACGCTTTGGcacaagt                                                                                                                                               >  2:306818/1‑245 (MQ=255)
                                                                                      ccTGCACGGGTGACCCGCTACGAAATTGCGGATACTATGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCAATCGGTGCTCAATGCCTGATGCGACGATTACCGCGTCTTATCAGACCTGGAG‑ATTCACAACGGATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCATTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTcccgtacg                                                                                                                                            >  1:8296/1‑245 (MQ=255)
                                                                                                                   ggATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGGGCGTTTCTTGTActc                                                                                                               >  1:58932/1‑251 (MQ=255)
                                                                                                                        cTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGt                                                                                                          >  1:207791/1‑251 (MQ=255)
                                                                                                                                                   ggtCACGCCGTAACCGCCATTCGATGCTCAATGCCTGATGCGACGCTTACCGCGTATTATCAGGCCT‑CAGGATTCACAACGTATTGGAATTTCT‑TGGTTTGTAGGCCGGCTAAACCCTTCACGCCGTAACCCGCATTTTACGCATTGCTTAGCAGTACGCTTGTCCTTCCTACGCCTCGGCTCCAGTGCGACGGCGACCAGAGTGGGTTTCTTGCACTATTCGTACTCgtacagatcagatacttagcact                                                                               >  2:128954/3‑229 (MQ=255)
                                                                                                                                                      cACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCaa                                                                            <  1:58808/251‑1 (MQ=255)
                                                                                                                                                                 cGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTc                                                                  >  1:44447/1‑250 (MQ=255)
                                                                                                                                                                     aTTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCt                                                             >  1:116194/1‑251 (MQ=255)
                                                                                                                                                                                               gCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGAt                                    >  1:201956/1‑250 (MQ=255)
                                                                                                                                                                                               gCTTACCGCGTCTTATCAGGCCTGGAG‑ATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGAt                                    >  1:186645/1‑250 (MQ=255)
                                                                                                                                                                                                                          gattcaCAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTACAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGTAACCAACGGTCGTGCgg         >  2:361071/6‑250 (MQ=255)
                                                                                                                                                                                                                                 aaCGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAAcc  >  2:86123/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                         
CAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAG‑TTTTACAACGTATTGAAATTGCT‑GGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACC  >  NC_000913/4628564‑4629038

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