Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,404,215	0	T	G	25.5%	102.7 / 14.6	51	G244G (GGT→GGG)	yhdH	acrylyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (7/31);  new base G (13/0);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 1.63e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCA  >  NC_000913/3404077‑3404350                                                                                                                                           |                                                                                                                                        tgatgaGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGggtg                                                                                                                                        >  1:1055846/1‑140 (MQ=255)      aGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGgttt                                                                                                                                   >  2:364586/1‑140 (MQ=255)      aGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGTCTGGCGGGGGGgttt                                                                                                                                   >  2:1055407/1‑140 (MQ=255)      aGTTTGCCGAATCCCGTCCTCTGG‑AAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTa                                                                                                                                  <  2:176566/140‑1 (MQ=255)          tGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGGCACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTAccc                                                                                                                               >  1:335497/1‑138 (MQ=255)          tGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGTGGTCTGGCGGGTGGGTTTtctc                                                                                                                               >  2:334645/1‑140 (MQ=255)           gCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGGTTTActct                                                                                                                              >  1:913825/1‑140 (MQ=255)              gAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGcc                                                                                                                           <  1:1386848/140‑1 (MQ=255)                 tCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTATCTGCCAAc                                                                                                                        >  2:877672/1‑140 (MQ=255)                      tCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGGGGCGGCCTGGGGTCTGGCGGGGGGGTTTTCTCTGCCAACCCCgg                                                                                                                   >  1:1149637/1‑140 (MQ=255)                       ccTCTGGAAAAACAAGTCTGGGCTGGGGC‑ATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTc                                                                                                                 <  1:1398864/140‑1 (MQ=255)                        ctctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTc                                                                                                                 <  1:984437/140‑1 (MQ=255)                          ctGGAAAAACAAGTCTGGGCTGGGGCAATTGGCACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACCCTGCCAACCACGGTCAt                                                                                                               <  2:335497/140‑1 (MQ=255)                          ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                               <  1:334645/140‑1 (MQ=255)                          ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                               <  1:877672/140‑1 (MQ=255)                          ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                               <  2:21012/140‑1 (MQ=255)                          ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                               <  1:170106/140‑1 (MQ=255)                          ctGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCAt                                                                                                               <  1:1078474/140‑1 (MQ=255)                           tGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGAAAAAGGGCTGGCGCAAATGAATTACGGCGGCTGCGTGGGGGCCTGTGGTCTGGCGGGGGGTTTTTCTTTTCCCACCACGGTGAta                                                                                                              >  2:622130/1‑139 (MQ=255)                                    aaGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTAtt                                                                                                     <  2:114758/140‑1 (MQ=255)                                     aGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGACGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTc                                                                                                    <  1:1080509/140‑1 (MQ=255)                                        ctgggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGc                                                                                                 <  2:844857/140‑1 (MQ=255)                                         tgggctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCg                                                                                                <  2:1025545/140‑1 (MQ=255)                                            gctgggGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTaa                                                                                             <  2:1028724/140‑1 (MQ=255)                                                gggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTc                                                                                         >  1:445425/1‑140 (MQ=255)                                                 ggCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTcc                                                                                        <  2:758330/140‑1 (MQ=255)                                                     aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGGCATGCCATTTATTCTTCGTAATGTCCGttt                                                                                    >  2:891496/1‑140 (MQ=255)                                                     aTTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGGCATGCCCTTTATTCTGCGTAATGTCCGttt                                                                                    >  2:1204617/1‑140 (MQ=255)                                                         acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                <  2:913825/140‑1 (MQ=255)                                                         acacCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCaa                                                                                <  2:1259663/140‑1 (MQ=255)                                                              gTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGt                                                                           <  1:364586/140‑1 (MQ=255)                                                               ttGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTg                                                                          >  1:985721/1‑140 (MQ=255)                                                                 ggCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGGCTGGCGGGGGGGTTTTATCTGCCAACCACGGTCATGCCATTTATTTTTCGGAATGTCCGTTTGCAAGGGGTggg                                                                        >  2:732621/1‑139 (MQ=255)                                                                  gCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTGCCAACCACGGTCATGCCCTTTTTTTTTCGTAATGTCCGGTTGCCAGgggggggt                                                                       >  2:823894/1‑138 (MQ=255)                                                                       aaaGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGt                                                                  <  1:622130/140‑1 (MQ=255)                                                                         aGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGGCATGCCCTTTATTTTGCGTAATGTCCGTTTGCAAGGGGGGGGTTCAGTaa                                                                >  1:1196692/1‑140 (MQ=255)                                                                           tGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATg                                                              <  1:823894/140‑1 (MQ=255)                                                                                 cAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCCACCACAGGCATGGCATTTTTTTTTCGTAATGTCCCGTTTCATGGGGGGGATTTAGTAATGatcgcc                                                        >  2:121087/1‑135 (MQ=255)                                                                                       tGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGGGGTTTTACTCTGCCAACCACGGGCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGGGGATTCAGTAATGACGCCACCAGaa                                                  >  1:64879/1‑140 (MQ=255)                                                                                              gcAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCcgcg                                           <  2:343020/140‑1 (MQ=255)                                                                                                aaaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGca                                         <  1:1409968/140‑1 (MQ=255)                                                                                                 aaTGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGcac                                        <  1:891496/140‑1 (MQ=255)                                                                                                     aaTTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGc                                    <  2:985721/140‑1 (MQ=255)                                                                                                            gcggcTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTTCTCTTCCAACCACGGGCATGCCCTTTATTTTGCGGAATGTCCGTTTGCAAGGGGGGGGTTTAGTAATGACGCCACCAGAACGCCGCGCACAAACCTGGCAg                             >  2:1289417/1‑140 (MQ=255)                                                                                                            gcggcTGCGTGGCGGCCTGTGGTCTGGCGGGGGGGTTTACTCTGCCAACCACGGGCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGGGGATTCAGTAATGACGCCCCCAGAACGCCGCGCACAAGCCTGGCAg                             >  1:1178793/1‑140 (MQ=255)                                                                                                             cggcTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGGGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGc                            >  1:518304/1‑140 (MQ=255)                                                                                                                          gcCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGa               >  1:631464/1‑140 (MQ=255)                                                                                                                                    tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa     <  2:518304/140‑1 (MQ=255)                                                                                                                                    tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa     <  2:1281795/140‑1 (MQ=255)                                                                                                                                    tGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGaa     >  2:100569/1‑140 (MQ=255)                                                                                                                                     ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt    >  1:185289/1‑140 (MQ=255)                                                                                                                                     ggCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAAt    <  2:1178793/140‑1 (MQ=255)                                                                                                                                      gCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc   <  2:1166647/140‑1 (MQ=255)                                                                                                                                      gCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATc   <  2:1055846/140‑1 (MQ=255)                                                                                                                                       cGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCa  <  1:420146/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        TGATGAGTTTGCCGAATCCCGTCCTCTGGAAAAACAAGTCTGGGCTGGGGCAATTGACACCGTTGGCGACAAAGTGCTGGCAAAAGTGCTGGCGCAAATGAATTACGGCGGCTGCGTGGCGGCCTGTGGTCTGGCGGGTGGTTTTACTCTGCCAACCACGGTCATGCCATTTATTCTGCGTAATGTCCGTTTGCAAGGGGTGGATTCAGTAATGACGCCACCAGAACGCCGCGCACAAGCCTGGCAGCGACTGGTCGCCGATTTACCGGAATCA  >  NC_000913/3404077‑3404350

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