Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,628,807	T→G	intergenic (+260/+48)	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_000913	4,628,807	0	T	G	93.8%	89.4 / ‑4.4	32	intergenic (+260/+48)	nadR/ettA	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (17/13);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGC‑TGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGC  >  NC_000913/4628564‑4629043                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             cAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAgg                                                                                                                                                                                                                                      >  1:204194/1‑251 (MQ=255)         aGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATa                                                                                                                                                                                                                               >  2:119861/1‑250 (MQ=255)                  aataagttttttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                                                                                       >  1:61836/7‑250 (MQ=255)                        ttttttGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                                                                                                >  1:25981/1‑251 (MQ=255)                                  aGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACACCGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtt                                                                                                                                                                                                      >  2:72598/1‑251 (MQ=255)                                   gTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTGTTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt                                                                                                                                                                                                     >  1:46514/1‑251 (MQ=255)                                   gTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAGGGCGTTCACGCCGCATCCGGCAttt                                                                                                                                                                                                     >  1:200834/1‑251 (MQ=255)                                                  ggaggagGCTGACGAGATAAACCAAGGTTTGAAAAAACCTCCACGGGTGCCCCGCTACGAAATTGCGGATACTAAGGCGTTTGTAGGCCGGAAAAGGCGTTCCCGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAGGGCGTTCACGCCGCATCCGGCATTTTACGCATTACtt                                                                                                                                                                                         <  2:200834/244‑1 (MQ=255)                                                        ggcgtcagagaATAAACAAGGTTTGAAAAACCCTGCACGGGGGACCCGCTACGAAATTGCGGAAACTAAGTCTGTGGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACg                                                                                                                                                                                <  2:121243/244‑1 (MQ=255)                                                          cTGTCGCGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGTCGCATCCGGCATTTTACGCATTACTTCGCAATACGCt                                                                                                                                                                              >  1:213514/1‑251 (MQ=255)                                                              cGAGATTAAACAAGGTTTAAAAAACCCTGCACGGGTGACCCGCTCCGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGt                                                                                                                                                                           <  1:72598/250‑1 (MQ=255)                                                                       aCAAGGTTTGAAAAAACCTACAGGGGTGACCCCCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACg                                                                                                                                                                 <  2:46514/251‑1 (MQ=255)                                                                                         tGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTc                                                                                                                                               >  2:73289/1‑251 (MQ=255)                                                                                          gCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCCGTGCGTCt                                                                                                                                              >  2:106752/1‑251 (MQ=255)                                                                                            aCGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAACGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATACGGCCTTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCCGTGCGTCTgc                                                                                                                                            >  2:95064/1‑251 (MQ=255)                                                                                                          ccaaaaTTGCGGATGTTAGGGCTGTGGAAGGCCGGAAAAGGCGTCCGCGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTTCCGCGGCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCGTGGTTT‑GTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCTTCTGCTCCCATCGTGCGtt                                                                                                                              <  2:194884/248‑1 (MQ=255)                                                                                                                      atacTAAGGTTGTTGTAGGCGGAAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCtt                                                                                                                  <  1:52066/248‑1 (MQ=255)                                                                                                                        acTAAGTCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTc                                                                                                                 <  2:236371/249‑1 (MQ=255)                                                                                                                           aaGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTAc                                                                                                             <  1:119861/251‑1 (MQ=255)                                                                                                                              gCTGTTGTAGGCGGGAAAAGCCGTTCCCGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCTGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTca                                                                                                          <  2:25981/251‑2 (MQ=255)                                                                                                                                     taggccGGAAAAGGCGTTCACGCCGTATCCGGAATTCAGTGCTCACTGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGATTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGACCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGt                                                                                                   >  2:17923/5‑251 (MQ=255)                                                                                                                                      aggccgaaaaaggcgtTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGAAGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATACACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGt                                                                                                   <  1:5586/236‑1 (MQ=255)                                                                                                                                               aaggcgtTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGCTTCACAACGTATTGAAATTTC‑TTGGTTTGTCGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCga                                                                                          >  1:163896/6‑250 (MQ=255)                                                                                                                                                         cgccgtATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTGTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACt                                                                               <  2:264393/251‑1 (MQ=255)                                                                                                                                                                  cGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTc                                                                       >  1:4745/1‑250 (MQ=255)                                                                                                                                                                             gCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTGCGAAGAACTCAACTTTACCTTCATCCTGGTagtc                                                           <  1:98574/251‑4 (MQ=255)                                                                                                                                                                                    gCCTGATGCGACGCTTACCGCGTATTATCAGGCCGGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAggat                                                     <  2:210024/250‑4 (MQ=255)                                                                                                                                                                                                gCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGAt                                         >  2:207907/1‑250 (MQ=255)                                                                                                                                                                                                         gTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTAcc                                                                          >  1:42728/1‑208 (MQ=255)                                                                                                                                                                                                         gTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGGTACCTTCGAAGAACTCAACTTTAcc                                                                          <  2:42728/208‑1 (MQ=255)                                                                                                                                                                                                                    cctggagattcaCAACGTATTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTc                    <  1:106752/240‑1 (MQ=255)                                                                                                                                                                                                                                      aTTGAAATTTC‑TTGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGc  >  1:153412/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             CAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGC‑TGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGC  >  NC_000913/4628564‑4629043

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