Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	458,742	0	T→A	19.4%	13.7 ‑0.9 ‑0.1	37	intergenic (+42/‑146)	clpX/lon	ATPase and specificity subunit of ClpX‑ClpP ATP‑dependent serine protease/DNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (15/14):  total (17/20)
Fisher's exact test strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.86e-01

CATCGTAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTG  >  NC_000913/458547‑458947                                                                                                                                                                                                    |                                                                                                                                                                                                              cATCGTAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGTTGAATAATTAACCATTCCCATACAATTAGGTAACCAAAAATGGGAGATAAAATCCCCCCtttt                                                                                                                                                                                                 >  1:72561/1‑209 (MQ=255)  ttCGTAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCttttt                                                                                                                                                                                                <  2:148899/209‑3 (MQ=255)  aTCGTAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCttttt                                                                                                                                                                                                <  1:55995/210‑3 (MQ=255)    cGTATAAACCGCACTGCTCAGTACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGTTTATCGACGAGTCGGTAATTGATGGTAAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAAtt                                                                                                                                                                                              <  2:315607/210‑1 (MQ=255)    cGTAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttgg                                                                                                                                                                                              <  1:289157/210‑5 (MQ=255)     gTAGAAGCCGAACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttggt                                                                                                                                                                                             <  2:72561/210‑6 (MQ=255)      tAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGAT‑GTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGAGATAAAATCCCCCCtttttggtta                                                                                                                                                                                           <  1:391712/210‑8 (MQ=255)                      cGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTGTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCg                                                                                                                                                                            >  1:22642/1‑210 (MQ=255)                            cATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAAt                                                                                                                                                                      >  2:321399/1‑210 (MQ=255)                            cATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGA‑TTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAAtg                                                                                                                                                                     <  2:148921/210‑1 (MQ=255)                                   gATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACAAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTggggg                                                                                                                                                               <  1:308073/210‑1 (MQ=255)                                              cATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCa                                                                                                                                                    <  1:293184/210‑1 (MQ=255)                                                   gactcgtcgataaccactttttCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATAc                                                                                                                                               >  2:80371/22‑210 (MQ=255)                                                    aGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCGAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACt                                                                                                                                              <  1:174854/210‑1 (MQ=255)                                                    aGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACt                                                                                                                                              <  2:293026/210‑1 (MQ=255)                                                                              aGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGa                                                                                                                    >  1:77033/1‑210 (MQ=255)                                                                               gTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGaa                                                                                                                   >  2:214678/1‑210 (MQ=255)                                                                                    tAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGcaca                                                                                                              >  1:134596/1‑210 (MQ=255)                                                                                         gATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGATTGTGGGGGAAACATCCCCATATACTGACGTACATGTTATTAGATGGCGTGAAGCAAAGTCgt                                                                                                         >  1:266815/1‑210 (MQ=255)                                                                                          aTGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtg                                                                                                        <  1:166570/210‑1 (MQ=255)                                                                                          aTGGTAAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCgtg                                                                                                        <  2:77033/210‑1 (MQ=255)                                                                                                   gCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGAt                                                                                               >  2:340468/1‑210 (MQ=255)                                                                                                            tGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCg                                                                                      <  1:41200/210‑1 (MQ=255)                                                                                                              cTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGa                                                                                    >  2:406829/1‑210 (MQ=255)                                                                                                                 aTTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGGGAATAATTAACCATTCCCATTCAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGAGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAAt                                                                                 <  1:95566/210‑1 (MQ=255)                                                                                                                   ttATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTa                                                                               <  1:406855/210‑1 (MQ=255)                                                                                                                      tGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAAc                                                                            <  1:415165/210‑1 (MQ=255)                                                                                                                          aaGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAag                                                                        <  2:30800/210‑1 (MQ=255)                                                                                                                                          aGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATGTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAAt                                                        <  1:340487/210‑1 (MQ=255)                                                                                                                                           ggCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATc                                                       >  1:360312/1‑210 (MQ=255)                                                                                                                                                 tGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGc                                                 >  1:437428/1‑210 (MQ=255)                                                                                                                                                                      cATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGACCGTCCTGAACGCATTTAAATcc                            >  2:398652/1‑210 (MQ=255)                                                                                                                                                                          cAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCCGt                        <  2:153800/210‑1 (MQ=255)                                                                                                                                                                                   aaCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCCGTATTGCCGCt               >  1:151375/1‑210 (MQ=255)                                                                                                                                                                                    aaattaaaGGGGAGATAAAATCCCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTg              >  1:247711/6‑210 (MQ=255)                                                                                                                                                                                    aCCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCACGTATCGACGcct              >  2:236896/1‑208 (MQ=255)                                                                                                                                                                                                ggATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATgtggtg  >  2:331068/1‑210 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                              CATCGTAGAAGCCGCACTGCTCGATACCATGTACGATCTGCCGTCCATGGAAGACGTCGAAAAAGTGGTTATCGACGAGTCGGTAATTGATGGTCAAAGCAAACCGTTGCTGATTTATGGCAAGCCGGAAGCGCAACAGGCATCTGGTGAATAATTAACCATTCCCATACAATTAGTTAACCAAAAAGGGGGGATTTTATCTCCCCTTTAATTTTTCCTCTATTCTCGGCGTTGAATGTGGGGGAAACATCCCCATATACTGACGTACATGTTAATAGATGGCGTGAAGCACAGTCGTGTCATCTGATTACCTGGCGGAAATTAAACTAAGAGAGAGCTCTATGAATCCTGAGCGTTCTGAACGCATTGAAATCCCCGTATTGCCGCTGCGCGATGTGGTG  >  NC_000913/458547‑458947

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.