Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,177	0	A	C	42.5%	63.6 / 11.6	41	G308G (GGT→GGG)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base A (20/3);  new base C (0/17);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 2.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACT  >  NC_000913/2921030‑2921319                                                                                                                                                    |                                                                                                                                               cgcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc                                                                                                                                             <  2:1194221/151‑1 (MQ=255) cgcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc                                                                                                                                             <  2:445901/151‑1 (MQ=255)  gcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc                                                                                                                                             >  2:176088/1‑150 (MQ=255)   ccATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAg                                                                                                                                           <  1:112573/151‑1 (MQ=255)     aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCGCCCAGca                                                                                                                                         <  1:1009991/151‑1 (MQ=255)     aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGc                                                                                                                                          >  2:361725/1‑150 (MQ=255)       aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGcac                                                                                                                                        >  2:1212429/1‑150 (MQ=255)           ccGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAcc                                                                                                                                     <  2:426210/149‑1 (MQ=255)                         caaaaCTTTGACGTAGTTGCAGACCCCCATCACCTTCGACACCAACATGCCCTTTCCGATCGGTTTTTTACGCGCAATGTCCAGCGATCCCGTGCGGCGCATCGTCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAaaaccgca                                                                                                                      <  2:893806/147‑1 (MQ=255)                                ttGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAAccgccg                                                                                                              >  1:197785/1‑151 (MQ=255)                                         ttGCAGAACCCCATCCCCATCGACAGCACCATCCCCATTACGATCGGTGTTTTCCGCGAAATTTTCAGCGACCCCGTGGGGCACATCACCCAGTCGGAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                         <  2:1249360/147‑1 (MQ=255)                                           gCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGcc                                                                                                   >  2:350571/1‑151 (MQ=255)                                             aGAACCCCACCACCATCGACAGCAACATGCCCTTTAGGATCGGTGTTTTACGCGCAATTTTCAGCGACCCCGTGGGGCGCATCACCGAGTGGAAAAAAACCCCACCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGcc                                                                                                   <  2:220419/149‑1 (MQ=255)                                               aaCCCCACCCCCATCAACAACAACACCCCCATTACGATGGGTTTTTTGCGCGAAATGTCCAGCAACCCCGGGCGGCGCACCACCAAGTCGAAATTAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGc                                                                                                    <  2:709891/146‑1 (MQ=255)                                                           atcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                    >  2:1300176/1‑150 (MQ=255)                                                             cGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                    >  2:1062824/1‑148 (MQ=255)                                                              gACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATCAAACCGCACACCGCCGGAACGGAAGCCACAAAGCCCGCttt                                                                                     >  1:685966/1‑146 (MQ=255)                                                                caccaccATCCCAATTACGATCGTTGTTTCACGCGCAATGTTCAGCGCCCCCGTGCGCCGCATCAGCCAGCCGAAAAAACCCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCa                                                                                   <  1:973826/140‑1 (MQ=255)                                                                      ctgccccATAAGGTTCGGTGTTTTACGCGAAATGTCCAGCGACCCCGGGCGGCCCACCAGCCAGTCGAAAAAAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGa                                                                            <  2:1048226/143‑1 (MQ=255)                                                                          ccaaTTACGAGCTGTTTGTTACGCGAAAGGTCCAGCGCCCCCGGGCGGCGCATCAGCCACTCGAAAATACTCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                    <  1:347784/148‑1 (MQ=255)                                                                                 cgtgtggTTTTTTACGCGAAATTTCCAGCCACCCCGGGCGGCCCATCAGCCAGCCGAAAAAAACCCCCCCCAGAACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                              <  2:1267627/145‑1 (MQ=255)                                                                                  gATGGGTTTTTTACGCGCAATGTCCAGCGTCCCCGGGCGGCCCATCAGCCGGTCGGAAAAAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTg                                                             <  2:721164/150‑1 (MQ=255)                                                                                        tttttccggcgaAATTTTCAGCGCCCCCGTGGGGCGCTGAAGCGAGTGGAAAATAACCCCCCCAAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCa                                                        <  2:664836/141‑1 (MQ=255)                                                                                         gtTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                         >  1:1162502/1‑147 (MQ=255)                                                                                          ttttCCGCGCATTGTCCGGCACCCCCGCGCGGCAAAACGGCCAGTGGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGa                                                      <  1:1201454/149‑1 (MQ=255)                                                                                          ttttCCGCGAAATGTCCAGCGATCCCGGGCGGCGCACCACCCGGTCGGAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAg                                                       <  2:187474/148‑1 (MQ=255)                                                                                                gcgcAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCgg                                              >  1:762383/1‑151 (MQ=255)                                                                                                   cAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa                                            >  2:103958/1‑150 (MQ=255)                                                                                                   cAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa                                            >  1:1051103/1‑150 (MQ=255)                                                                                                    aaTGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAc                                          >  2:433997/1‑151 (MQ=255)                                                                                                          caccGACCCCGTGCGGCGCATCACCCAGTCGAAAATAACCCCCCCCAGCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAgg                                      <  1:366607/146‑1 (MQ=255)                                                                                                            gCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                                  >  2:1016881/1‑151 (MQ=255)                                                                                                             cgcccccGGGCGGCGCACCGCCCAGGCGAAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg                                       <  1:193070/141‑1 (MQ=255)                                                                                                                 ccggggcggcgCACCAGCCAGACGAAAAAACCCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCa                                        <  1:159429/135‑1 (MQ=255)                                                                                                                   cGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaaga                             >  1:896231/1‑149 (MQ=255)                                                                                                                     tGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                          >  1:36800/1‑150 (MQ=255)                                                                                                                        ggcgCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATGAACCGGGAACCAGGTAATAAAGaagtaagt                      >  2:720018/1‑151 (MQ=255)                                                                                                                          cgcATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                     >  1:743110/1‑150 (MQ=255)                                                                                                                           gcATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                     >  1:479854/1‑149 (MQ=255)                                                                                                                             aTCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAggg                 >  1:1037493/1‑151 (MQ=255)                                                                                                                                        gaaaaTACCCCCCCCAAGCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACa       <  2:762413/148‑1 (MQ=255)                                                                                                                                         aaaaaaaCCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt     <  2:814719/150‑1 (MQ=255)                                                                                                                                          aaaTAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt     >  2:983704/1‑150 (MQ=255)                                                                                                                                            ataacCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACt  <  2:1027831/151‑1 (MQ=255)                                                                                                                                                 cccccccaCCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACt  <  2:432252/146‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               CGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACT  >  NC_000913/2921030‑2921319

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