Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,853	0	A	C	49.0%	67.9 / 15.2	49	S725A (TCA→GCA)	arcB	aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (21/4);  new base C (0/24);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 3.76e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTC  >  NC_000913/3350708‑3350993                                                                                                                                                  |                                                                                                                                             gCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccg                                                                                                                                         >  1:247410/1‑151 (MQ=255)        acGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCAcc                                                                                                                                  >  1:537551/1‑151 (MQ=255)            tttCAGCCCTTCTCCGTCGTGCCCCCCCTCTTCTTTCACCTCTTCAACCCATCCCCCGACTTTTTCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCCCTGACCCCGCCGCACCt                                                                                                                                 <  2:831396/148‑1 (MQ=255)                     ttCTCCGTCGTGACCCCCCTCTTCTTTCACCCCTTAAATCCATCCCCCGACGTTATCTCCCCAGGCGGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTAt                                                                                                                       <  1:871716/149‑1 (MQ=255)                       ctgcgtggtGCCCCCACTCTTCTTCCATCTCTTCAACCCTTTCCCCGACGTTATCTTCCCAGGCCGGAGGGTCAGGAGACTGATTTTGCTCCCCTAGCTTTTGGTGAGGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATg                                                                                                                      <  2:735307/145‑1 (MQ=255)                       cTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGtcc                                                                                                                   >  1:200254/1‑151 (MQ=255)                          cgtcgtGACCCCCCTCTTCTTTCATCTCTTAAATCCATTCCCCGACGTTATCTTCCCAGGCCGAAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTTTTGCAGATGGCGTAACCCCCCTGCCCCCGCCGCACCTTTAATTTTATGtc                                                                                                                    <  1:823175/147‑1 (MQ=255)                               tGACCCCACTCTTCTTCCATCTCTCCAACCCATCCCCCGACTTTATCTTCCCAGGCGGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCCCTGACCCCGCCGCACCTTTAATTTTATGtc                                                                                                                    <  1:180065/142‑1 (MQ=255)                                  cccccctctTCTTTCATCTCTCCAACCCATTCCCCTACGTTATCTCCCCAGGCCGAAAGGTCAGGATATTGAATTTGCCCCCCCAGGTTTTGCAGATGCCCTAACCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaacaa                                                                                                        <  2:531182/146‑1 (MQ=255)                                  cGCCACTCTTCTTTCATCTTTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                         >  1:978685/1‑150 (MQ=255)                                   gCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                                       >  2:958440/1‑151 (MQ=255)                                    cctctctTCTTTCATCTCTTCAACCCATCCCCCGACGTTATCTCCCCAGCCGGGAGGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                                       <  1:382373/147‑1 (MQ=255)                                             tttCACCTCTACAACCCATTCCCCGACGTTATCTCCCCAGGCCGAAGGGTCAGGAGACTGAATTTGCTGCCCCAGTTTTTGCAGATGCCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                     <  2:91195/143‑1 (MQ=255)                                              ttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                            >  1:136349/1‑151 (MQ=255)                                               tCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGTCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATc                                                                                           >  1:1047969/1‑151 (MQ=255)                                                  tctcTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                                        >  2:41907/1‑151 (MQ=255)                                                      ttCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGcc                                                                                    >  1:174200/1‑151 (MQ=255)                                                          aTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                 >  1:485835/1‑150 (MQ=255)                                                             cactcccccactttATCTCCCGGGGCCGGAGGGTCGGGGGACTAAATTTGCTGCCCTTGTTTTTGCAGTTGGCGTACCCCCTCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgcg                                                                                      <  2:22202/130‑1 (MQ=255)                                                               tcccccgcctttATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                            <  1:882768/140‑1 (MQ=255)                                                                  cccGACTTTCTCTCCCCGGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCGGTTGTTGCAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTc                                                                           <  2:853138/147‑1 (MQ=255)                                                                   ccGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGTTTTTGCAGTTGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                             <  2:865363/145‑1 (MQ=255)                                                                       ctttttctcCCCAGGCCGGAGGGTCGGGAGACTGAATTGGCTGACCCAGTTGTGGCGGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                                    <  2:963228/145‑1 (MQ=255)                                                                         ttATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTTTTGAAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                  <  2:932896/150‑1 (MQ=255)                                                                         ttATCTCCCCAGGCCGGAAGTTCAGGAGACTGATTTTGCTGACCCGTCTGTTGCAGATGGCGTAACCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                                    <  1:547876/148‑1 (MQ=255)                                                                               tCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                            >  2:1075750/1‑150 (MQ=255)                                                                                  cAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                        >  2:695944/1‑151 (MQ=255)                                                                                    ggCCGGAAGGTCAGGAGACTGATTTTCCTGCCCCAGCTTTTGCAGAGGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa                                                       <  1:661585/150‑1 (MQ=255)                                                                                    gccgggAAGGCCAGGACACAAATTTTGCCCACCCGGCTTTTCCAGAGGCCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAg                                                      <  2:982791/147‑1 (MQ=255)                                                                                     gCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGc                                                     >  1:387078/1‑151 (MQ=255)                                                                                       cGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGccc                                                   <  1:929890/151‑1 (MQ=255)                                                                                             gTCGGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatc                                             <  1:641870/151‑1 (MQ=255)                                                                                                   agaCTGATTTTGCTGCCCCAGCTGTTGAAGATGGCGTACCCCCACGCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                                       <  2:1021249/151‑1 (MQ=255)                                                                                                    gaCTGATTTGGCTGACCCAGCTGTTGCAGAGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                                       <  1:861339/150‑1 (MQ=255)                                                                                                     aCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTtctc                                     >  1:1021234/1‑151 (MQ=255)                                                                                                           tttccgcccccATCTGTTGAGGTTGGCTAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAaca                                <  2:668138/142‑1 (MQ=255)                                                                                                            ttGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                              >  2:236218/1‑151 (MQ=255)                                                                                                                tGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                           >  1:762343/1‑150 (MQ=255)                                                                                                                  aCCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAc                        >  1:825112/1‑151 (MQ=255)                                                                                                                    ccAGCTGTTGCAGTTGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAccc                      <  1:898216/151‑1 (MQ=255)                                                                                                                    ccAGCTGTTGAAGATGGCCCACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAAccc                      <  2:90332/151‑1 (MQ=255)                                                                                                                        cTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCg                  >  1:938047/1‑151 (MQ=255)                                                                                                                            tgcagccgttgTACCCCCACAGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc                   <  2:221580/136‑1 (MQ=255)                                                                                                                              cAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc            >  1:697946/1‑151 (MQ=255)                                                                                                                               agaggccttacccccccTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa           <  2:1057292/140‑1 (MQ=255)                                                                                                                                 ctGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGGCGTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTAAAACACCGCTAACCCGTCGGTGATCAGc         >  2:845337/2‑151 (MQ=255)                                                                                                                                  tGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa           <  2:87277/148‑1 (MQ=255)                                                                                                                                        acccccACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg     <  2:926972/146‑1 (MQ=255)                                                                                                                                        aaCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTc  >  2:113578/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                             GCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTC  >  NC_000913/3350708‑3350993

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