Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,740,909	0	T	G	26.5%	89.9 / 10.1	49	V402G (GTG→GGG)	avtA	valine‑‑pyruvate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (11/25);  new base G (13/0);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 9.51e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGCCTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATA  >  NC_000913/3740774‑3741039                                                                                                                                        |                                                                                                                                   aGCGCCTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAg                                                                                                                                >  1:405924/1‑140 (MQ=255)    gcCTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGAtt                                                                                                                             <  2:405924/140‑1 (MQ=255)    gcCTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGAtt                                                                                                                             >  1:652359/1‑140 (MQ=255)     ccTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTc                                                                                                                            >  2:520819/1‑140 (MQ=255)       tGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTg                                                                                                                          >  1:471675/1‑140 (MQ=255)           aGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCgg                                                                                                                      <  2:341117/140‑1 (MQ=255)           aGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCgg                                                                                                                      <  1:470248/140‑1 (MQ=255)             cACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGaa                                                                                                                    <  2:969386/140‑1 (MQ=255)             cACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGaa                                                                                                                    >  1:764231/1‑140 (MQ=255)                gcggcgTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAagag                                                                                                                 >  1:591264/1‑140 (MQ=255)                gcggcgTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAagag                                                                                                                 >  2:19031/1‑140 (MQ=255)                 cggcgTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAgaga                                                                                                                <  2:309616/140‑1 (MQ=255)                    cgTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAg                                                                                                             >  2:945603/1‑140 (MQ=255)                           aTGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGc                                                                                                      <  2:764231/140‑1 (MQ=255)                             ggTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCt                                                                                                    <  1:471066/140‑1 (MQ=255)                              gTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATAGAAAGAGCCTg                                                                                                   >  2:41896/1‑140 (MQ=255)                                 ccGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGc                                                                                                >  1:753427/1‑140 (MQ=255)                                  cGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCt                                                                                               >  1:696724/1‑140 (MQ=255)                                      gCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGaaa                                                                                           <  2:807257/140‑1 (MQ=255)                                         caACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTTTGGCGGAAGAGATAGAAAGAGGCTGGGCTGAAAggc                                                                                        >  2:471066/1‑138 (MQ=255)                                               cttcCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaaa                                                                                  >  2:924098/1‑140 (MQ=255)                                                ttcCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaaaa                                                                                 >  2:424277/1‑140 (MQ=255)                                                     aGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGc                                                                            <  1:984848/140‑1 (MQ=255)                                                       ggCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCt                                                                          >  2:885454/1‑140 (MQ=255)                                                        gCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTg                                                                         <  2:471675/140‑1 (MQ=255)                                                           ggATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATg                                                                      <  1:424277/140‑1 (MQ=255)                                                            gATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgc                                                                     >  2:541197/1‑140 (MQ=255)                                                            gATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCCTGGGTGAAAATCACTAAAAATTGCCTGAttc                                                                     >  2:984848/1‑138 (MQ=255)                                                             aTAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcg                                                                    <  1:520819/140‑1 (MQ=255)                                                             aTAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcg                                                                    <  2:892354/140‑1 (MQ=255)                                                             aTAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcg                                                                    <  2:696724/140‑1 (MQ=255)                                                                aaCCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTa                                                                 <  1:19031/140‑1 (MQ=255)                                                                  ccGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACg                                                               <  1:677327/140‑1 (MQ=255)                                                                  ccGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACg                                                               <  1:924098/140‑1 (MQ=255)                                                                      ggCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGTGCTGAAAGTCACTAAAAATTGCCTGATgcg                                                                    <  1:64032/131‑1 (MQ=255)                                                                      ggCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGTGCTGAAAGTCACTAAAAATTGCCTGATgcg                                                                    >  2:64032/1‑131 (MQ=255)                                                                      ggCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTa                                                           >  1:121968/1‑140 (MQ=255)                                                                        ccGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATc                                                         <  2:753427/140‑1 (MQ=255)                                                                        ccGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATc                                                         <  2:121968/140‑1 (MQ=255)                                                                                       tGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAAATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATTCGCTCCGCTTATCAGGCCCACATGGttt                                          >  2:287513/1‑140 (MQ=255)                                                                                            gcgcATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCa                                     <  1:945603/140‑1 (MQ=255)                                                                                            gcgcATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCa                                     <  2:51992/140‑1 (MQ=255)                                                                                            gcgcATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCa                                     >  1:980868/1‑140 (MQ=255)                                                                                               cATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAata                                  >  2:187570/1‑140 (MQ=255)                                                                                                        cGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCt                                        >  1:27061/1‑125 (MQ=255)                                                                                                        cGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCt                                        <  2:27061/125‑1 (MQ=255)                                                                                                        cGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAAtt                         >  2:471617/1‑140 (MQ=255)                                                                                                        cGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAAtt                         >  2:301782/1‑140 (MQ=255)                                                                                                          tACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTg                       <  2:652359/140‑1 (MQ=255)                                                                                                                agCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGgaagaa                 >  2:83626/1‑140 (MQ=255)                                                                                                                agCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATAAAAAAAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGgaagaa                 >  2:960599/1‑140 (MQ=255)                                                                                                                     gagaAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAAAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGt            >  2:675878/1‑140 (MQ=255)                                                                                                                         aaaTTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGc        <  2:980868/140‑1 (MQ=255)                                                                                                                               aGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATa  >  2:956018/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                   AGCGCCTGAAAGCACGCGGCGTGCTGATGGTGCCGGGGCACAACTTCTTCCCAGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATTTCTGCAATATATTGAATTTGGAAGAATTTGTAGGCCGGATA  >  NC_000913/3740774‑3741039

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: ‑