Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	44.0%	77.4 / 10.3	50	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (27/1);  new base C (0/22);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 2.59e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993776‑994063                                                                                                                                                  |                                                                                                                                                   cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCca                                                                                                                                           >  2:352752/1‑151 (MQ=255)             aaCATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                  >  2:21662/1‑148 (MQ=255)                     tCGCGGTTCTT‑CTGAATTTT‑CACCCGGGGGTG‑GGTGCCCCGTTAACACCCCCCCTGGGTTGGGGGTTCCCAGCCCGCCCAGAGGCCCCAGCAGGTTTTTTTTCCCC‑CCCCCGCGGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAata                                                                                                                       <  1:1234114/151‑1 (MQ=255)                              ttCCCGAA‑TTTCCACCCGGGGGTGGGGGC‑CCCGCAAACACCCCCCCGGGGGTGGGCTTTCCCAGCCCTCCCAGCAGGCGCAGCAGGGTATTTTTCCCA‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                               <  1:106040/151‑1 (MQ=255)                                ccTGAA‑TTT‑CAGCCAGCGGTGTGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGCAAATCCAg                                                                                                            >  2:178335/1‑151 (MQ=255)                                     aaatcccagcgcgcGGTGGGGTG‑CCCCAAAACACCTCGCCTGCGGTGGGCGTTCCCAGACCGGCCGGCGGGGGCAGCAGGGTAC‑TTTGCCACCCCCCGAGCCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                         <  2:558154/139‑1 (MQ=255)                                     a‑TTT‑CAGCCAGCGGTGTGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       >  2:1189727/1‑151 (MQ=255)                                     a‑TTT‑CAGCCAGCGGTGTGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       >  1:1006687/1‑151 (MQ=255)                                                    gtgtGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                          >  1:314377/1‑151 (MQ=255)                                                    gtgtGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                          >  1:471227/1‑151 (MQ=255)                                                    gtgtGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                          >  1:631040/1‑151 (MQ=255)                                                                 cTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCATCAGGGCACCTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGATCAGGACGATATTTTCCGCGTAATg                                                                              >  2:365743/1‑151 (MQ=255)                                                                 cTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                              >  1:689333/1‑151 (MQ=255)                                                                   aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            >  2:559909/1‑151 (MQ=255)                                                                     cacaTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          >  1:844527/1‑151 (MQ=255)                                                                          cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                      >  2:572413/1‑150 (MQ=255)                                                                            ccgggggTTGGCGTTCCCAGCCCCCCCAGCAGGCGCGGCGGGTTATTTTCCCCC‑CCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                   <  1:1034668/146‑1 (MQ=255)                                                                                 ggTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              >  2:744612/1‑151 (MQ=255)                                                                                         ttcccaccccccccAGCGGGCGCGGCAGGTTATTTTGCCCC‑CCACCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                      <  1:214300/139‑1 (MQ=255)                                                                                         tccccaccccTGCCAGCGGGCCCGGGGGGGTTCTTTCCCCA‑CCCCCGCGGCGGCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                      <  1:883643/143‑1 (MQ=255)                                                                                            caactccagCCAGCAGCGGCAGCTGGGTCCTTGTCCCC‑CCACGGGCGCTGCCCACCGCCGCCACAAACTGTACTGCCTGAAGATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATTTTTGCTGATTGCAGTGAGCAACAAggg                                                   <  2:250082/143‑1 (MQ=255)                                                                                            caaccccgcccAGCGGGCGCAGGGGGGTTTTTTTGCCC‑CCCCCGTGTCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  2:774422/142‑1 (MQ=255)                                                                                              accccccccAGCGGGCGCAGCGGGTTCTTTTTCCCC‑CCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                 <  1:744928/144‑1 (MQ=255)                                                                                                    gCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                             >  1:1072968/1‑149 (MQ=255)                                                                                                          gggCGAAGCAGGATTTTTCCGCCG‑CCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      <  1:663196/149‑1 (MQ=255)                                                                                                            gcggAGCAGGGTTTTTTCGCCC‑CCCCCGCTGGGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      <  1:254058/148‑1 (MQ=255)                                                                                                            gcgcAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                   >  1:449523/1‑151 (MQ=255)                                                                                                             cgcgGCAGGTTACTTTCCCCC‑CCCCCGCGGGGGCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  <  2:841933/151‑1 (MQ=255)                                                                                                              ccagcggGGTCCTTTGGCCC‑CCCCCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 <  2:489599/150‑1 (MQ=255)                                                                                                                 gcagGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                               >  1:260276/1‑150 (MQ=255)                                                                                                                    gggTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                           >  2:299937/1‑151 (MQ=255)                                                                                                                      gTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                         >  2:1255556/1‑151 (MQ=255)                                                                                                                        ttttttGCCC‑CCCCGGCGGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  2:862413/149‑1 (MQ=255)                                                                                                                        ttttttCCCC‑CCCCCGCGGGGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  1:883567/149‑1 (MQ=255)                                                                                                                        tttctcccca‑cccCCGCCGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  2:468940/144‑1 (MQ=255)                                                                                                                         cTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                      >  2:188206/1‑151 (MQ=255)                                                                                                                           tttgcccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                     <  1:118842/140‑1 (MQ=255)                                                                                                                            ttGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa                    >  1:499752/1‑150 (MQ=255)                                                                                                                             ttcccccccccccggcgcCCCCCCACCGCCACAAACTGCCCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  1:954592/137‑1 (MQ=255)                                                                                                                             ttcccccccccGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:529554/142‑1 (MQ=255)                                                                                                                             tccccccccccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  1:643243/142‑1 (MQ=255)                                                                                                                             tccccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  1:175585/142‑1 (MQ=255)                                                                                                                             tccccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  1:21667/142‑1 (MQ=255)                                                                                                                             tcccaccccgcggtGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:1044262/138‑1 (MQ=255)                                                                                                                             tGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   >  2:1207516/1‑150 (MQ=255)                                                                                                                                    ccagggcgcgccacacaacCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACGGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc            <  2:489110/138‑1 (MQ=255)                                                                                                                                        gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATTCTGTGGcctc          >  2:567393/1‑149 (MQ=255)                                                                                                                                             ggCCCACCACCGCCACACACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaaa   >  1:1082056/1‑151 (MQ=255)                                                                                                                                              gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  1:881116/1‑151 (MQ=255)                                                                                                                                              gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  2:1208076/1‑151 (MQ=255)                                                                                                                                              gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  1:1160966/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTG‑CCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993776‑994063

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