Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,405,417	(C)6→7	intergenic (‑345/+286)	nuoA ← / ← lrhA	NADH:quinone oxidoreductase subunit A/DNA‑binding transcriptional dual regulator LrhA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,405,411	1	.	C	89.2%	108.5 / 7.3	37	intergenic (‑339/+292)	nuoA/lrhA	NADH:quinone oxidoreductase subunit A/DNA‑binding transcriptional dual regulator LrhA
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base C (14/19);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CATTAGCATTTAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCT‑CCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTAACATCATTGTAGCAGGTGATTTTTTTCAGGCGAT  >  NC_000913/2405280‑2405532                                                                                                                                     |                                                                                                                          cATTAGCATTTAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATtcct‑c                                                                                                                          <  2:33576/133‑1 (MQ=255)  aTTAGCATTTAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATtcct‑cc                                                                                                                         >  1:387304/1‑133 (MQ=255)   ttAGCATTTAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCT‑ccc                                                                                                                        >  2:24930/1‑133 (MQ=255)   ttAGCATTTAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCT‑ccc                                                                                                                        >  1:48199/1‑133 (MQ=255)   ttAGCATTTAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCT‑ccc                                                                                                                        >  2:11362/1‑133 (MQ=255)           tAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCT‑CCCCCCcaaaa                                                                                                                >  2:494706/1‑132 (MQ=255)           tAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCT‑CCCCCCcaaaa                                                                                                                >  2:103321/1‑132 (MQ=255)           tAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCT‑CCCCCCcaaaa                                                                                                                >  1:75420/1‑132 (MQ=255)             aTGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAaga                                                                                                               <  1:334098/133‑1 (MQ=255)                  cTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTa                                                                                                          <  2:383658/133‑1 (MQ=255)                  cTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTa                                                                                                          <  1:21959/133‑1 (MQ=255)                  cTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTa                                                                                                          <  1:11362/133‑1 (MQ=255)                          ccAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTg                                                                                                  <  1:103321/133‑1 (MQ=255)                                   tGACGACTGCAGGATGATCAGAGGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGac                                                                                         <  1:220967/133‑1 (MQ=255)                                   tGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGac                                                                                         >  1:82324/1‑133 (MQ=255)                                   tGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGac                                                                                         <  1:494706/133‑1 (MQ=255)                                    gacgacTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATTCGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGaca                                                                                        <  2:48199/133‑1 (MQ=255)                                     acgacTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGacac                                                                                       <  2:501361/133‑1 (MQ=255)                                          tCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCttt                                                                                  >  1:506500/1‑133 (MQ=255)                                          tCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCttt                                                                                  >  1:427908/1‑133 (MQ=255)                                          tCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCttt                                                                                  >  1:263931/1‑133 (MQ=255)                                          tCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATGTTCTTGATCTGTGCCACGCttt                                                                                  >  1:201134/1‑133 (MQ=255)                                             aGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGt                                                                               <  2:75420/133‑1 (MQ=255)                                                  aaTTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTc                                                                          <  1:24930/133‑1 (MQ=255)                                                   aTTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTcc                                                                         <  2:263931/133‑1 (MQ=255)                                                    ttAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCa                                                                        >  2:380190/1‑133 (MQ=255)                                                         tgttgtTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGAGACGTTTATGTCATTCCATAACa                                                                   >  1:146540/1‑133 (MQ=255)                                                         tgttgtTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTGGATCTGTGACACGCTTGTGTCTTTCCATAACa                                                                   >  2:222205/1‑133 (MQ=255)                                                          gttgttGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACaa                                                                  >  2:41228/1‑133 (MQ=255)                                                            tgttgAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACaaaa                                                                >  1:313251/1‑133 (MQ=255)                                                                     gtgtCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGcaa                                                       <  1:418365/133‑1 (MQ=255)                                                                      tgtCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGcaac                                                      >  1:331359/1‑133 (MQ=255)                                                                                 aaGCAATGGAAAAAAAACGGGTCTATTTTGACGGGAATTTGGGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTa                                           <  2:201134/133‑1 (MQ=255)                                                                                       tGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACTAAAACGCGGCAGCAAATTTACTTAtttttt                                     >  1:179782/1‑133 (MQ=255)                                                                                         gAAAAAATCCGGGTCGATCTTGGCAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTaa                                   <  1:205137/133‑1 (MQ=255)                                                                                         gAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTaa                                   <  1:451853/133‑1 (MQ=255)                                                                                                               aCAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTAACATCATTGTAGCAGGTGAtttt             >  1:347654/1‑133 (MQ=255)                                                                                                               aCAGGAATTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTAACATCATTGTAGCAGGTGAtttt             >  1:54531/1‑133 (MQ=255)                                                                                                                  ggAATTTGTGTCGATTCCT‑CCCCCCAAATGAGCGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTAACATCGTTGTAGCAGGTGATTTTTTTc         <  2:174982/133‑1 (MQ=255)                                                                                                                     aTTTGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTAACATCATTGTAGCAGGTGATTTTTTTCAg       <  1:380190/133‑1 (MQ=255)                                                                                                                      tttGTGTCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTAACATCATTGTAGCAGGTGATTTTTTTCAgg      <  2:506500/133‑1 (MQ=255)                                                                                                                          tgtCGATTCCTCCCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTAACATCATTGTAGCAGGTGATTTTTTTCAGGCGat  <  2:387304/133‑1 (MQ=255)                                                                                                                                     |                                                                                                                          CATTAGCATTTAATGCCCTTTCACCCCAGATCCTTGACGACTCCAGGATAATTAGATGTTGTTGAATCGTGTCCGTTGTGAAGCAATGGAAAAAATACGGGTCTATTTTGACAGGAATTTGTGTCGATTCCT‑CCCCCCAAAAGAGAGTATTTTCTTGATCTGTGACACGCTTTTGTCATTCCATAACAAAAACGCAGCAACAAATTTACGTATTTTTTAACATCATTGTAGCAGGTGATTTTTTTCAGGCGAT  >  NC_000913/2405280‑2405532

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