Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	4,060,264	0	A→C	15.0%	13.6 ‑0.4 ‑0.6	40	intergenic (+34/‑183)	typA/yihL	GTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (16/18):  total (18/22)
Fisher's exact test strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.81e-01

ACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACC  >  NC_000913/4060057‑4060467                                                                                                                                                                                                                |                                                                                                                                                                                                            aCTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTTCCAGCGATg                                                                                                                                                                                                           >  1:271622/1‑210 (MQ=255)      cGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                                                                                      >  1:265745/1‑210 (MQ=255)               tgttCTGGTTCCGCCTATCCGAAGGGCTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                                                                                                                             <  2:440618/210‑4 (MQ=255)                  tCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                                                                                                                                          >  2:73168/1‑210 (MQ=255)                  tCTAGTTCCGCCTATCCGCATGCCTCTGTAACAAGCTCTGGAGTTCATCAACGCTGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGTAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgact                                                                                                                                                                                          <  1:35310/210‑7 (MQ=255)                           gCCTATCCGCATGACTCTGTAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGAACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                                                                                                                                 <  1:38736/210‑1 (MQ=255)                                         ctctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTTAGCGCACACCAATGCg                                                                                                                                                                   >  2:359675/1‑210 (MQ=255)                                           ctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                                                                                                                                 <  1:38203/210‑1 (MQ=255)                                                    gCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                                                                                                                        >  1:132330/1‑210 (MQ=255)                                                                aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                                                                                                                            >  1:418318/1‑210 (MQ=255)                                                                      gatgaCGCACTGGTAGAAGTCGCTCCTACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                                                                                                                      <  2:10764/210‑1 (MQ=255)                                                                         gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                                                                                                                                   <  1:377505/210‑1 (MQ=255)                                                                           ttccctacgcgacgctcttCCGATCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                                                                                 <  2:263976/192‑1 (MQ=255)                                                                                  gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                                                                                                                          <  1:361002/210‑1 (MQ=255)                                                                                            cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAATTCTTTGGTTACGTCATTAGAGACCTTTTTCACaa                                                                                                                >  1:259938/1‑210 (MQ=255)                                                                                                        tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAg                                                                                                    <  2:319518/210‑1 (MQ=255)                                                                                                          cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTc                                                                                                  >  2:386059/1‑210 (MQ=255)                                                                                                           gTATTCGTAAACGTCACCTGCCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                                                                                 <  1:120094/210‑1 (MQ=255)                                                                                                           gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                                                                                 >  1:239434/1‑210 (MQ=255)                                                                                                                        tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAtt                                                                                    <  1:386085/210‑1 (MQ=255)                                                                                                                                      aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg                                                                      <  2:20705/210‑1 (MQ=255)                                                                                                                                      aCGATCGTCGCCGCGCAAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTgg                                                                      <  2:418291/210‑1 (MQ=255)                                                                                                                                           cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaa                                                                 >  1:284444/1‑210 (MQ=255)                                                                                                                                                                acGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGt                                            <  1:255097/209‑1 (MQ=255)                                                                                                                                                                     aGACGATTAATTTCGTTTTTTAGTCATAAAAGACCTGACAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTaaa                                       <  2:371326/210‑1 (MQ=255)                                                                                                                                                                        cgcttAATTTCTTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTGTTTCCTGTCGTGCGCGCACACCAATGAGAGTATCTCCCTTATAAGTCTGTGGTTGCGTCAATAGAGAGCTTTGTCACAATTATCTGCACAATCTTACGCACTTAATTGCTTTCTTTTTTTGCGTAATCGTAAGATGCTTCTTCTGGTCTAATCCa                                    >  1:68929/4‑210 (MQ=255)                                                                                                                                                                        cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCACAGTCATACGCCGTGAATTGCTTTCTTTTTTTGCGTAAGCGTAAGATGCTTCATCTGGTTCAATCCa                                    >  1:68936/1‑210 (MQ=255)                                                                                                                                                                        cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTAAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAATCATACGCCGTTAATTGCTTTCTTTTATGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCa                                    >  1:87288/1‑210 (MQ=255)                                                                                                                                                                                       tttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAAt                     <  2:471614/210‑1 (MQ=255)                                                                                                                                                                                        ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATg                    <  1:359698/210‑1 (MQ=255)                                                                                                                                                                                             cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGga               >  1:314561/1‑210 (MQ=255)                                                                                                                                                                                             cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGga               >  2:87769/1‑210 (MQ=255)                                                                                                                                                                                              aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGgag              >  1:413280/3‑210 (MQ=255)                                                                                                                                                                                                aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTGTGTGGTTACGGCAATAGAGAGCTTTGTCACAATTATCTGCAAAGGCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGGTTAAACCAAAAGGATTAAACAATGGCGGAGaa            >  1:48464/1‑210 (MQ=255)                                                                                                                                                                                                 aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGaat           <  2:265731/210‑1 (MQ=255)                                                                                                                                                                                                    aaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGaatcaa        <  2:413253/210‑1 (MQ=255)                                                                                                                                                                                                     cccTGCCATCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGaatcaat       <  2:314546/209‑1 (MQ=255)                                                                                                                                                                                                     aCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGaatcaat       <  1:87770/210‑1 (MQ=255)                                                                                                                                                                                                        tGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTTGTTTAAACCAAAAGGATTAAACAATGGCGGAGAAACAATACa    >  1:117710/1‑210 (MQ=255)                                                                                                                                                                                                          ccATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCAcc  <  1:14622/210‑1 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                            ACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACC  >  NC_000913/4060057‑4060467

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 28 ≤ 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.