Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,480,371	0	T	G	46.0%	62.6 / 12.7	50	G390G (GGT→GGG)	dsdA	D‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base T (6/21);  new base G (23/0);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 5.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATC  >  NC_000913/2480229‑2480513                                                                                                                                               |                                                                                                                                                gTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑gggccg                                                                                                                                        >  2:128168/1‑147 (MQ=255)   cgcgCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGa                                                                                                                                       <  2:119281/150‑1 (MQ=255)          cAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGCCGGGCCTCCgc                                                                                                                                >  2:707468/1‑150 (MQ=255)           aGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAgcgc                                                                                                                              <  1:98039/151‑1 (MQ=255)             gCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGGCCCCAGcgcg                                                                                                                             >  2:67562/1‑150 (MQ=255)                         gtgcAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGt                                                                                                                <  1:571474/149‑1 (MQ=255)                           gcAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTaa                                                                                                              <  2:377946/151‑1 (MQ=255)                              aTGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGGCGGGCCCCAGCGCGGGTGTGGGTCAGGAAAtt                                                                                                           >  2:22807/1‑151 (MQ=255)                               tGGAGCGTCTGCTGGATGGCTTCTATACCCTTTGCGATCCAATCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGCACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGt                                                                                                            <  1:482192/149‑1 (MQ=255)                               tGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGGCCCCCACGCGTGTGTGCATCAGTa                                                                                                               >  2:843972/1‑146 (MQ=255)                                ggAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTAc                                                                                                         <  2:267524/151‑1 (MQ=255)                                    cGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAG‑TGCCGGGCCCCCGCGCGTGTGGGGATCAGGAAGGTACCaa                                                                                                      >  2:768203/1‑150 (MQ=255)                                     gTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTTCCaa                                                                                                      <  1:732192/149‑1 (MQ=255)                                     gTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGGCGGAGCCCACCGCGGGTGTGCCTCCGttaagtaccca                                                                                                      >  1:758201/1‑139 (MQ=255)                                         gctgGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAg                                                                                                   <  1:361553/148‑1 (MQ=255)                                         gctgGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGAt                                                                                                 <  1:332421/150‑1 (MQ=255)                                           tgGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑TGGCGGGCCCCCGCGCGGGTGTGGATCAGTAAGTTAACAACAAATGc                                                                                               >  2:762528/1‑150 (MQ=255)                                            ggATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTCT‑GGGCGGGCCGCAGCGCGGGGGTGGCGCAGAAAGGTACCAAACGATGGAc                                                                                             >  1:180876/1‑151 (MQ=255)                                              aTGGCTTCTATACCCGTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACgg                                                                                           <  1:823330/151‑1 (MQ=255)                                                   ttCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGGCGGAACCCCGCGCGGGGGTGGCTCCGTAAGTTACCAACAAAAGCACggggt                                                                                        >  1:668722/1‑146 (MQ=255)                                                       ataCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAgcg                                                                                   <  2:766479/150‑1 (MQ=255)                                                       ataCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAA‑TGGCGGGCCTCCACGCGTGTGTGCGTCAGCAAAATAACCACAAAATGACGGTTTCAg                                                                                     >  1:766422/1‑148 (MQ=255)                                                              tAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑TGGCGGAACCCAGCGCGGGTGGGCATCCGTAAAGTACCCACAAATGCACGGGTTTCGCGCCGaa                                                                              >  2:687053/1‑148 (MQ=255)                                                                   aTCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGGT‑GGGCGGGGCTCAGCGCGGGTGTGCCTCCGGAAGTTTACCACAGATACACGGGTTTCGCGCAGaaaaa                                                                           >  1:267656/1‑146 (MQ=255)                                                                   aTCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGGCGGGCCTCAGCGGGGGGGTGGGTCATGAAGGTACCAAACAAGGCACGGGTTTAGCGCAGaaaaa                                                                           >  2:43758/1‑146 (MQ=255)                                                                   aTCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGCCGGGCCCCCGCGGCGGGGTGCATCAGTAAATTAACCACAGAAGCACGGGTTTAGCGCAGAAAAAATGc                                                                       >  1:531628/1‑150 (MQ=255)                                                                    tCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCCGGGCC‑CAGCGCGTGTGTGGAACAGAAAATTACCCAAAGATGGACGGGTTTTGCGCAGAACAAATGCg                                                                      >  2:585053/1‑150 (MQ=255)                                                                    tCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGGCCGGCCCCAGCGCGGGGGTGGCACAGAAAGTTACCAACAGATGCGCGGGTTTAGGGCAGaaaaa                                                                           >  2:357246/1‑145 (MQ=255)                                                                      aaaCCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGGCGGGCCCCCGCGCGTGTGTGCATCAGTAAGTTACCAACCGATGCACGGGTTTAGCGCAGAAAAACAGCGTa                                                                    >  1:281585/1‑150 (MQ=255)                                                                          cATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGCCGGACCCCAGGGGGTGGGTGCATTAGTAAGTTAACAACAGAAGCACGGGTTCCGCGCCGAAAAAATGAGTAATAcc                                                               >  2:237167/1‑151 (MQ=255)                                                                           aTGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCa                                                              <  1:762471/151‑1 (MQ=255)                                                                             gTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCAc                                                             <  2:101937/150‑1 (MQ=255)                                                                                tGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAGCCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTc                                                           <  1:323129/149‑1 (MQ=255)                                                                                  aCATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCt                                                       <  1:743162/151‑1 (MQ=255)                                                                                       cTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGGCGGGCCCCAGCGCGTGTGGGCATCAGTAAGTTACCAACAAAAGCACGGGTTTAGCGCAGAACAACTGCGGAATAACACCCCTCTGgtg                                                   >  2:292310/1‑150 (MQ=255)                                                                                       cTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGGT‑GGGCGGGGCCCCGCGCGGGGGTGCCTCAGTAAGTTACCAACAGATGGACGGGTTTCGCGCAGAAAAAATGCGTAATAACACCCCTCTGgtg                                                   >  1:10191/1‑150 (MQ=255)                                                                                        ttGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGGCCTCCGCGCGTGTGTGCATCAGTAAGATACCAACAAAAGCAACGGTTCAGCGCCGAACAACTGCGGAAAACCACTCATCTGgggt                                                  >  2:214674/1‑150 (MQ=255)                                                                                         tggctggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAA‑GGGCGGGCCCCCGGGCGCGTGTGCCTTACTAAAGTACCCACACAAGTCCGGGTTCAGCGCCGAAAAAATGcccaaaacc                                                               >  1:178950/1‑128 (MQ=255)                                                                                          ggctggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTggg                                               <  1:704421/151‑1 (MQ=255)                                                                                           gctggctggcGCAGGAAGAAGGTATTCGTCTTGATCCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATCCCACTCATCTGGTGTggg                                               <  1:737280/150‑1 (MQ=255)                                                                                              ggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGCCGGGCCCCAGCGCGGGTGTGCATCAGTAAGTTACCAACAAAAGCACGGGTTTAGCGCAGAACAAAAGCGGAATACCACTCATCTGGGGTggggg                                             >  2:696739/1‑147 (MQ=255)                                                                                                ctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTGT‑TGCCGGGCCCCAGCGCGGGTGTGGATTCGTAAAGTAACAACAAATTCACGGGTTTAGCGGAGAACAAATGCCGAATACCAATCCTCCGGTGTGGGGGGCgg                                         >  2:266537/1‑151 (MQ=255)                                                                                                           agaagGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAT‑GGGCGGGCCTCAGCGCGTGTGTGGATCAGTAAGGTACCAAAAGATGGACGGGTTCAGCGCCGAAAAAATGCGGAATAACACTCATCCGGTGTGGGGGGCGgggggggg                                  >  2:391946/1‑147 (MQ=255)                                                                                                            gaagGTATTCGTCTTGAACCTTCGGCACTGGCGGGGAA‑GGGCGGGCCTCAGCGCGTGGGGGCAGAAGTAAGATAACAACAAATGCACGGGTTTAGGGGAGAACAAAAGCGTAATAACACTCATCTGGGGGGGGGGAAGgggggggg                                  >  2:329850/1‑146 (MQ=255)                                                                                                                     cGTCTTGAACCTTCGGCACTGGCGGGGAT‑TGCCGGGCCTCCGCGCGGGTGTGCATCAGTAAAATACCAACAAAAGGACGGGTTTAGCGCAGAAAAAATGCGTAATAACACTCATCTTGTGGGGGGGGCGGGGGGGGGGAAGGGGCCgg                       >  2:700675/1‑148 (MQ=255)                                                                                                                          tGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAagaaga                <  2:306244/150‑1 (MQ=255)                                                                                                                          tGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAagaaga                <  1:40220/150‑1 (MQ=255)                                                                                                                          tGAACCTTCGGCACTGGCGGGGAA‑TGGCGGACCTCCCCGCGTGTGTGCCTCAGGCAGGTACCAAAAAACGACCGGGTTCTGCGCACAACCACCACGGAATAATACTCCTCTGGTGTGGGCGGCGggggggggggaggggccgg                       >  1:547962/1‑132 (MQ=255)                                                                                                                                      aCTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAatat   <  1:825138/151‑1 (MQ=255)                                                                                                                                         ggCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATc  <  2:359656/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                GTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATC  >  NC_000913/2480229‑2480513

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