Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,717,014	0	T	G	43.2%	71.9 / 13.0	45	G156G (GGT→GGG)	yiaD	multicopy suppressor of bamB; outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (8/17);  new base G (19/0);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 1.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTA  >  NC_000913/3716865‑3717162                                                                                                                                                         |                                                                                                                                                       aTCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGggtg                                                                                                                                                       <  2:492014/151‑1 (MQ=255)      cAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGGTT‑GGGTTTACCCACAGCACGGGGGGGCAc                                                                                                                                                  >  1:124377/1‑151 (MQ=255)      cAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCa                                                                                                                                                   <  1:482093/150‑1 (MQ=255)       aaTATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AAAGGGGTT‑TGTTTTAACGACAGCCCgggggg                                                                                                                                                      >  2:864967/1‑146 (MQ=255)               gAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGAT‑AACGTGCCT‑GCGTATAACGTCAGCACCGGCGGGGGCGA‑CCCGaa                                                                                                                                         >  1:822153/1‑150 (MQ=255)                aaCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGT‑TAAGTGATT‑TGTTATACCCACAGCACGGGGGGGCCCCA‑CCCGaa                                                                                                                                         >  2:796679/1‑149 (MQ=255)                   aaTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGc                                                                                                                                    <  2:210340/151‑1 (MQ=255)                             tCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTcc                                                                                                                            <  2:12696/149‑1 (MQ=255)                              cGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAg                                                                                                                         <  2:724531/151‑1 (MQ=255)                                  agcagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑TACGTGATT‑GGTTATACCCACAGCACGGGGGGGCACGA‑CCTGAAAATGCGGCTCTCCCCGCaa                                                                                                                      >  1:952236/1‑150 (MQ=255)                                       cagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGT‑AAAGGGGTT‑GGTTATACCGACAACACGGGGGGGGCCGA‑CCCGAACATGCGTCTCcccccccaaactgcg                                                                                                                >  2:369726/1‑140 (MQ=255)                                          cgcgACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGT‑AACCGGGTT‑GGGTTTAACCACAGCACCGGGGGGCACGA‑CCCGAACATGGCGCTCTCCCAGCAACG‑GGCgg                                                                                                              >  1:849333/1‑149 (MQ=255)                                            cgACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTc                                                                                                          <  1:683927/151‑1 (MQ=255)                                             gACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTT‑AAC‑TGATT‑GGTTTTACCCACAGCACGGGGGGGCACGA‑CCCGAACATGGGTCTCTCCCAGCAAAC‑GGCGGATTcc                                                                                                         >  1:761766/1‑151 (MQ=255)                                              aCCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAACGGGATT‑GGTT‑TACCCACAGCACCGGGGGGCCCGA‑CCCGAACATGCGGCTCTCCCCGCCAAG‑GGCgggatcc                                                                                                         >  1:584779/1‑145 (MQ=255)                                               cccTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGt                                                                                                       <  2:589767/151‑1 (MQ=255)                                                     aaCCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGGTTGGGTTGACCCGAC‑GCACGGGGGGTCAAGACCCTGAA‑ATTCGCCTCTTCCCGCAAAGTTGCGGA‑TCCCTTGCCac                                                                                                 >  1:233118/1‑150 (MQ=255)                                                       ccGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAgcg                                                                                               <  1:172443/151‑1 (MQ=255)                                                        cggcggGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAgcgc                                                                                              <  1:931306/151‑1 (MQ=255)                                                                gcTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGGGATT‑GGGTATAACGACAGCACGGGGGGGCACGA‑CCTGAACATGCGTATCTCCCAGCAACG‑GGCGGATTCCGTTGCCAGCGCGTTGATc                                                                                       >  2:309829/1‑150 (MQ=255)                                                                  tAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGGTT‑GGTTTTACCGACAACACGGGGGGTCACGA‑ACTGAAAATGCGGCTCTCCCAGCAAAC‑TGGGGATTTCGGTGCCAgcgccgttgtcacc                                                                                    >  1:210017/1‑141 (MQ=255)                                                                  tAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGGTATACCGACAGCACGGGGGGGCACGA‑CCTGAACATGCGTCTCTCCCCGCAACG‑GGCGGATTCCGGTGCCAGCGCGTTTATCAcc                                                                                    >  1:252398/1‑151 (MQ=255)                                                                    acacCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCa                                                                                  <  1:475940/151‑1 (MQ=255)                                                                    acacCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGGGATT‑GGTTATACCGACAGCACGGGGGGGCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCAccc                                                                                   >  2:758375/1‑150 (MQ=255)                                                                           gACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTg                                                                           <  2:838393/151‑1 (MQ=255)                                                                            aCCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AAAGTGATT‑GGGTATACCGACAGCACGGGGGGGCACGA‑CCTGAACATGCGTCTCTCCCCGCAACG‑TGCGGGTTCCCTTGCCCGCGCGTTGGTCACCCAgggggggg                                                                          >  2:453190/1‑146 (MQ=255)                                                                             ccGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGGTT‑GGTTATACCGACAGCACGGGGGGTGCCGA‑CCTGAAAATGCGTCTCTCCCCGCAACG‑TGCGGATTTCCCTGGCACCGCGTTGATCACCCAggggctgg                                                                          >  2:611579/1‑145 (MQ=255)                                                                                  gTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AAAGTGATT‑GGTTATAACGACAGCACGGGGGGGCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGGCGcc                                                                     >  2:201915/1‑150 (MQ=255)                                                                                   tGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AAAGTGATT‑GGGTATACCGACAGCACGGGGGGTCACGA‑CCTGAACATGCGTCTCCCCCCGCAACG‑TGGGGATTTCGGTGCCAGGGCGGTGGTCACCCCGGGGGTGGACGcc                                                                     >  2:495082/1‑149 (MQ=255)                                                                                     gCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGcc                                                                 <  1:886227/151‑1 (MQ=255)                                                                                       aaTGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGc                                                               <  2:876700/151‑1 (MQ=255)                                                                                         tGGTACTGAAAGAGTATCCGAAAACGGC‑GGGT‑AACGGGATT‑GGTTTTACCGACAGCACGGGGGGGCACGA‑CCTGAAAATGCGTCTCTCCCAGCCACG‑TGCGGATTTCGTTGCCCACGCGGTGATCACCCCCGGCGGGGGCGcccacccc                                                               >  1:961141/1‑143 (MQ=255)                                                                                           gTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑TGGTATAACGACAGCACGGGGGGGCACGA‑CCTGAAAATGCGTCTCTCCCAGCAAAG‑GGCGGATTCCGTTGCCAGCGCGGTGATCACCCAGGGCGTGGACGCCAGCCGCATcc                                                           >  2:538579/1‑151 (MQ=255)                                                                                                    gagTATCCGAAAACGGC‑GGTT‑AAAGTGATT‑GGTTATATCGACAGCACGGGGGGGCACGA‑CCCGAACATGCGTCTCTCCCAGCAACG‑TGCGGGTTTCGTTGCCCGCGCGTTGATTACCCCGGGCGTGGGCGCCAGCCGCCTCCGCTCTCAgg                                                  >  1:772938/1‑151 (MQ=255)                                                                                                          ccGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGGCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCCTCCGTACTCAGGGCCt                                              >  2:247456/1‑149 (MQ=255)                                                                                                              aaaCGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGccc                                        >  1:134216/1‑151 (MQ=255)                                                                                                              aaaCGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGccc                                        <  1:897332/151‑1 (MQ=255)                                                                                                                 cggc‑ggTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCgg                                      >  2:784681/1‑150 (MQ=255)                                                                                                                  ggc‑ggTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGc                                     >  1:640061/1‑150 (MQ=255)                                                                                                                  ggc‑ggTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCt                                    <  2:757646/151‑1 (MQ=255)                                                                                                                    c‑ggTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTaa                                  <  2:952408/151‑1 (MQ=255)                                                                                                                         t‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAAccc                               >  1:858771/1‑150 (MQ=255)                                                                                                                                      gTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCa                    >  2:852145/1‑150 (MQ=255)                                                                                                                                          taCCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGc               >  1:564049/1‑151 (MQ=255)                                                                                                                                           aCCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCa              <  2:451589/151‑1 (MQ=255)                                                                                                                                            ccGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAc             >  2:504757/1‑151 (MQ=255)                                                                                                                                                   cACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGaa      <  1:796530/151‑1 (MQ=255)                                                                                                                                                       ggtggtCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTa  <  1:864811/151‑1 (MQ=255)                                                                                                                                                        gtggtCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTa  >  2:522888/1‑150 (MQ=255)                                                                                                                                                        gtggtCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTa  >  1:178835/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                       ATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTT‑AACGTGATT‑GGTTATACCGACAGCACGGGTGGTCACGA‑CCTGAACATGCGTCTCTCCCAGCAACG‑TGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTA  >  NC_000913/3716865‑3717162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑