Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,355,457	0	A	C	46.2%	47.7 / 13.0	39	T245P (ACC→CCC)	gltB	glutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base A (16/5);  new base C (0/18);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 9.98e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGCGCTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGG  >  NC_000913/3355308‑3355606                                                                                                                                                      |                                                                                                                                                      tGCCGCGCTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCAc                                                                                                                                                      >  2:111789/1‑151 (MQ=255)    cgcgCTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCgg                                                                                                                                                   <  1:564669/151‑1 (MQ=255)        cTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAAc                                                                                                                                               >  2:328653/1‑151 (MQ=255)          tttATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAAcc                                                                                                                                              <  1:45903/150‑1 (MQ=255)                        gCGGCCCTGCGTCTGGAATCGCCCTTTTGCCTTTCCCCCCAGCGCTTCTCCACTACCCCCGTCCCGCGCGGGCCGCGGGCCCACCCTTCCCGCTATCGGGCGCATAACGGTAAAATCACCCCCATCCCCGGTAACCGCCAATGGgcgc                                                                                                                                  <  1:93574/148‑1 (MQ=255)                            aCCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg                                                                                                                           >  1:382107/1‑151 (MQ=255)                               tGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTAc                                                                                                                        >  2:498034/1‑151 (MQ=255)                                gCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACt                                                                                                                       >  2:354472/1‑151 (MQ=255)                                 cGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACtt                                                                                                                      >  1:839075/1‑151 (MQ=255)                                      ggAATCGGCCTTTTGCCTGTTCCACCCGCCCTTCTCCACTAACCCCGTCCCGCGCTGGCCGCTGGCCCAACCTTTCCGCTGTCTGGCCCATAGGGGGGAAATCAACCCCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTtat                                                                                                                    <  2:303343/148‑1 (MQ=255)                                        aaTCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAAtt                                                                                                               >  2:442802/1‑151 (MQ=255)                                          tCGCCCTTTTGCCTGTCCCCCCAGCGCTTCCCCCCAACCCCCGTCCCGCGCTGCCCGCTGGCGCACCCTTCCCGCTATCGGGCGCATAACGGGAAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATaaa                                                                                                                 <  2:41455/147‑1 (MQ=255)                                           cGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTcc                                                                                                             >  2:608923/1‑150 (MQ=255)                                              ccttttGCCTTTCCCCCCGCCGCTTCCCCACTACCCCCGTCCCGCGCTGGCCGCGGGCGCACCCGTCCCGCTTTCTGGCGCATAACGGTAAAATCAACCCCATCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTcc                                                                                                             <  2:752493/144‑1 (MQ=255)                                                 tttGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACAcc                                                                                                      >  1:32888/1‑151 (MQ=255)                                                 tttGCCTGTCCCACCAGCGCTTCCCCACTAACACCGTCCCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCTTAACGGTGAAACCACCACCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac                                                                                                       <  2:765969/150‑1 (MQ=255)                                                   tGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGc                                                                                                    >  2:684548/1‑151 (MQ=255)                                                     ccTGTCCCACCCGCTCTTCCCCACAACCCCCGTCCCGCGCTGCCCGCTGGCCCACCCGTCCCGCTATCGGGCGAATAAGGGTGAAACAAACACCATCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCt                                                                                                   <  2:852465/150‑1 (MQ=255)                                                           cccccggcgtTTCCCCACAACCCCCGTCCCGCGCTGGCCGTGGGCCAACCCGTCCCGCTATCTGGCGCATAACGGTGAAATCAACCCCCCCCCGGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGc                                                                                                    <  2:226676/137‑1 (MQ=255)                                                            caccaGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTg                                                                                           >  1:24311/1‑151 (MQ=255)                                                             acccGCGCTTCTCCACGAACCCCGTCCCGCGTTTGCCGCGGGCGAACCCTTCCCGCGTTCTGGCGAATAAAATAAAATTCAGCACCCCCACGTGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATcc                                                                                              <  2:547004/147‑1 (MQ=255)                                                                  cgcTTCCCCACAACCCCCGTCCCGCGCGGGCCGCTGGCGCACCCTTCCCGCTATCTGGCGCATAAGGGTAAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGc                                                                                     <  2:636660/151‑1 (MQ=255)                                                                   gcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGc                                                                                     >  2:186508/1‑150 (MQ=255)                                                                     ttCTCCACAACCCCCGTCCCGCGCTGGCCGCTGGCGCAACCTTCCCGTTTTCGGGCGCATAACGGTAAAATAAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCt                                                                                       <  2:123224/146‑1 (MQ=255)                                                                       ctcCACAACCCCCGTCCCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATACGGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCacg                                                                                  <  2:206552/149‑1 (MQ=255)                                                                                 acCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGt                                                                      >  2:80257/1‑151 (MQ=255)                                                                                              gcccgcgggcGCACCCGTTCCGCTATCTGGCGCATAACGGTAAAACAACCCCCATCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAAc                                                              <  2:168413/144‑1 (MQ=255)                                                                                                 cgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTCACCGCCAATGGGCGCGAGCGCGTACTTATACATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCgg                                                      >  1:312739/1‑151 (MQ=255)                                                                                                            ccGTCCCGCTTTCTGGCGAATACCGGTAAAACCACCACCACCCCCGTAACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTc                                             <  1:147307/149‑1 (MQ=255)                                                                                                            ccGTCCCGCTATCTGGCGCATAACGGTAAAATCAACCCCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCa                                           <  2:9481/151‑1 (MQ=255)                                                                                                             cGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTcc                                            >  2:115849/1‑149 (MQ=255)                                                                                                               tccgcctattggcggcATAGCGGAAAAATCCCCCCCACCCCGGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGt                                         <  1:201726/136‑1 (MQ=255)                                                                                                                 ccGCTTTCGGGCGCATAACGGTAAAACAAACCCCACCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTc                                       <  1:236262/150‑1 (MQ=255)                                                                                                                  ctcTATCTGGCGCAAAACGGTAAAACCACCACCACCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTc                                       <  2:189872/147‑1 (MQ=255)                                                                                                                                         aaTCAACCCCATCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAActgctg              <  2:285886/151‑1 (MQ=255)                                                                                                                                          accacccccATCCCCCGCAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAActgctg              <  2:506236/143‑1 (MQ=255)                                                                                                                                                 ccATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAgg      <  2:248730/151‑1 (MQ=255)                                                                                                                                                 ccATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAg       <  2:855838/150‑1 (MQ=255)                                                                                                                                                     cACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAActgctgct            >  2:837033/1‑141 (MQ=255)                                                                                                                                                      aCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCggg  >  2:315909/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TGCCGCGCTTTTATCTGGATCTTGCGGACCTGCGTCTGGAATCGGCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGG  >  NC_000913/3355308‑3355606

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