Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,834,329	G→A	D140N (GAT→AAT)	norW →	NADH:flavorubredoxin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,834,329	0	G	A	93.5%	88.8 / NA	31	D140N (GAT→AAT)	norW	NADH:flavorubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/18);  minor base C (2/0);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCGCCAGACGGCGGGTGAATTTGCCGAGCAGTTTAATCTGCACCTGTTTCCACAAACCTGGGTGACGGATATCGATGCCGAAGCCCGTGTGGTGAAAAGCCAGAATAATCAGTGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGGATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAATCTCAGTTACAGGGGCTGGAAAAAACGGATTCTGGCATTCAGGCAACGCTGGAC  >  NC_000913/2834089‑2834574                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      aCCCGCCAGACGGCGGGTGAATTTGCCGAGCAGTTTAATCTGCACCTGTTTCCACAAACCTGGGTGACGGATATCGATGCCGAAGCCCGTGTGGTGAAAAGCCAGAATAATCAGTGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGAc                                                                                                                                                                                                                                              <  1:72401/250‑1 (MQ=255)                     aTTTGCCGAGCAGTTTAATATGCACCTGTTTCCACAAACCTGGGTGACGGATATCGATGCCGAAGCCCGTGTGGTGAAAAGCCAGAATAATCAGTGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGGGTTGATTGTTGGCggtg                                                                                                                                                                                                                         >  2:30873/1‑251 (MQ=255)                            gAGCAGTTTAATATGCCCCTGTTTCCACAAACCTGGGTGATGGATATCGATGCCGACGCCCGTGTGGTGAAAAGCCAGAATAATCAGTGGCAATACGACAAGCTAGTCATAGCAAACGGTGCCAGGGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGCCGTTCAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCACCTGCGGCATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGAt                                                                                                                                                                                                                  >  2:101543/1‑251 (MQ=255)                                 gTTTAATCTGCACCTGGTTCCACAAACCGGGGTGACGGATAGCGATGCCGAAGCCCGTGTGGTGAAAAGCCAGAATAAGCAGTGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCGTTGGCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGACTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTa                                                                                                                                                                                                             <  1:113969/251‑1 (MQ=255)                                     aaTCTGCACCTGTTTCCAAAAACATGGGTGACGGATAGCGATGCCGAAGCCCGTGTGGTGAAAAGCCAGAATAATCAGGGGCAATACGACAAGCTAGTACTGGCAACCGGTTCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCg                                                                                                                                                                                                          <  2:208853/250‑1 (MQ=255)                                      aTCTGCACCTGTTTCCACAAACCTGGGTGACGGATATCGATGCCGAAGCCCGTGTGGTGAAAAGCCAGAATAATCAGTGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCACGAGTATCGCGCCTGTGAAACGAAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGCGGTTTGATTGGTATCGa                                                                                                                                                                                                         >  2:77501/1‑250 (MQ=255)                                      aTCTGCACCTGTTTCCACAAACCTGGGGGACGGATATCGATGCCGAAGCCCGTGTGGTAAAAAGCCAGAATAATCAGAGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCTTTGTACCGCCTGTGCCTGGGCGTGGGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGAGGAATGACCGACCCGTGTTGATTGTTGGCGGAGGTTTGTTTGGTAGCGa                                                                                                                                                                                                         >  2:23347/1‑250 (MQ=255)                                                gttgcaccaAACCTGGGTGACGGATATCGATGCCGAAGCCCGTGGGGTGAAAAGCCAGAATAACCAGGGGCAATACGACAAGCTAGTACTGGCAACCGGTGCAAGTGCCTTGGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATgg                                                                                                                                                                                              <  2:194091/244‑1 (MQ=255)                                                   cccACAAACCAGGGGTACAGATAGCGATGCGGAAGCCCGCGTGGTGAAAAGCCAGAATATTTAGTGGCAATACGACAAGCTAGTACTGCCAACCGGGGCCCGTGCCTGTGGCCCGCATGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTTCGGAATGCCCGACGCGTTTTGATTGGTGGCGGTGGTTGGATTGGTCGCTAACTGGCGATGGAtt                                                                                                                                                                                           <  2:100221/250‑1 (MQ=255)                                                                               gCCGAAGCCCGTGTGGTAAAAAGCCAGAAAAATCAGTGGCACTACGACAAGCTAGTACGGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTTACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCt                                                                                                                                                               <  2:205170/251‑1 (MQ=255)                                                                                          tgtgGTGAAAAGCCAGAATAATCAGTGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACg                                                                                                                                                    >  1:157584/1‑251 (MQ=255)                                                                                            gggtgAAAAGCCAGAATAATCAGTGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGGTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACGGCGGAATGCCCGTCGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCt                                                                                                                                                  <  2:122833/250‑1 (MQ=255)                                                                                                       cAGAATAATCAGTGGCAATACGACAAGCTAGAACTGGCAACCGGAGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAACTAGTCAGCAAGAGCATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGGTTGGTCGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAATCGGTCACGCTAATCGGCAACGCTGCCAATATTcg                                                                                                                                       >  2:214567/1‑250 (MQ=255)                                                                                                          aataatCAGTGGCAATACGACAAGCTAGTACTGGCACCCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGc                                                                                                                                    <  2:38968/251‑1 (MQ=255)                                                                                                                     gCAATACGACAAGCTAGTACTGGCACCCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGATCTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGc                                                                                                                         <  2:238132/251‑1 (MQ=255)                                                                                                                                   tAGTACTGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGCACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAagc                                                                                                           >  1:97538/1‑251 (MQ=255)                                                                                                                                         tGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGTCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGTTTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTACTCGACAACGCTGCCAGTCTTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGc                                                                                                     >  2:66267/1‑251 (MQ=255)                                                                                                                                                  gTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCAt                                                                                            >  1:71090/1‑251 (MQ=255)                                                                                                                                                  gTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGAGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCAt                                                                                            >  1:88951/1‑251 (MQ=255)                                                                                                                                                                         gtgCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCTGCTTGCAGCATCGGTTGACGGAGATGGGCGtt                                                                       <  2:189559/249‑1 (MQ=255)                                                                                                                                                                           gCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCAt                                                                    <  2:71090/250‑1 (MQ=255)                                                                                                                                                                                 gcgcgAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGt                                                             <  2:245219/247‑1 (MQ=255)                                                                                                                                                                                         tAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAAtct                                                     >  1:8672/1‑251 (MQ=255)                                                                                                                                                                                             gCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGTTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAATCTCAGt                                                 <  1:178005/251‑1 (MQ=255)                                                                                                                                                                                                   gTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGCATGCCAGACGCGTGTTGATTGTTAGCGGAGGTTTGATTGGTAGAGAACTGGCGATGGATTTTTGTAGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAATCTCAGTTACAgg                                           >  1:166060/1‑251 (MQ=255)                                                                                                                                                                                                          aGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAATCTCAGTTACAGGGGCTGGa                                    <  1:117019/251‑1 (MQ=255)                                                                                                                                                                                                           gTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCTTGTTGATTGGTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGGATTCTTGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAATCTCAGTTACAGGGGCTGGaa                                   <  1:157655/251‑1 (MQ=255)                                                                                                                                                                                                                        tCGCGCCTGTGAAACGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAATCTCAGTTACAGGGGCTGGAAAAAACGGATTCTg                      >  2:39224/1‑251 (MQ=255)                                                                                                                                                                                                                              atgtgAAACGCAACTGCGGAATGCCCGACGCGTGTTGCTTGTTGGCGGAGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAATCTCACTTACATGGGCTGGAAAAAACGGATTCTGGCATTc                <  1:69794/250‑1 (MQ=255)                                                                                                                                                                                                                                      cGCAACTGCGGAATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCAGCTGCTGTTGAAATCTCAGTTACAGGGGCTGGAAAAAACGGATTCTGGCATTCAGGCAACg        <  2:97538/251‑1 (MQ=255)                                                                                                                                                                                                                                             gCGGAATGCCCGAGGCGTGTTGATTGTTGGCGGTGGGTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCAGCTGCTGTTGAAATCTCAGTTACAGGGGCTGGAAAAAACGGATTCTGTCATTCAGGCAACGCTGGAc  <  2:234009/250‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      ACCCGCCAGACGGCGGGTGAATTTGCCGAGCAGTTTAATCTGCACCTGTTTCCACAAACCTGGGTGACGGATATCGATGCCGAAGCCCGTGTGGTGAAAAGCCAGAATAATCAGTGGCAATACGACAAGCTAGTACTGGCAACCGGTGCCAGTGCCTTTGTCCCGCCTGTGCCTGGGCGTGAGTTAATGCTGACGTTAAATAGTCAGCAAGAGTATCGCGCCTGTGAAACGCAACTGCGGGATGCCCGACGCGTGTTGATTGTTGGCGGTGGTTTGATTGGTAGCGAACTGGCGATGGATTTTTGTCGTGCAGGCAAAGCGGTCACGCTAATCGACAACGCTGCCAGTATTCTGGCGTCGTTAATGCCACCGGAAGTAAGCAGCCGCTTGCAGCATCGGTTGACGGAGATGGGCGTTCATCTGCTGTTGAAATCTCAGTTACAGGGGCTGGAAAAAACGGATTCTGGCATTCAGGCAACGCTGGAC  >  NC_000913/2834089‑2834574

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