Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	380,013	0	G→.	12.8%	15.7 0.0 ‑0.1	47	intergenic (‑132/+56)	frmR/yaiO	regulator protein that represses frmRAB operon/outer membrane protein
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (32/9):  total (37/10)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.38e-01

TGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGA  >  NC_000913/379875‑380202                                                                                                                                            |                                                                                                                                                                                               tGGGCATTTCGCGCCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑gg‑g                                                                                                                                                                                             <  1:140767/141‑1 (MQ=255)       tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCgg                                                                                                                                                                              >  4:654779/1‑150 (MQ=255)             aCCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGa                                                                                                                                                                        <  4:220970/150‑1 (MQ=255)                 catcatCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑gg‑gggggtgg                                                                                                                                                                                      >  3:420221/1‑132 (MQ=255)                       cTGCATGCAATATACTATAGGGGGGTATTCTATATGTAAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                                                                                                                              <  4:576421/150‑1 (MQ=255)                                  ataCTATAGGGGGGTATTCTATATGTCAATGCATACCTCCCTATAGCATATAACACTGGAGAATAAAATTTATCCGGTGCATGTGGTCGG‑AAACAAAGTGGAAA‑GC‑GGGGGGGGCTAATCGGCAGCGAAggg                                                                                                                                                                    >  3:85893/1‑131 (MQ=255)                                  ataCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                                                                    >  2:310335/1‑149 (MQ=255)                                  ataCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                                                                    <  4:310304/149‑1 (MQ=255)                                     cTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGgccg                                                                                                                                                                  >  3:96739/1‑132 (MQ=255)                                       ataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgc                                                                                                                                                       >  1:569772/1‑141 (MQ=255)                                       ataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgc                                                                                                                                                       >  1:217216/1‑141 (MQ=255)                                            ggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTCTCCGGTGAATGTGGTCGGTAAACACAGAGGAAC‑GA‑GGGGGGGGATACTCGGAAGGGAAGGCCGCCCCGGATAGAGGGCGGCAGAGGGa                                                                                                                                         >  4:558864/1‑150 (MQ=255)                                            ggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGTGGCTAATCGGCAGGGAAGGCCGCCCAGGATAGCGGGCGGCAGAAGGa                                                                                                                                         >  4:496471/1‑150 (MQ=255)                                                 aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑gg                                                                                                                                                                                               <  2:664010/92‑1 (MQ=255)                                                 aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑gg                                                                                                                                                                                               >  4:663943/1‑92 (MQ=255)                                                  ttCTATATGTCAAAGCATAACCCCCTATAGTATAGAAAACTGGAGAATAAAATTTATTCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGAGGGGGGGCCTAATCGGCAGGGAAGgccgaaagggattgcg                                                                                                                                                      >  3:315268/1‑120 (MQ=255)                                                    cTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgg                                                                               >  2:561695/1‑200 (MQ=255)                                                              aTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑G‑GGGGGGGGCTAATCGGCAGGGAAGCCCGCCCCGGATAGCGGGCGGCAGAAGGCATCAGAAttt                                                                                                                               >  1:36622/1‑141 (MQ=255)                                                                 cATACCCCCCTATAGGATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                                                                             <  1:96749/141‑1 (MQ=255)                                                                     ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATTGGCAGGGAAGGCCGCCCCGGATAGCGGTCGGCAGAAGGAATCAGAATTTACAGGTCAGACGGGATGCAAGTTGCAGCCCGTTAAAATCATCGGTTAGGGTGTCGTACGACACtt                                                              >  2:533499/1‑200 (MQ=255)                                                                     ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                                                                                                                                  >  3:379487/1‑132 (MQ=255)                                                                     ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGt                                                                                                                         >  1:198408/1‑141 (MQ=255)                                                                      gccccTATAGTATGTAACCCTGGCGAATAAAGTTTACCCAGTGACTATGATCGGAAAACCAAGAAGAAG‑GA‑GGGGGGGACTAATCGGCAGGGCAAACCGCCCGAGATAGCGGACGGCAGAAAGAATCAGAATTTCCAGGTc                                                                                                                        <  1:44926/140‑1 (MQ=255)                                                                       ccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATATCGGGCGGCAGACGGAATCAGAGTGTCCAGGTCAGAAGGGCTGCAAGTTTCAGACCGTTAAAATAATCGGTTGGGGTGTCGTACCACACTTTa                                                            >  2:694030/1‑200 (MQ=255)                                                                       ccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                                                                                >  3:480944/1‑132 (MQ=255)                                                                       ccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                                                                                >  3:6281/1‑132 (MQ=255)                                                                       ccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAACCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTa                                                            >  2:110205/1‑200 (MQ=255)                                                                       ccccTATAGTATATAACACTGCAGAATAAAATCTATCTGGCGAATGCGGTCGGAAAACAAAGAGGATA‑GG‑GGGGGGGGTTAATCGGTAGGGAAGGTCGCTCCGGATAGTGCGCGGCTGAAGGAATCAGAAtt                                                                                                                                >  3:44920/1‑132 (MQ=255)                                                                       caccTATAGTATATAAAACTGGAGAATAAAATTTATCCGGTGCATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGTGCTAATCAGCGGGGAAGGCCGCTCCGGAAAGCGGGCGGAAGAAGGAACCAGAATTTCAAGGTCAGATGGGCTg                                                                                                              >  4:138708/3‑150 (MQ=255)                                                                         ccTATAGTAGATAACGCTGGATAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCGGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAcc                                                          <  2:432030/200‑1 (MQ=255)                                                                          cTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGCAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                                                                             >  3:516846/1‑132 (MQ=255)                                                                          cTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt                                                         >  2:270944/1‑200 (MQ=255)                                                                          cTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt                                                         >  2:602459/1‑200 (MQ=255)                                                                          cTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑G‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                                                                          >  4:231458/1‑150 (MQ=255)                                                                            ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑G‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAgg                                                                                                                          >  3:542922/1‑132 (MQ=255)                                                                               gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                                                                                                                        >  3:356672/1‑132 (MQ=255)                                                                               gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt                                                    >  2:360257/1‑200 (MQ=255)                                                                               gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATAGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt                                                    >  2:21798/1‑200 (MQ=255)                                                                               gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑G‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc                                                   >  2:191947/1‑200 (MQ=255)                                                                               gTATATAACACGGGAGAATAAAATTTATCCGGTGAAAGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                                                                                      <  4:54532/150‑1 (MQ=255)                                                                                 atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAGA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGCCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCAAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACAACACTTTACCTGCCGTCa                                                  >  2:282310/1‑200 (MQ=255)                                                                                 atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATCGCTGGCGGCAGAAGGAATCAGACTCTCCAGGCCAg                                                                                                                      >  3:294639/1‑132 (MQ=255)                                                                                 atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                                                                             >  1:109277/1‑141 (MQ=255)                                                                                 atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg                                                                                                    >  4:345745/1‑150 (MQ=255)                                                                                              gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                                                                         >  3:251097/1‑132 (MQ=255)                                                                                                   aaaTTTATCCGGTGACTGTGGTCGGAAAACAAAGAGGAAA‑‑G‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGCATTTCCAGGTCAGACGGGCTGCCAGTTGCATACCGTTaaaa                                                                                          >  1:80588/1‑141 (MQ=255)                                                                                                     aTTTATCCGGTGAATGTGGTCGGAAAACAAAGATGAAA‑GG‑GGGGGGGGCTAATCGGCATGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                                                                             <  3:553970/105‑1 (MQ=255)                                                                                                     aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                                                                             >  1:554003/1‑105 (MQ=255)                                                                                                                          ggAAAAAAAAGAGGAAA‑‑G‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAcc                                                          <  4:454290/150‑1 (MQ=255)                                                                                                                                cAAAGAGGAAAGGG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACt                                                               >  1:209607/1‑141 (MQ=255)                                                                                                                                 aaaGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGa  <  2:461988/200‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                                                               TGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GG‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGA  >  NC_000913/379875‑380202

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