Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,550,663	0	A	C	25.5%	85.0 / 13.9	47	N68T (AAC→ACC)	fdnH	formate dehydrogenase N subunit beta
Reads supporting (aligned to +/- strand):  ref base A (29/6);  new base C (0/12);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 3.55e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCT  >  NC_000913/1550526‑1550802                                                                                                                                          |                                                                                                                                            ggTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGGCGTTCCCCCCTGTATGGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACACCGGTGATGAGGTGGGCCACTGCGTCGGGGTTTACGATAccc                                                                                                                                           <  2:332824/140‑1 (MQ=25)   tGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAAcccc                                                                                                                                         >  1:672158/1‑140 (MQ=255)       tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                                     >  1:1194617/1‑140 (MQ=255)       tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                                     >  1:332824/1‑140 (MQ=255)       tGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGCGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCccgccg                                                                                                                                     >  1:282262/1‑140 (MQ=255)                aGCAGAAGTCGCAAAAGTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAgcg                                                                                                                            >  2:658664/1‑140 (MQ=255)                 gCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAgcgc                                                                                                                           >  1:752346/1‑140 (MQ=255)                  cAGAAGTCAAAAAGCTTACCGCCTTTTCCACTTGTACGGGCTTTAAAGCCTGTCGGGTGGCGGGTGCGGAGGGAAACCACCCCCGTGATAAGGGGGGCCACTGCGTGGGGGTTTCCAATACCCCCGCCGATCTGAGCGcc                                                                                                                          <  2:1283238/140‑1 (MQ=255)                  aaGAAGTCGCAAACTTTATCGACGTTTCCCCCTGTATCGGCTGTAAAGCCTGTCGGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAGGGGGGGCACTGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGcc                                                                                                                          <  2:445272/139‑1 (MQ=25)                  aaGAAGTCGCAAAACTTATCGACGTTTCCCCCTGTATGGGCTGTAAAGCTTGTCGGGTGGCGTGTTCGGAGTGGAACGACATCCGTGAGGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGcc                                                                                                                          <  2:1194617/139‑1 (MQ=255)                            aaaaCTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTGGGGTGGGGTGTTCGGAGTGGAACGACATCCGTGATGAGGTGGGGAACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGa                                                                                                                <  2:1195452/140‑1 (MQ=25)                               aCTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCGGGGGGGGTGTTGGGAGTGAAACGACTTCGGTGATGAGGGGGGGCCCGGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACgg                                                                                                             <  2:278636/140‑1 (MQ=25)                                   aTCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGAt                                                                                                         >  1:857132/1‑140 (MQ=255)                                              ccccTGTCTCGGCTGTAAGGCTTGTCAGGGGGGGTGTTCGGAGTGGAACGACATCGGTGATGAAGGGGGGCACGGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCg                                                                                              <  2:711077/138‑1 (MQ=255)                                                   gTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAAcc                                                                                         >  1:670280/1‑140 (MQ=255)                                                    tATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGCTGAAGTGGGGCACTGCGTCGGGGTTTACGCTAACCCCGCCGATCTGAGCGCCAAGTCCCGGACGGCGATGCGCTTTAGCGAAACCg                                                                                        >  2:194879/1‑140 (MQ=255)                                                     aTCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGa                                                                                       >  2:121427/1‑140 (MQ=255)                                                     aTCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGa                                                                                       >  1:373518/1‑140 (MQ=255)                                                           tGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGaa                                                                                 >  1:804737/1‑140 (MQ=255)                                                             tAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACg                                                                               >  1:1137447/1‑140 (MQ=255)                                                             tAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACg                                                                               >  1:273605/1‑140 (MQ=255)                                                             aaaaGCTTGTCAGGTGGCGTGTTGGGAGTGGAACGACATCCGTGATGAAGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACg                                                                               <  2:245189/139‑1 (MQ=255)                                                                aGCCGGTCGGGTGGGGTGTTCGGAGTGGAACGACATCCGTGATGAAGGGGGGCACTGCGGGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCa                                                                            <  2:273605/140‑1 (MQ=255)                                                                    tGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCt                                                                        >  2:584640/1‑140 (MQ=255)                                                                       cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                                     >  2:156135/1‑140 (MQ=255)                                                                       cAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGa                                                                     >  1:1119118/1‑140 (MQ=255)                                                                           tGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTgg                                                                 >  2:336525/1‑140 (MQ=255)                                                                                  ttgggAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                          <  1:837981/137‑1 (MQ=255)                                                                                  ttgggAGTGGAACGACATCCGTGATGAAGGGGGGCAGGGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                          <  2:1074414/137‑1 (MQ=255)                                                                                  ttCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATcc                                                          >  1:548497/1‑140 (MQ=255)                                                                                           gAACGACATCCGTGATGAGGGGGGGCCCTGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                                 <  2:670280/140‑1 (MQ=255)                                                                                           gAACGACATCCGTGATGAGGGGGGGCACTGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                                 <  2:1137447/140‑1 (MQ=255)                                                                                           gAACGACATCCGTGATGAGGGGGGGCACTGCGTGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACg                                                 <  2:373518/140‑1 (MQ=255)                                                                                            aaCGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACgg                                                >  2:748547/1‑140 (MQ=255)                                                                                             aCGACATCCGTGATGAGGGGGGGCCCGGGGGGGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGc                                               <  2:1035832/140‑1 (MQ=255)                                                                                             aCGACATCCGTGATGAAGGGGGGCACTGCGTCGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGc                                               <  1:133075/140‑1 (MQ=255)                                                                                             aCGACATCCGGGAGGAGGGGGGGCACTGGGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGc                                               <  2:1119118/140‑1 (MQ=255)                                                                                                aCATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGt                                            <  1:658664/140‑1 (MQ=255)                                                                                                     cGTGATGAAGTGGGGCACTGCGTGGGGGTTTACGTTAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                         <  1:503403/138‑1 (MQ=255)                                                                                                     cGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGTTAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATg                                         >  2:503403/1‑138 (MQ=255)                                                                                                     cGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCa                                       <  1:898737/140‑1 (MQ=255)                                                                                                         atgaAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACtgt                                   >  1:566411/1‑140 (MQ=255)                                                                                                                 gggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATcc                           >  1:24938/1‑140 (MQ=255)                                                                                                                 gggCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATcc                           >  2:196829/1‑140 (MQ=255)                                                                                                                    cACTGCGTCGGGGTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                        <  1:657002/140‑1 (MQ=255)                                                                                                                    cACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                        <  1:196829/140‑1 (MQ=255)                                                                                                                    cACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                        <  1:584640/140‑1 (MQ=255)                                                                                                                    cACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAATCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                        >  1:839968/1‑140 (MQ=255)                                                                                                                    cACGGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                        <  2:857132/140‑1 (MQ=255)                                                                                                                    cACGGCGTGGGGTTTTACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                        <  2:752346/140‑1 (MQ=255)                                                                                                                       tGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTg                     >  2:569232/1‑140 (MQ=255)                                                                                                                           tCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                        <  1:1190267/133‑1 (MQ=255)                                                                                                                           tCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCgg                        >  2:1190267/1‑133 (MQ=255)                                                                                                                              gggTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAg              >  1:1181419/1‑140 (MQ=255)                                                                                                                                 tttACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCg           <  1:241866/140‑1 (MQ=255)                                                                                                                                   tACGATACCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCg                         <  1:600760/124‑1 (MQ=255)                                                                                                                                   tACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCg                         >  2:600760/1‑124 (MQ=255)                                                                                                                                        tAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGt    <  1:748547/140‑1 (MQ=255)                                                                                                                                          aCCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTct  >  1:75878/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                            GGTGCGTGATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGACATCCGTGATGAAGTGGGGCACTGCGTCGGGGTTTACGATAACCCCGCCGATCTGAGCGCCAAGTCCTGGACGGTGATGCGCTTTAGCGAAACCGAACAGAACGGCAAGCTGGAGTGGCTGATCCGTAAAGACGGCTGTATGCACTGTGAAGATCCCGGCTGCCTGAAGGCGTGCCCGTCT  >  NC_000913/1550526‑1550802

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