Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,853	0	A	C	51.5%	36.5 / 48.9	33	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 (12/4);  new base C (0/17);  total (12/21)
Fisher's exact test for biased strand distribution p-value = 5.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-05

GCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350717‑3351003                                                                                                                                            |                                                                                                                                                       gCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCtt                                                                                                                                          >  2:224349/1‑151 (MQ=255)       aGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAAttt                                                                                                                                    >  1:3477101/1‑151 (MQ=255)            ttCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATg                                                                                                                                >  2:3675601/1‑150 (MQ=255)            ttCTACGTCGTGACCCCCCTCTTCTTTCATCTCTACAACCCATCCACCGCCTTTATCTCCCCAGGCCGGATGGTCGGGAGACTGAATTTGCTCCCCCATTTT‑TTGAG‑GATGGGCCACCCCCACTC‑CCCCCGCCGCACCTTTAAtttt                                                                                                                                   <  1:1682000/147‑1 (MQ=255)                    tgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                       >  1:1808724/2‑151 (MQ=255)                     gtGACGCCCCTTTTCTTTCATCTCTTCAATCCATTCCCCGACTTTATCTCCCCAGGCCGGAAGGTCAGGAGATTTATTTTGCT‑ACCCCACTGCTTGCA‑AGTGGGCTACCCCCGCCT‑CCCCCGCCGCACCTTTAATTTTATGtc                                                                                                                              <  2:791179/143‑1 (MQ=255)                      tgccccccctctTCTTTCATCTCTTCAACCCATTCCCCGACGTTATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTT‑TTGCA‑GAGGGCGTACCCCCACTG‑CCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                           <  1:839743/137‑1 (MQ=255)                         cccccctctTCTTTCACCTCTTCAACCCTTTCCCCGACTTTATCTCCCCAGCCCGGAAGGTCAGGAGACTGAATTTCCTTCCCCAGCTC‑TTGAA‑AAGGCCTAAACCCCACTC‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaac                                                                                                                    <  2:2071358/144‑1 (MQ=255)                           cccctctTTTTTCATCTTTCAAACCCATTCCCCGCCTTTATCTCCCCAGGCGGGAAGGTCGGAAGACTGATTTTGTTTTCCCGGCTG‑TAGGA‑GGGGGTAAACCCCAACTG‑CCGCCGCCGCACCTTTAATTTTATGTCCTTCCTcaa                                                                                                                     <  2:2117848/143‑1 (MQ=255)                                     ttCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTAt                                                                                                      >  1:3602211/1‑151 (MQ=255)                                      tCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATc                                                                                                     >  1:3482747/1‑151 (MQ=255)                                              tCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCg                                                                                             >  1:858640/1‑151 (MQ=255)                                                        ccccacggTTTTCTTCCAAGGCCGGAAGGTCGGGCGCCTGATTTGCCTGACCCAGTTG‑TTGAA‑GATGCCGTACCCCCACGG‑CCACCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAg                                                                                         <  1:469741/138‑1 (MQ=255)                                                         cccGACTTTATCTCCCCGGGCCGGAAGGTCAGGAGACTAATTTTGCTACCCCAGTTT‑TTGAA‑GATGGCGTACCCCCACTC‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    <  1:3002428/148‑1 (MQ=255)                                                             actttATCTCCCCAGGCCGGAAGGTCAGGAGACTGAATTGGCTGCCCCAGCTT‑TTGCA‑GATGGCGTACCCCCTCTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                                              <  1:3437701/148‑1 (MQ=255)                                                             aCGTTATCTCCCCAGGCGGGAAGGTCAGGAGACTAATTTTGCTCACCCATTTT‑TTGAA‑GATGGCGAACCCCCACTC‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                                              <  1:1057753/151‑1 (MQ=255)                                                              cGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTcc                                                                               >  2:2530267/1‑149 (MQ=255)                                                                   tCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGc                                                                        >  2:944291/1‑151 (MQ=255)                                                                    cTTCCCAGGGCGGAGGGTCAGGACACTAATTTTGCTGCCCCAGGTG‑TTGCA‑AATGGCTTACCCCCACGG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACg                                                                         <  2:1857499/149‑1 (MQ=255)                                                                     ttCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                                      >  2:2647551/1‑151 (MQ=255)                                                                       cccAGGCCGGAGGGTCGGGGGACGGAATTTTCTCCCCCAGTTT‑TTGAA‑GATGGCAAAACCCCACTG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGa                                                                     <  1:2227/150‑1 (MQ=255)                                                                       cccAGCCGGGAAGGTAGGGAGATTAATTTTGCGCCCCCTGCTG‑TTGCA‑GTTGGCGAACCCCCACTC‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCa                                                                              <  1:2062184/141‑1 (MQ=255)                                                                           ggCCGGAAGGTCAGGAGACTAAATTTGCTGACCCAGCTG‑TTGCA‑GAGGGCGTACCCCCACTG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACAt                                                                  <  2:1622063/149‑1 (MQ=255)                                                                               ggAAGGTCAGGAGACTGAATTTGCTGCCCCAGCTG‑TTGCA‑GATGGCGTACCCCCACTG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCg                                                            <  2:524184/151‑1 (MQ=255)                                                                                      cAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGcatcat                                                     >  2:957242/1‑151 (MQ=255)                                                                                        ggAGACTAAATTTGCTGACCTAG‑GG‑TTGGAGGATGGCGCACCCCCACGG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGca                                                         <  1:3830516/145‑1 (MQ=255)                                                                                                  tttGCTGCCCCAGCTG‑TTGCA‑GATGGCGTAACCCCACGG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCaaa                                            <  2:2182895/148‑1 (MQ=255)                                                                                                     gCTCACCCGTCTT‑TTGAA‑GATGGCGTACCCCCACTG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                                        <  1:2712207/149‑1 (MQ=255)                                                                                                      cTGCCCCAGTTG‑TTGCA‑GAGGGCGTACCCCCACTG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACAcc                                        <  1:2873272/148‑1 (MQ=255)                                                                                                          cccAGTTT‑TTGCA‑GTTGCCATAACCCC‑CTGCACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCg                                       <  1:1111847/145‑1 (MQ=255)                                                                                                                            ggCGTACCCCCCCTG‑CCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa                     <  1:1609333/147‑1 (MQ=255)                                                                                                                            ggCGTAACCCCACGG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCa                     <  2:1888573/147‑1 (MQ=255)                                                                                                                              ctgcacccctcctcaCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCCCGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAGCCCGTCGGTGATCAGCtt                 <  1:1151086/136‑1 (MQ=255)                                                                                                                                    cccACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt          >  2:1702252/1‑150 (MQ=255)                                                                                                                                            aCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTc  >  1:2998514/1‑151 (MQ=255)                                                                                                                                            |                                                                                                                                                       GCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAAGGTCAGGAGACTGAATTTGCTGACCCAGCTG‑TTGCA‑GATGGCGTAACCCCACTG‑ACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350717‑3351003

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