Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,844	0	A	C	52.3%	54.1 / 70.2	44	L728V (TTA→GTA)	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 (18/3);  new base C (0/23);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.29e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-08

TTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCA  >  NC_000913/3350695‑3350984                                                                                                                                                      |                                                                                                                                             tttAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTaa                                                                                                                                              <  1:517509/150‑1 (MQ=255)   tAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAAccc                                                                                                                                           >  1:534316/1‑151 (MQ=255)        gCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTg                                                                                                                                      >  1:307846/1‑151 (MQ=255)                cacccacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCccgcc                                                                                                                              <  1:450119/151‑1 (MQ=255)                    cacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCAc                                                                                                                          >  1:581303/1‑151 (MQ=255)                           tCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAt                                                                                                                   >  1:472181/1‑151 (MQ=255)                             aGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAttt                                                                                                                 >  2:388852/1‑151 (MQ=255)                              gCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                >  2:154321/1‑151 (MQ=255)                              gCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAAtttt                                                                                                                >  2:540392/1‑151 (MQ=255)                                  ttctcctccgtGCCCCCCCTCTTTTTTCATCTCTCCAACCCATCCCCCGACTTTCTCTTCCCAGGCGGGAAGGCCAGGAGACTAATTTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCCCTGACCCCGCCGCACCTTTaa                                                                                                                    <  1:400311/135‑1 (MQ=255)                                  ttCTGCGTCGTGCCCCCTCTCTTCTTTCATCTCTTCAACCCATCCCCCGCCTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTAAATTTGCTGCCCCAGTTTTTGCAGAGGGCGAACCCCCACTGCCCCCGCCGCACCTTTAATCTTAt                                                                                                              <  1:310630/149‑1 (MQ=255)                                     tccgtcgtGACCCCACTCTTCTTTCATCCCTTCAATCCATTCCCCGCCGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCGGCCCCAGCTGTTGCAGAGGGCGTAACCCCACTGCCCCCGCCGCACCTTTAATTTTATg                                                                                                             <  1:379401/145‑1 (MQ=255)                                             gACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTCGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaa                                                                                                  >  1:381242/1‑150 (MQ=255)                                              accccACTCTTCTTTCATCTCTTCAATCCATCCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaa                                                                                                  <  2:221491/146‑1 (MQ=255)                                               cGCCCCTCTTTTTTCTTCTCTTAATTCCATCCCCCGACGTTATCTCCCCAGGCGGGAAGGTCGGGAGACTAATTTTGTTGCCCCAGTTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                 <  2:239762/149‑1 (MQ=255)                                                 cccctctTCTTACATCTCTCAAACCCTTCCCCCGACTTTATCTTCCCGGGCGGGAGGGTCGGGAGACTTATTTTGCTCCCCAATTTGTTGAAGAGGGCCAACCCCCACGGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAAt                                                                                              <  1:444825/147‑1 (MQ=255)                                                 cccctctCCTTTAATCTTTCCAATCCATCCCCCGCGTTTTTCTCCCCAGGCCGGAAGGTCAGGAGTCTAAATTTGTTTACCCAGCTGTTGCAGAGGGCGAAACCCCCCTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTca                                                                                                   <  2:667898/142‑1 (MQ=255)                                                 ccACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATg                                                                                             >  1:122688/1‑151 (MQ=255)                                                       ttcttTCACCTCTTCAACCCATCCACCGCCTTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGATTTTCCTGACCCAGTTTTTGCAGTTGGCGTACCCCCACTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCt                                                                                          <  1:338872/148‑1 (MQ=255)                                                          tttCATCTCTTCAACCCTTCCCCCGACGTTATCTCCCCAGGCGGGAGGGTCAGGAGACTGAATTTGCTGACCCAGCTTTTGCAGATGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttt                                                                                        <  2:582304/147‑1 (MQ=255)                                                          tttCATCTCTTAAACCCATCCCCCGCCGTTATCTCCCCAGGCGGGAGGGTCAGGAGACTAATTTTGCTGCCCCAGTTTTTGAAGATGGCGTACCCCCCCTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt                                                                                       <  1:265199/148‑1 (MQ=255)                                                           tccATCTCTTAAACCCATCCCCCGACTTTATCTCCCCAGGCCGGAGGGTCAGGAGACTGAATTTGCTGCCCCAGTTTTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                     <  1:412143/147‑1 (MQ=255)                                                                  ctccAACCCATCCCCCGATGTTATCTCCCCGGGCCGGAAGGTCAGGAGACTAAATTTGCTGCCCCATTTTTTGAAGTTGGCTAACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTg                                                                               <  1:160172/145‑1 (MQ=255)                                                                       aTCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                        >  1:567183/1‑150 (MQ=255)                                                                            tcccccGACTTTATCTCCCCAGGCCGGAGGGTCGGGAGACTGAATTTGCTGACCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                   <  2:618564/146‑1 (MQ=255)                                                                              ccccGACGTTACCTTCCCAGGCCGGAGGTTCGGGAAGCTGATTTGGCTACCCCAGCTGTTGCAGAGGCCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGa                                                                <  2:534319/149‑1 (MQ=255)                                                                              cACCGACTTTATCTCCCCAGGCGGGAAGGTCGGGATACTGAATTTGCTGCCCCAGTTGTTGCAGATGGCGTACCCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                    <  2:550690/147‑1 (MQ=255)                                                                                ccGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACt                                                              >  2:470118/1‑151 (MQ=255)                                                                                      ttATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                         <  2:639478/150‑1 (MQ=255)                                                                                       tATCTCCCCAGGCAGGCGGGTCAGGACAATGACTTCGCTTCCCAACATTTTGTGGAGGGCCTACCTCCCCAGTCCCCGCCGCACCTTTAAGTTTATGTCCTTCCTCAACAATGCCTTCTTTATCCTGCGCCGCCAGATTCGACTcaatc                                                         <  1:400201/149‑5 (MQ=255)                                                                                         tCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGc                                                     >  1:282374/1‑151 (MQ=255)                                                                                           ttCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                   >  2:421093/1‑151 (MQ=255)                                                                                             cccAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTTTTGCAGATGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCCACTCCAGCACGCTg                                                   <  2:235974/149‑1 (MQ=255)                                                                                                 ggCCGGAAGGTCAGGAGACTGATTTTGCGGCCCCAGTTGTTGCAGATGGCTTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATa                                              <  2:545337/150‑1 (MQ=255)                                                                                                          gTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatc                                    >  2:39752/1‑151 (MQ=255)                                                                                                              ggAGACTAATTTTCCTGCCCCAGCTGTTGCAGAGGGCATACCCCCCCTCCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatca                                   <  1:665807/148‑1 (MQ=255)                                                                                                              ggAGACGAAATTTGCTGACCCAGCTGTTGCAGAGGGCGAACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatca                                   <  2:291450/148‑1 (MQ=255)                                                                                                                          tGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                    >  1:418745/1‑151 (MQ=255)                                                                                                                            cTGCCCCAGCTGTTGCAGATGGCGTACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGc                   <  2:300507/150‑1 (MQ=255)                                                                                                                                cccAGTTGTTGCAGAGGGCGTACCCCCACTGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTaa                <  2:340466/149‑1 (MQ=255)                                                                                                                                    gCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTc          >  2:181728/1‑151 (MQ=255)                                                                                                                                          gCAGATGGCGTACCCCCCCGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGAt    <  2:392554/151‑1 (MQ=255)                                                                                                                                          gCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGAt    >  1:463810/1‑151 (MQ=255)                                                                                                                                           cAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATc   >  2:71678/1‑151 (MQ=255)                                                                                                                                             gagggCGAACCCCCACGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa  <  2:69604/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             TTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAACCCCACTGACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCA  >  NC_000913/3350695‑3350984

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