Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,085,828	G→A	8.7%	intergenic (‑6/‑188)	fdoG ← / → fdhD	formate dehydrogenase‑O, large subunit/formate dehydrogenase formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,085,828	0	G	A	8.7%	135.2 / 3.1	46	intergenic (‑6/‑188)	fdoG/fdhD	formate dehydrogenase‑O, large subunit/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base G (13/29);  new base A (2/2);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

TTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACA  >  NC_000913/4085680‑4085969                                                                                                                                                     |                                                                                                                                              ttACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTg                                                                                                                                               <  4:164974/149‑1 (MQ=255) ttACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTg                                                                                                                                               <  7:164065/149‑1 (MQ=255) ttACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTg                                                                                                                                               >  6:120836/1‑149 (MQ=255)       gTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGTGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCtt                                                                                                                                         <  3:187381/149‑1 (MQ=255)       gTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCtt                                                                                                                                         <  6:34616/149‑1 (MQ=255)              gggTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGa                                                                                                                                  <  7:247176/149‑1 (MQ=255)                   cgcAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACAtt                                                                                                                              <  4:97275/148‑1 (MQ=255)                   cgcAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTg                                                                                                                             <  6:55653/149‑1 (MQ=255)                    gcagcaGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGC‑AAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTc                                                                                                                           <  1:122760/149‑1 (MQ=255)                     cagcagTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAAGCCAGTGCCGCTGCCGTGGTGCGTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTc                                                                                                                           <  2:29393/149‑1 (MQ=255)                      agcagTTTATCCTGCCGGGGTTCCGCGAGTGCTACGCTGGGTGCAAACCCCAGTGCCGCTGCCGTGGGGCCTGCCATCCCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACATGCATGGATTGCTCCTTCTTTCGACATTGTCa                                                                                                                          <  6:182841/149‑1 (MQ=255)                            ttATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTg                                                                                                                                        >  5:161588/1‑129 (MQ=255)                            ttATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTg                                                                                                                                        <  6:161588/129‑1 (MQ=255)                                    ccGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACAt                                                                                                            <  6:215176/149‑1 (MQ=255)                                     cGGGTTTCAGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCAATTACAtt                                                                                                           >  1:235509/1‑149 (MQ=255)                                             cgcgAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTtgct                                                                                                   <  7:79528/149‑1 (MQ=255)                                               cgAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTtgctgc                                                                                                 <  2:67216/149‑1 (MQ=255)                                                    gCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGGGGTGCCTGCCATACAGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCgtgt                                                                                             <  5:33324/148‑1 (MQ=255)                                                     cTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGtgtg                                                                                            >  6:113082/1‑148 (MQ=255)                                                     cTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGtgtg                                                                                            >  8:116953/1‑148 (MQ=255)                                                      tACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGACCCATTTACATTCGCTTGCTGCGtgtg                                                                                            >  7:132263/1‑147 (MQ=255)                                                        cGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAgg                                                                                        <  4:112690/149‑1 (MQ=255)                                                                     aCCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTAt                                                                           <  8:197613/149‑1 (MQ=255)                                                                         aGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTAt                                                                           >  7:102582/1‑145 (MQ=255)                                                                             ccGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCtt                                                                    <  7:116953/148‑1 (MQ=255)                                                                               gCTGCCGAGGAGCCAGCCATACCGCCAGCGCAGATCGTACACAACTGCCTTCGGCTGACCTCCATGGATAGCTCCTTGATTCGACATTGTCACCTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGc                                                                 <  2:235509/149‑1 (MQ=255)                                                                                   ccGTTGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGat                                                             <  3:169485/149‑2 (MQ=255)                                                                                      tggtgCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAggg                                                          <  5:120836/149‑1 (MQ=255)                                                                                         tgcctgccATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGAGTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGcc                                                       <  5:113082/149‑1 (MQ=255)                                                                                              gccATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTACTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCa                                                    >  1:60047/1‑147 (MQ=255)                                                                                              gccATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTACTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCa                                                    >  1:122806/1‑147 (MQ=255)                                                                                                 aTACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACCTTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGt                                               >  2:30566/1‑149 (MQ=255)                                                                                                 aTACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTACTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGt                                               <  2:60047/149‑1 (MQ=255)                                                                                                 aTACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTACTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGt                                               <  2:122806/149‑1 (MQ=255)                                                                                                  tACGGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTc                                              <  3:95962/149‑1 (MQ=255)                                                                                                  tACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTGCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTc                                              <  3:92508/149‑1 (MQ=255)                                                                                                            gcAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGc                                    >  2:108174/1‑149 (MQ=255)                                                                                                                     aaaGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAgg                                                                                        >  2:98986/1‑88 (MQ=255)                                                                                                                     aaaGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAgg                                                                                        <  1:98986/88‑1 (MQ=255)                                                                                                                     aaaGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTcc                           >  8:71214/1‑149 (MQ=255)                                                                                                                       aGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCt                         >  4:198153/1‑149 (MQ=255)                                                                                                                              ccTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCt                                                                     <  4:196755/98‑1 (MQ=255)                                                                                                                              ccTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCt                                                                     >  3:196755/1‑98 (MQ=255)                                                                                                                              ccTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCTGCAAAGGTCCAGAATTGGCTCAATTTCCCtccatc                    >  7:68724/1‑147 (MQ=255)                                                                                                                               cTTCTGCTGAACTGCATGGATTGCGCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAg                 <  5:133030/149‑1 (MQ=255)                                                                                                                                        aCCTGCATGGATTGCTCCTTGTTTCTACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGCCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGtt         <  5:204014/148‑1 (MQ=255)                                                                                                                                          cTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGGGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGt      <  2:235585/149‑1 (MQ=255)                                                                                                                                             cATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTaa    <  3:198153/148‑1 (MQ=255)                                                                                                                                              aTGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACa  <  4:230073/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              TTACGGGTTTCGCGGGTGCGCAGCAGTTTATACTGCCGGGTTTCCGCGAGTGCTACGCTGGGTGCAAAACCCAGTGCCGCTGCCGTGGTGCCTGCCATACCGCCAGCGCAGATCTTAAAGAACTGCCTTCTGCTGACCTGCATGGATTGCTCCTTGTTTCGACATTGTCACGTCCCATTTACATTCGCTTGCTGCGTGTGCAGGGAGTGGGAGTTATTTTTCTTTGCGGAAGGGGCCGCAAAGGTCCAGAATTGGCTCAATTTCCCTCCATCCAGGAAGGGTTTGTAACA  >  NC_000913/4085680‑4085969

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