Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	298,960	0	A	C	48.4%	24.5 / 13.3	31	Y659D (TAT→GAT)	paoC	PaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (9/7);  new base C (0/15);  total (9/22)
Fisher's exact test for biased strand distribution p-value = 8.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACA  >  NC_000913/298816‑299104                                                                                                                                                  |                                                                                                                                                 acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCg                                                                                                                                             <  1:277200/150‑1 (MQ=255) acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCg                                                                                                                                             <  2:500582/150‑1 (MQ=255)    ccGCGTTGGCGATAACCGCGCTCACGCCGCACAGGCCCAGCTCCCCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGG‑CACCTCATACCCCGCCat                                                                                                                                         <  2:177258/151‑1 (MQ=255)          gggCGATAGCCGCGCTCACGCCGCACAGCCCCAGCTCACCGACCCCTTGGGCCTTCATGGGGGAGGATATGGGGTCGGTATCACCCAGAAAATCCCCCTCCTTTTTTGGATTGCCCGATGGACCCGC‑CACCCCATACCCCGCCATATc                                                                                                                                      <  2:556972/147‑1 (MQ=255)              gATCGCCGCGCTCACGCCGCACAGGCCCAGCTCCCCGACACCTTGGCCCTTCACGGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAATCCCCCCCTTTTTGGGAAAACCCGCAGGAACCGG‑CACCTCATCCCCCGCCATATCGTGAt                                                                                                                                 <  2:3523/149‑1 (MQ=255)               agaGCCGCGCTCACGCCGCACAGGCCCAGCTCCCCGACCCCTTTGCCCTTCATCGGGGAGGATAGGGGGTCGGTATCTCCCAGAAAAATCACCCCCTGTTTTGGGATACCCGCATGAACCGG‑CCCCTCATACCCCGCCATATCGTGATTaa                                                                                                                              <  1:447228/149‑1 (MQ=255)                tggCCGCGCTCACGCGGAACAGGCCCAGCTCCCCGACCCCTTGGGCCTTCATGGGGGAGGATATGGGGTCGGTATCATCCAGGAAAACCCCCCCCTTTTTGGGGATACCCGCATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTaa                                                                                                                              <  2:493056/148‑1 (MQ=255)                   ccGCGCTCCCGCCGCACGGCCCCAGCTCACCGACACCTTTGGCCTTCAGGGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAATCACCTCCGGTTTTGGAATATCCGCATGAACCGG‑CCCCTAATACCCCGCCATATCGTGATTAACGa                                                                                                                           <  2:166216/150‑1 (MQ=255)                      cgcTCCCGCCGCAAGGGCCCAGCTCACCGACCCCTTGGGCTTTCACGGGGGGGGATAGGGGGTGGGTATCACCCAGGAAAACCCCCTCCTTTTTTGGGATATCGGAATGAACGGG‑CACCTCATCCCCCGCCATATCGTGATTAACGAAGTa                                                                                                                       <  2:515497/151‑1 (MQ=255)                              cctcagaGGCCAATCTCCCGGACCCTTTTGCCCTTCTTGGGGAAATATATGGGGTCGTTCTCACCCAGAAAAATCCCCCCCTTTTTGGGTATATCCCAAAAACCCGG‑CCCCGCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                  <  1:436971/142‑1 (MQ=255)                                    aGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt                                                                                                         >  2:247308/1‑151 (MQ=255)                                       cccAGCTCCCCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCa                                                                                                      <  1:158086/151‑1 (MQ=255)                                            cTCCCCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAATCCCCCCCTGTTTTGGAATACCCGCATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc                                                                                                  <  1:276847/150‑1 (MQ=255)                                              ccccacccccTTGGGCCTTCAGGGGGGAGGATATGGGTCCGGTATCATCCAGAAAAATCCCCTCTTTTTTTGGATTATCCGAATCACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg                                                                                                   <  2:290116/139‑1 (MQ=255)                                                 cGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctc                                                                                            >  1:308631/1‑151 (MQ=255)                                                     acCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAATGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc                                                                                        >  2:391855/1‑151 (MQ=255)                                                        tttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                     >  2:228814/1‑151 (MQ=255)                                                         ttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCAGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                     >  2:562677/1‑150 (MQ=255)                                                                                ggggTGGGTATCA‑CCAAGAAAAACCCCCCCCTTTTTGGGAATACCGCAAGCACCGGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAg                                                               <  2:331022/149‑1 (MQ=255)                                                                                ggggTCGGTATCACCCAGGAAAACCCCCCCCTGTTTGGGAATTCCCGCATGACCCGC‑CCCCCCTTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATa                                                                <  1:155134/148‑1 (MQ=255)                                                                                 gggTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                             >  1:54807/1‑150 (MQ=255)                                                                                      ggTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg                                                       >  2:384846/1‑151 (MQ=255)                                                                                             cccAGAAAATTCCCCCCCTGTTTGGGGATCCCCCCATGAACCGC‑CCCCTCCTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                  <  1:270059/148‑1 (MQ=255)                                                                                                   aaacccccctcctTTTTGGGGATATCCGCATGACCCGC‑CACCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa                                                <  1:523370/138‑1 (MQ=255)                                                                                                          cctcctGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                   >  2:251079/1‑151 (MQ=255)                                                                                                                 ttttGGGATACCCGCATAACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                             <  1:533234/150‑1 (MQ=255)                                                                                                                     gggATACCCGCAGGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                          <  1:444786/149‑1 (MQ=255)                                                                                                                         tatCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                      >  1:118610/1‑149 (MQ=255)                                                                                                                          acccGCATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAg                    <  1:500504/148‑1 (MQ=255)                                                                                                                                        gccACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc               <  1:678326/140‑1 (MQ=255)                                                                                                                                            ccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcaca  <  1:139450/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 ACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACA  >  NC_000913/298816‑299104

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