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

AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCG  >  NC_000913/2921028‑2921324                                                                                                                                                        |                                                                                                                                                    aGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATcccac                                                                                                                                                    <  2:288534/151‑1 (MQ=255)  gcgcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATcccac                                                                                                                                                    <  2:714980/150‑1 (MQ=255)      cATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGc                                                                                                                                               >  1:300842/1‑151 (MQ=255)                atgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACAccgccg                                                                                                                                      >  1:1007639/1‑150 (MQ=255)                atgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGa                                                                                                                                     <  2:409860/151‑1 (MQ=255)                   atcatcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATa                                                                                                                                   >  2:133172/1‑150 (MQ=255)                                  tggACGTAGTTCCAGACCCCCATCCCCACCGACAGCAACATTCCCATTACGTTGGGTGTTTTACGCGAAAGGTTCAGCGA‑TCCCGCGCGGCGCATAAGCCAGGCGGAAAT‑AACCCCACCCACCACACCGCCGATAaaaccgcaaaccgc                                                                                                                     <  2:150291/147‑1 (MQ=255)                                    gACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcgg                                                                                                                  >  2:182793/1‑150 (MQ=255)                                      cGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGGCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaac                                                                                                               >  1:804693/1‑151 (MQ=255)                                        tAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgg                                                                                                             >  1:644919/1‑151 (MQ=255)                                           tggCAGAACCCCACCCCCATCGACAGAAACGTGCCCATTCCGACCGTTGTTTCACGCAAAATTTCCAGCGA‑CCCCGGGCGGGGCATAACCCAGCGGAAAAT‑AACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacggaa                                                                                                           <  1:630478/148‑1 (MQ=255)                                           tggCAAAACCCCACCCCCACCGACACCAACATCCCCATTACGATCGGTTTTTTCCGCGAAATTTCCAGCGA‑CCCCGGGCGGCACACAAGCCGGTCGAAAAT‑AACCCCCCCCAGCACCCCGCCGATAAAACCGCAAAccgccg                                                                                                                   <  1:770725/140‑1 (MQ=255)                                            tGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAg                                                                                                          >  1:589314/1‑150 (MQ=255)                                                gAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCcac                                                                                                      >  2:969988/1‑150 (MQ=255)                                                 aaCCCCACCCCCACCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGAAAAT‑AACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAg                                                                                                          <  2:404782/145‑1 (MQ=255)                                                   ccccatccccatcGACAGCAACATGCCCTTTCCGATCGGTGTTTTACGCGAAATTTTCAGCGA‑CCCGGTGCGGCGCATCAGCCAGTCGAAAAT‑AACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAg                                                                                                  <  2:1020297/149‑1 (MQ=255)                                                                gACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGa                                                                                      >  1:737895/1‑150 (MQ=255)                                                                  cAGCAACATGCCATTTACGCTCGGGTTTTTAGGCGCATTTTGCAGCGT‑CCCGGTGCGGGCCATCGCCCCGT‑GGAAATAAACCCCCCCCAGCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAg                                                                                       <  2:246965/147‑1 (MQ=255)                                                                       aCATGCCCATAACGACCGGTGTTTGACGCGCAATGTTCGGCGA‑TCCCGTGGGGCGCACCGCCCAGTCGAAAAT‑ATCCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGa                                                                                 <  2:406459/148‑1 (MQ=255)                                                                         aTGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATcc                                                                            >  2:31198/1‑151 (MQ=255)                                                                         aTGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATcc                                                                            >  1:500526/1‑151 (MQ=255)                                                                            cccttaaCGACCGGTTTTTTCCGCGCAATGTACAGCGA‑CCCGGTGCGGCGCACCACCCAGCGGGAAAA‑AATCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                                <  1:925392/138‑1 (MQ=255)                                                                            cccATTCCGATGGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                         <  2:721400/151‑1 (MQ=255)                                                                            cccATTACGTCGGGTGTTTTACGCGCAAT‑TTCCGCCACCCCCGTGCGGCGCATCGGCCAGTGGGAAAA‑AACCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATc                                                                             <  1:746675/147‑1 (MQ=255)                                                                             ccATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGTCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                        >  1:773535/1‑151 (MQ=255)                                                                              cATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgcg                                                                        >  1:640426/1‑150 (MQ=255)                                                                              cATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGt                                                                       >  1:240424/1‑151 (MQ=255)                                                                                 tACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGcc                                                                    >  1:214706/1‑151 (MQ=255)                                                                                    gacgggTGTTTACCGCGAATGGTCCAGCGA‑CCCCGGGCGGCGCACCAGCCAGTGGAAAAA‑ACCCCCCCCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGc                                                                     <  1:508699/143‑1 (MQ=255)                                                                                         gtgtTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                              >  2:660040/1‑149 (MQ=255)                                                                                         gtgtTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                              >  2:120316/1‑149 (MQ=255)                                                                                            ttttACGCGAAATGTCCAGCGA‑CCCCGGGGGGCCAATCAGCCAGTGGGAAAT‑AATCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGa                                                           <  1:354440/149‑1 (MQ=255)                                                                                             tttACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGAt                                                          >  2:444779/1‑149 (MQ=255)                                                                                               tACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAc                                                      >  1:335262/1‑151 (MQ=255)                                                                                                  gcgaAATGTTCAGCGA‑CCCCGGGCGGCCCACCAGCCGGCCGAAAAA‑AACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCgg                                                   <  1:552289/151‑1 (MQ=255)                                                                                                    gcAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCggg                                                  >  1:20105/1‑150 (MQ=255)                                                                                                      aaTGTTCAGCGA‑CCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGaa                                                <  2:117829/150‑1 (MQ=255)                                                                                                      aaTGTCCAGCGA‑TCCCGTGCGGCGCACCAGCCAGTCGAAAAT‑AACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCggg                                                  <  2:66400/148‑1 (MQ=255)                                                                                                        tGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAcc                                              >  1:738662/1‑150 (MQ=255)                                                                                                          ttCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAgg                                           >  2:144242/1‑151 (MQ=255)                                                                                                          tccAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAg                                            <  2:976645/148‑1 (MQ=255)                                                                                                                gA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaata                                      >  1:287897/1‑150 (MQ=255)                                                                                                                    cccGTGCGGCGCATCAGCCAGTGGAAAAT‑AACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGt                                          <  2:768581/143‑1 (MQ=255)                                                                                                                    cccGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaataa                                     >  2:625057/1‑148 (MQ=255)                                                                                                                    cccGGGCGGCGCACCACCCAGTCGAAAAT‑AACCCCCCCCGGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAgg                                           <  1:615546/142‑1 (MQ=255)                                                                                                                         gcggcgCATCAGCCAGTCGGAAAT‑CATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGACGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCCGAAGCAGATAAAACGGGAACGAGGAAATAAACaagta                              >  1:62569/1‑150 (MQ=255)                                                                                                                           ggcgCTTCACCCAGTCGAAAAT‑AACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaa                             <  1:603577/149‑1 (MQ=255)                                                                                                                           ggcgCATCAGCCAGTCGAAAAT‑AATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                           <  2:683307/151‑1 (MQ=255)                                                                                                                                accAGCCAGTGGAAAAA‑AACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAggg                      <  2:840618/149‑1 (MQ=255)                                                                                                                                     ccAGTCGAAAAT‑AACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCa                         <  2:553723/143‑1 (MQ=255)                                                                                                                                        gTCGAAAAT‑AACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGAt               <  2:667876/150‑1 (MQ=255)                                                                                                                                           ggAAAT‑AACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAg           <  2:495252/151‑1 (MQ=255)                                                                                                                                                  aaccccacccaCCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACt       <  1:167801/146‑1 (MQ=255)                                                                                                                                                     cccccccaCCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCg  <  1:900452/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                    AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGA‑TCCCGTGCGGCGCATCAGCCAGTCGGAAAT‑AATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCG  >  NC_000913/2921028‑2921324

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