Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	1,344,036	0	G→C	13.2%	10.2 0.0 ‑0.2	38	T104S (ACT→AGT)	yciT	global regulator of transcription; DeoR family
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (17/16):  total (19/19)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.77e-01

CCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCG  >  NC_000913/1343837‑1344244                                                                                                                                                                                                         |                                                                                                                                                                                                                 ccagccaTCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAg                                                                                                                                                                                                        <  2:169758/210‑1 (MQ=255)   agacaTCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCGGACTCTGGGCGAgca                                                                                                                                                                                                      <  1:318197/207‑3 (MQ=255)     ccaTCAATACCAATAATTGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTACTCGCCCGGACTCTGGGCGAGCAga                                                                                                                                                                                                    <  2:91954/210‑2 (MQ=255)     ccaTCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCGGACTCTGGGCGAGCAga                                                                                                                                                                                                    <  1:106890/210‑2 (MQ=255)            tACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGACGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGTGCATtgc                                                                                                                                                                                             >  2:268612/1‑210 (MQ=255)              ccAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctg                                                                                                                                                                                           >  2:319770/1‑210 (MQ=255)                            gCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGAt                                                                                                                                                                             <  2:432546/210‑1 (MQ=255)                                             ggATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCg                                                                                                                                                            <  2:34213/210‑1 (MQ=255)                                                   aCTGGCGTGACAAAGTGCCAACCATACTTTCGCTTTTTTTCTGATAACCGCCAACGAGAAAAATAACTTCACAAGGCGCGTCTTTCAGAAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGaa                                                                                                                                                      <  2:429572/210‑1 (MQ=255)                                                             cAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGAcgc                                                                                                                                            >  1:468502/1‑210 (MQ=255)                                                               aaGGGCCAAACATACTTTCGCTTTTTTTCTGATACACGCCACCGAACAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCgg                                                                                                                                          <  2:266766/210‑1 (MQ=255)                                                                        ccATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGCCATTTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGACAAAGATGGTTTCGCCCGGTTGAACCATCGACGCGGCAAACTCg                                                                                                                                  >  1:503272/1‑210 (MQ=255)                                                                          aTACTTTCGCTTTTTTTATGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCg                                                                                                                               <  1:199143/210‑1 (MQ=255)                                                                                   ttgtttatcagataaacacgcccGAGAAAAATCAATTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGc                                                                                                                      <  1:41568/189‑1 (MQ=255)                                                                                   cTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGc                                                                                                                      >  2:21966/1‑210 (MQ=255)                                                                                       ttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCGTTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCa                                                                                                                  >  1:81477/1‑210 (MQ=255)                                                                                           cTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGGGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCCAGCTTCAGCGt                                                                                                              >  2:509678/1‑210 (MQ=255)                                                                                               tACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAg                                                                                                          >  2:226383/1‑210 (MQ=255)                                                                                                 cacGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGtt                                                                                                        >  2:75265/1‑210 (MQ=255)                                                                                                    gCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCt                                                                                                     <  1:412972/210‑1 (MQ=255)                                                                                                      cACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCtca                                                                                                   >  2:447670/1‑210 (MQ=255)                                                                                                             cAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATCCg                                                                                            >  1:145215/1‑210 (MQ=255)                                                                                                              aaaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACgg                                                                                           <  2:99960/210‑1 (MQ=255)                                                                                                                         acaAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtc                                                                                <  2:468533/210‑1 (MQ=255)                                                                                                                           aaGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcat                                                                              <  1:463936/210‑1 (MQ=255)                                                                                                                                      ttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGc                                                                   <  1:463516/210‑1 (MQ=255)                                                                                                                                        cAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGa                                                                 >  1:475099/1‑210 (MQ=255)                                                                                                                                          gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGaaa                                                               >  1:537688/1‑210 (MQ=255)                                                                                                                                                       tcttccgatctccGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCAtg                                                  <  1:423496/199‑1 (MQ=255)                                                                                                                                                            tgctgACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCAc                                             >  1:158671/1‑210 (MQ=255)                                                                                                                                                                  ccGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCAtgagatcggaaga                                       >  2:423525/1‑199 (MQ=255)                                                                                                                                                                                 cA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGttttt                        >  1:33340/1‑210 (MQ=255)                                                                                                                                                                                    tttttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCg                      <  1:12046/210‑1 (MQ=255)                                                                                                                                                                                    tttttCTGCGGCTCGCCCAGGGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCg                      <  1:85976/210‑1 (MQ=255)                                                                                                                                                                                      tttCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAg                    <  1:75263/210‑1 (MQ=255)                                                                                                                                                                                          ccctacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAgtgtg                >  1:213718/6‑210 (MQ=255)                                                                                                                                                                                           cctacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAgtgtgt               >  1:319534/5‑210 (MQ=255)                                                                                                                                                                                              gCTCGACCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTCCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTga            >  1:329610/1‑210 (MQ=255)                                                                                                                                                                                                        agTCCGGGCGGGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCg  <  1:226376/210‑1 (MQ=255)                                                                                                                                                                                                         |                                                                                                                                                                                                                 CCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACA‑TTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCG  >  NC_000913/1343837‑1344244

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