Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	46.8%	60.6 / 15.4	47	*559E (TAA→GAA)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (21/4);  new base C (0/22);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.47e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATC  >  NC_000913/501416‑501709                                                                                                                                                       |                                                                                                                                                  aCTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑Acc                                                                                                                                                 <  1:36980/151‑1 (MQ=255)    gggATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCg                                                                                                                                              >  1:902387/1‑151 (MQ=255)       aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTca                                                                                                                                           >  2:684532/1‑151 (MQ=255)       aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTc                                                                                                                                            <  2:693338/150‑1 (MQ=255)        tCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTca                                                                                                                                           >  1:742485/1‑150 (MQ=255)         cGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTcac                                                                                                                                          >  2:329248/1‑150 (MQ=255)          gTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTcacc                                                                                                                                         >  1:502540/1‑150 (MQ=255)                ttATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCAc                                                                                                                                  <  2:546860/151‑1 (MQ=255)                           tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑Accc                                                                                                                                                <  1:519724/126‑1 (MQ=255)                           tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑Accc                                                                                                                                                >  2:519951/1‑126 (MQ=255)                           tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTgg                                                                                                                        >  2:497549/1‑150 (MQ=255)                           tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTgg                                                                                                                        >  1:44904/1‑150 (MQ=255)                                ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                                                                                   >  1:606334/1‑150 (MQ=255)                                ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                                                                                   >  1:289126/1‑150 (MQ=255)                                     gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                             >  1:391959/1‑151 (MQ=255)                                     gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                             >  2:575475/1‑151 (MQ=255)                                     gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                             >  2:241770/1‑151 (MQ=255)                                                        ttttggtttttGGATCTCGCTGCTGAATGGGCCGTGGGCGCGTTTTTTCAGGCCTAAATTTTTTCCATTTTTTTCATCCGGCAACCCCCTTTT‑‑CCCCTGCACCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                                                                                          <  2:169221/143‑1 (MQ=255)                                                            ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                      >  1:713094/1‑151 (MQ=255)                                                                             gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                      >  1:73642/1‑150 (MQ=255)                                                                              cttGAGGCGACGCGGGCGCGTCTTTCCATGCTTCAATTTTTTTCATTTTTTCCA‑CCGGCCACCCCCGTTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                     <  1:151535/149‑1 (MQ=255)                                                                                  ttGCGCCGCTGGCGCGTCTTACCAGGCCTACATTTTTTTCATATTTTCCACCCGGCACCCCCCTTTT‑‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                 <  2:91441/149‑1 (MQ=255)                                                                                  aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                  >  2:293778/1‑149 (MQ=255)                                                                                     cgccgttGGCGCGTCTTACCATCCCTACATATTTTTCATATTTTACACCCGCCACCCCCCGTTT‑‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                             <  1:853748/145‑1 (MQ=255)                                                                                      caCGCGGGCGCGTTTTATGACACCTACATTTTTTTCATATTTTACACCCGGAAACCCCCTTTA‑‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                            <  1:60681/150‑1 (MQ=255)                                                                                           cggCGCTTTTTATCATGCCTACTTATTTTTCATATTTTACACCCGCCAACCACCTTTT‑‑CCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                       <  1:264731/150‑1 (MQ=255)                                                                                                 tttttATCAGCCCTACATTTTTTTCATTTTTTCCACCCGCCACCCCCCTTTC‑‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                 <  1:497338/148‑1 (MQ=255)                                                                                                    ttATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGCTTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                              >  2:954545/1‑151 (MQ=255)                                                                                                     tcccttcccTACATTTTTTTCATTTTTCCCACCCGCCACCCCCCGTTT‑‑CCCCCCCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                             <  2:998140/144‑1 (MQ=255)                                                                                                     taccATGCCTACTTATTTTTCATTTTTTACTGCCGGCAACCGCCGTTT‑‑CCCCGGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                             <  1:347247/148‑1 (MQ=255)                                                                                                     tATCTTCCCTACATTTTTTTCATTTTTTACCCCCGCCACCCACCGTTT‑‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                             <  1:1049286/151‑1 (MQ=255)                                                                                                     tATCATGCCTACATATTTTTCATATTTTACATCCGGTAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATGTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                             >  1:1035787/1‑151 (MQ=255)                                                                                                        cATGCCTACTTTTTTTCCATTTTTTCCCCCCGCCACCCCCCTTTT‑‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                          <  1:116432/151‑1 (MQ=255)                                                                                                        cATGCCTACATTTTTTTAATTTTTTACACCCGGCACCCACCGTTT‑‑CCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                          <  1:879498/151‑1 (MQ=255)                                                                                                          tGCCTCAATTTTTTTCATTTTTTCCACCCGGCACCCCCCGTTT‑‑CCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                        <  2:567499/151‑1 (MQ=255)                                                                                                          tGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                        >  2:464847/1‑151 (MQ=255)                                                                                                            ccTACTTATTTTTCATTTTTTACACCCGGCACCCCCCTTTT‑‑CCCCCGCCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                      <  1:584491/151‑1 (MQ=255)                                                                                                            ccTACTTATTTTACTTATTTTACATCCGGCACCCCCCGTTT‑‑CCCCCTCCCCCACCTCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                      <  1:716383/151‑1 (MQ=255)                                                                                                             cTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGGGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                     >  2:548946/1‑151 (MQ=255)                                                                                                              tACATATTTTTATTATTTTACACCGGGCACCCCCCGTTT‑‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                    <  2:459803/151‑1 (MQ=255)                                                                                                                        tCATTTTTTCCACCCGCCACCCACCGTTT‑‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                          <  1:306579/151‑1 (MQ=255)                                                                                                                                           cccccctttt‑‑CCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaa       <  2:906464/144‑1 (MQ=255)                                                                                                                                               cctttt‑‑CCCCCGCCCCCTCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat   <  2:502756/148‑1 (MQ=255)                                                                                                                                               ccGTTT‑‑CCCCCGTCCCCACCACACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat   <  1:386283/151‑1 (MQ=255)                                                                                                                                               ccGTTC‑‑CCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat   <  1:399207/151‑1 (MQ=255)                                                                                                                                                cGTTT‑‑CCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc  <  1:426277/151‑1 (MQ=255)                                                                                                                                                cGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCCCGCCCATCACTACCTGATTCGCACCACGTTCGGCGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat   >  2:559085/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                  ACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATC  >  NC_000913/501416‑501709

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