Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	380,012	1	.→G	13.5%	19.8 ‑3.2 ‑0.3	37	intergenic (‑131/+57)	frmR/yaiO	regulator protein that represses frmRAB operon/outer membrane protein
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (1/4):  ref base (23/9):  total (24/13)
Fisher's exact test strand distribution p-value = 4.23e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.59e-01

CCAGAGCATCAATCTGCCCCCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTT‑GAATACGTTGCAGACTGACGCTCTTCATGCTGC  >  NC_000913/379806‑380217                                                                                                                                                                                                                |                                                                                                                                                                                                                ccAGAGCATCAATCTGCCCCCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑gg‑g                                                                                                                                                                                                             >  1:319281/1‑210 (MQ=255) ccAGAGCATCAATCTGCCCCCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑gg‑g                                                                                                                                                                                                             <  1:314278/210‑1 (MQ=255)                 cccccGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGc                                                                                                                                                                                             >  1:548787/1‑210 (MQ=255)                     cGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATACGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGTGGCATAATTggacagg                                                                                                                                                                                         >  2:446032/1‑205 (MQ=255)                                 aCTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGGCGGGAAGGCCGGCCCg                                                                                                                                                                             >  2:672838/1‑210 (MQ=255)                                  cccGAGTAAGGACCTTTTTCTTCCCTTCCGGAGTACTGGGACTTTCGCCCCTCTACATCTGAATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTTTAGTATATAAAACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCg                                                                                                                                                                             <  2:319267/208‑1 (MQ=255)                                  cTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGCGGTCGGAAAACAAAGAGGAAA‑GGTGGGGGGGGCTAATGGGCAGGGAAGGCAGCCCCg                                                                                                                                                                             >  2:9743/1‑210 (MQ=255)                                          aGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGCATGTGGTCGGAAAACAAAGAGGAAAGGG‑GGGGGGGGCTAATCGGCAGGGAAGGCGTCCCAGGATAgcgg                                                                                                                                                                     >  2:319545/1‑210 (MQ=255)                                           ggACCTTTTTCTTCTCCTCCGGAGTACTGGGCATTTCGAACCTCATAATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGACCCGGATAgcgggc                                                                                                                                                                   <  1:298998/210‑1 (MQ=255)                                             aCCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg                                                                                                                                                                 >  1:405441/1‑210 (MQ=255)                                                tttttCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCGGGATAGCGGGCGGCAg                                                                                                                                                              >  2:177797/1‑210 (MQ=255)                                                         tcttcCGGAGTACTGGGCATTTCGCACCTCATCCTCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCa                                                                                                                                                     <  2:500873/210‑1 (MQ=255)                                                             ccGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAAACGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAt                                                                                                                                                 >  2:54297/1‑210 (MQ=255)                                                                     cTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGt                                                                                                                                         <  1:560838/210‑1 (MQ=255)                                                                       gggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCa                                                                                                                                       >  1:308187/1‑210 (MQ=255)                                                                       gggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCGAATCGGCAGGGAAGGCCGCCCCGGACAGGGGGGGGCAGAAGGACACAAAATTTCCAgggaa                                                                                                                                       >  2:444083/1‑207 (MQ=255)                                                                                gCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGcgg                                                                                                                              >  1:331413/1‑208 (MQ=255)                                                                                        tcatcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATAAAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                                                                                                      <  2:634790/210‑1 (MQ=255)                                                                                         catcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                                                                                                     >  2:365496/1‑210 (MQ=255)                                                                                           tcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                                                                                                   >  1:241696/1‑210 (MQ=255)                                                                                           tcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                                                                                                   >  2:153982/1‑210 (MQ=255)                                                                                            cttgaagcaATATACTTAAGGGGGGTATTTTATATGTCAATGCATACCCCCCTATAGTAAATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGAc                                                                                                                  <  2:548759/204‑1 (MQ=255)                                                                                                    aaTATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAGCAGAATTTCCAGGTCAGACGGGCAGCAAGTTGCAGACCGTTaaaa                                                                                                          >  1:254140/1‑210 (MQ=255)                                                                                                                 ggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg                                                                                             <  1:106305/210‑1 (MQ=255)                                                                                                                      aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCg                                                                                        >  1:259173/1‑210 (MQ=255)                                                                                                                        tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCGGCCAGTTGCAGCCCGTTAAAACCATCGGTTGGGGTTTCGTa                                                                                      >  1:376170/1‑210 (MQ=255)                                                                                                                           ttagggtaaatcaaaCCCCCCTAAAGTAAAAAACAATGGAAAATAAAAATTATCCGGGGAATGCGGGCGGAAAACAAAGGGGAAAGGG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAcc                                                                                    <  1:9744/196‑1 (MQ=255)                                                                                                                                          ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGTCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGCCGGCCGGGCTGGCGGTTCCAGGCCGTTAAAATCATCGGGTGGGGGGCCGTACCAGACTTAACCTGcctt                                                                    >  2:349532/1‑208 (MQ=255)                                                                                                                                                   aGTATATAACACTGAGGAATAAAATTTATCCGGTGAATGTGGGCGGAAAACAAAAAGAAAG‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCgag                                                           <  1:672875/210‑1 (MQ=255)                                                                                                                                                    gTATATAACACTGGAGAATAAAATTTATCCGGGGGATGTGGTCGGAAAACAAAGAGGAAAGGG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCgag                                                           <  1:375196/210‑1 (MQ=255)                                                                                                                                                      aaaaaaCACTGGAAAATAAAATTTATCCGGTGAATGTGGTCGGAAAAAAAAGGGGAAAGGG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGAt                                                         <  1:319559/206‑1 (MQ=255)                                                                                                                                                         taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAg                                                     <  2:254131/210‑1 (MQ=255)                                                                                                                                                                                             ccggtctaCAAAGAGGGAA‑GG‑GGGGGGGGCTAATTGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTT‑GAATACGTTGCAGACTGAc                <  1:498256/203‑1 (MQ=255)                                                                                                                                                                                                    aCAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGGGTCGTACCACACTTTACCTGCCGCCAGCCCGAGATTACGTTGCTCAGTCAGCGTTT‑GAATACGTTGCAGACTGAcagaccgg         >  2:498230/1‑203 (MQ=255)                                                                                                                                                                                                     caaaaggggaa‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTT‑GAATACGTTGCAGACTGACGCTCTTCa        <  2:81035/201‑1 (MQ=255)                                                                                                                                                                                                     cAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGGGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTT‑GAATACGTTGCAGACTGACGCTCTTCa        >  2:266402/1‑210 (MQ=255)                                                                                                                                                                                                     cAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAAAGCGGGGGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGAAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTGGTACCACACTTTCCTTGCGGTCAGCCCGAGATAAGGTTGCTCAG‑CAGGGGTTGGAAAAGGTTGCGGCCTGACGccctcaa        >  1:81039/1‑207 (MQ=255)                                                                                                                                                                                                         gagGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTT‑GAATACGTTGCAGACTGACGCTCTTCGtgct    >  2:511177/1‑210 (MQ=255)                                                                                                                                                                                                           ggAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGATATTAAGTTGCTCACTCAGCGTTT‑GAATACGTTTCAGACTGACGctctcaattcggc  >  2:52294/1‑201 (MQ=255)                                                                                                                                                                                                           ggAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTT‑GAATACGTTGCAGACTGACGCTCTTCAttctgc  >  2:165662/1‑210 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                CCAGAGCATCAATCTGCCCCCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTT‑GAATACGTTGCAGACTGACGCTCTTCATGCTGC  >  NC_000913/379806‑380217

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