Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,480,371	0	T	G	52.1%	52.6 / 73.1	48	G390G (GGT→GGG)	dsdA	D‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base T (4/19);  new base G (25/0);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 7.67e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-07

TGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATA  >  NC_000913/2480227‑2480511                                                                                                                                                 |                                                                                                                                              tGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGc                                                                                                                                          <  2:365643/150‑1 (MQ=255) tGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATgg                                                                                                                                           >  1:701319/1‑149 (MQ=255)   gTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGcc                                                                                                                                         <  2:585251/149‑1 (MQ=255)     cgcgCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGGTGGcc                                                                                                                                         >  2:438132/1‑147 (MQ=255)       cgcATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCgg                                                                                                                                       >  2:473674/1‑147 (MQ=255)            cAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAgc                                                                                                                               <  1:293545/150‑1 (MQ=255)            cAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATTGGCgggcatc                                                                                                                                  >  1:170680/1‑142 (MQ=255)                    gTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGCGCTAGCCGGATCTCatctcccgtg                                                                                                                        >  1:540800/1‑140 (MQ=255)                       gggcgggcAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGGACCCCGCGCgggtgggggt                                                                                                                   >  2:343886/1‑146 (MQ=255)                         gcgggcAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCa                                                                                                                 <  1:651526/151‑1 (MQ=255)                           gggcAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGCCGGACCTCCGCGCGGGGGTGGAt                                                                                                                   >  1:106403/1‑147 (MQ=255)                               aaTGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGGCCTCAGCGCGGGGGTGCATCAg                                                                                                                >  2:25500/1‑146 (MQ=255)                                       gTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGGCCCCTGCGCGGGTGGGGGTTAGTAAGGTACCaa                                                                                                     >  2:563376/1‑149 (MQ=255)                                           gctgGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGa                                                                                                 <  2:675704/149‑1 (MQ=255)                                                aTGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGAACCCCGCGCGGGGGGGCGTTAGGTAATTACCAACACATA‑CACgg                                                                                         >  2:668374/1‑151 (MQ=255)                                                     ttCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGCCGGGGCCCACCGCGGGTGTGTATCCGGAAAGTAACACCAAATG‑GACGGgtt                                                                                      >  2:822878/1‑149 (MQ=255)                                                           aCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAg                                                                              <  2:845031/151‑1 (MQ=255)                                                           aCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAg                                                                              <  2:116687/151‑1 (MQ=255)                                                            cccTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAg                                                                              <  1:20597/150‑1 (MQ=255)                                                              cTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGCGGGCCGGGCCTCCGCGCGTGGGTGCGTCAGTAAGTTACCAAACGATG‑GACGGGTTCAGCGCAGaa                                                                            >  2:114143/1‑150 (MQ=255)                                                               ttAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGGCCTCCGCGCGGGTGGGCATCAGTAAGGTACCAACAGATG‑CACGGGTTTAGCGCAGaa                                                                            >  1:480681/1‑149 (MQ=255)                                                                    gATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGCCGGAACTCAACGCGGGCGGTCATCTATAAGGTACCAACAGAAGTCACGGttccgggcagaaacaa                                                                         >  2:507611/1‑134 (MQ=255)                                                                       cAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGCCGGACCCCCACGCGGGTGTGGATCCGTAAGTTACCCACCAAAT‑CAAGGGTTTAGCGCAGAAAAAATg                                                                      >  2:754037/1‑147 (MQ=255)                                                                        aaaCCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGt                                                                   <  2:797120/149‑1 (MQ=255)                                                                           ccATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATAc                                                              <  1:352177/151‑1 (MQ=255)                                                                           ccATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGGCCTCAGCGCGCGGGTGCATCCGTAAATTAACAAAAAAAG‑CACGGGTTTAGCGCAGAACAAACGCGGAATa                                                               >  2:82667/1‑150 (MQ=255)                                                                                 aTGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGGCCCCCGCCCGGGTGTGCGTCAGGAAGGTACCAACAAATG‑CACGGGTTTCGCGCCGAAAAAATCCGGAATaaccatc                                                         >  2:641642/1‑144 (MQ=255)                                                                                   gACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGCCGGGCCTCAACGCGTGGGTGGATCAGTAAGTTAACAACAAAAG‑CACGGGTTTAGCGCAGAAAAAATGCGGAATa                                                               >  1:904661/1‑142 (MQ=255)                                                                                      aTGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTAAGGACGGACCCCAGCGCCGGCGTGCATCATGAAAATACCCACAGATG‑CACCGTTTCCGCGCCAAACAACAGCGCCATAcccatca                                                        >  1:885526/1‑141 (MQ=255)                                                                                          ttGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGtgtg                                               >  2:213777/1‑151 (MQ=255)                                                                                            ggctggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTgg                                              <  1:286214/150‑1 (MQ=255)                                                                                               tggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGCCGGGCCTCAGCGCGGGGGTGGGTCAGGAAGGTAACCACAAATG‑CACGGGTTTAGCGCAGAAAAACTGCGTAAAAACACCCATCTGGTGgggg                                             >  1:22348/1‑148 (MQ=255)                                                                                                ggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCg                                           <  2:755965/149‑1 (MQ=255)                                                                                                 gctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGACCTCAGCGCGGGTGTGCATCAGTAAGGTACCAACAGATG‑CACGGGTTCAGCGCAGAACAACTGCGGAATACCCCTCATCTGGGGTGGGCGCCg                                        >  1:74729/1‑151 (MQ=255)                                                                                                    ggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACggg                                      <  1:801849/150‑1 (MQ=255)                                                                                                      cgcAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAACACCACTCATCTGGTGTGGGCGACGGGAgg                                   <  2:872530/151‑1 (MQ=255)                                                                                                             agaagGTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGCCGGACCCCCGCGCGTGTGTGGATCCGTAAGTTACCAAAAGATG‑CCCGGGTTCTGCGCAGAACAACTGCGTCATACCACTCAACTGGGGTGGGCGACGGGGGGGGGaa                              >  1:680681/1‑149 (MQ=255)                                                                                                                agGTATTCGTCTTGAACCTTCGGCACTGGCGGGGTTGGGCGGGCCTCAGCGCGTGTGTGGCTTAGTAAGGTACCAACAGAAG‑CACGGGTTCAGGGCAGAAAAAACGGGGAAAAACCCCCATCTGGTGGGGGGGGCGGGGGGGGGGATgg                           >  2:154016/1‑149 (MQ=255)                                                                                                                 ggTATTCGTCTTGAACCTTCGGCACTGGCGGGGATGGGCGGGCCCCCGCGCGTGTGTGGATCCGTAAGTTACCCACAAATG‑GACGGGTTTCGCGCAGAACAAATGCGTAATACCACTCCTCTGGTGTGGGGGGCGgggggggggaagg                           >  1:484079/1‑143 (MQ=255)                                                                                                                  gTATTCGTCTTGAACCTTCGGCACTGGCGGGGATTGCCGGGCCCCACCGCGGGTCTGCCACACTAAATTAACAAAAATCT‑CACCGGTTCAACGGCAAAAAACCGCGCAATAGCACTTCACTGGGGTGGGCGGCCGGGGGTGGGACGGGGcc                       >  1:442127/1‑151 (MQ=255)                                                                                                                   tATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGGCCTCCGCGCGGGGGGGGATCAGTAAGGTAACAACAAATG‑CACGGGTTTAGCGCAGAACAAACGCGTAATACCACCCCTCTTGTGGGGGGGACGGGGGGGGGAATGGGGCCg                      >  1:846153/1‑151 (MQ=255)                                                                                                                       cGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGgaag                  <  2:251553/151‑1 (MQ=255)                                                                                                                        gTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGgaaga                 <  2:592724/151‑1 (MQ=255)                                                                                                                         tCTTGAACCTTCGGCACTGGCGGGTATTGCCCGGACTCAGCGCGTGTGTGGCTCAGTAAGTTACCCACAGATG‑CACGGGTTTAGCGCAGGACAACAGCGGAATATCCCTCATCTGGGGTGGGGGACGgggggggggaagggggcgggag                  >  1:506820/1‑135 (MQ=255)                                                                                                                            tGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGaagaag               <  1:292588/149‑1 (MQ=255)                                                                                                                                   tCGGCACTGGCGGGGATGGCCGGACCCCAGCGCGGGTGTGCATCAGTAAGTTACCAACAGATG‑GACGGGTTTAGCGCAGAACAAATGCGGAATAACACTCATCTGGTGTGGGGGACGGGGGGGGGGATGGTGGCGGGAGAAGAGATGaa        >  2:315420/1‑149 (MQ=255)                                                                                                                                    cGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATc      <  1:563395/150‑1 (MQ=255)                                                                                                                                       cACTGGCGGGTATGGCCGGACCTCAGCGTGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAata  <  2:479706/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                              TGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATG‑CACGGTTTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATA  >  NC_000913/2480227‑2480511

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