Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,480,380	0	A	G	46.0%	52.4 / 14.4	50	G393G (GGA→GGG)	dsdA	D‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base A (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

CGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATCTGGCAAAAGGCCGTT  >  NC_000913/2480231‑2480528                                                                                                                                                       |                                                                                                                                                     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)                                                                                                                                                    cGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATCTGGCAAAAGGCCGt   <  1:683447/151‑1 (MQ=255)                                                                                                                                                      gACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATCTGGCAAAAGGCCGtt  >  2:172165/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                     CGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAT‑GGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGCACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATCTGGCAAAAGGCCGTT  >  NC_000913/2480231‑2480528

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