Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,785	0	A	G	45.3%	74.8 / 14.3	53	F25S (TTC→TCC)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base A (3/26);  new base G (24/0);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 4.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAA  >  NC_000913/1535645‑1535935                                                                                                                                              |                                                                                                                                                         aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACAcc                                                                                                                                                <  1:89681/151‑1 (MQ=255)  aaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGaaaaaaaa                                                                                                                                                  >  2:600322/1‑146 (MQ=255)            tCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAAtt                                                                                                                                      <  2:874121/149‑1 (MQ=255)              aaCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑CGaagttgc                                                                                                                                   >  1:587123/1‑145 (MQ=255)                 cGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑a                                                                                                                              <  1:423644/151‑1 (MQ=255)                    tGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGC‑CGAAATTCCGCG‑gaa                                                                                                                            >  2:99372/1‑150 (MQ=255)                    tGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑aaa                                                                                                                            <  1:910116/150‑1 (MQ=255)                        gTAAAGTAGCGGATGCGCACATCGCCGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑AGaaattgcgcgcaaacg                                                                                                                          >  1:120300/1‑135 (MQ=255)                          aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑AGAAATGCCGCG‑CAACgggt                                                                                                                       >  1:265086/1‑147 (MQ=255)                          aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGtt                                                                                                                       <  2:265159/149‑1 (MQ=255)                               aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑AGGATTGCCGCG‑CAAAGGGTGTgggg                                                                                                                 >  1:522115/1‑150 (MQ=255)                                       cgcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCGGC‑AGAAATGTCGCG‑GAACGGGTGTGGGGGaaaa                                                                                                            >  2:788132/1‑147 (MQ=255)                                         cacaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCg                                                                                                      <  2:558054/151‑1 (MQ=255)                                            aTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGCGGGAATTGCCGGG‑AAACGGTGTTGGGGAaaaacatc                                                                                                        >  2:527536/1‑143 (MQ=255)                                            aTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCGGC‑AGAATTTCCGCG‑GAAAGGGTGTgggggggaaagaacga                                                                                                     >  1:546705/1‑137 (MQ=255)                                              cGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGt                                                                                                  <  1:706551/150‑1 (MQ=255)                                                  gtcgtcGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑AGAAATGTCGCG‑AAAAGGGTGTGGGGGGAAAGCATCGaagggg                                                                                                >  2:402758/1‑143 (MQ=255)                                                       cgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCt                                                                                         <  2:426387/150‑1 (MQ=255)                                                       cgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCt                                                                                         <  2:189817/150‑1 (MQ=255)                                                        gCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCCGC‑AGGAATGCCGCG‑CAAAAGGTGTGGGGGAAAAACAACGACGCGGGACACCtcc                                                                                       >  1:197879/1‑149 (MQ=255)                                                            tgCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGt                                                                                    <  1:284113/150‑1 (MQ=255)                                                               agcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑AGAAATTCCGCG‑GAAAGGGTGTGGGGGAAAAAAAACGACGGGGTACACCCGCGGgtt                                                                                  >  2:261023/1‑149 (MQ=255)                                                               agcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCa                                                                                <  2:245238/151‑1 (MQ=255)                                                               agcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCa                                                                                <  1:668485/151‑1 (MQ=255)                                                                gcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCAt                                                                               <  1:610648/151‑1 (MQ=255)                                                                 cagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑AGAAATGCCCCG‑GAACGGGTGTGGGGGAAAAACAACGACGTGGTACAACTGTGGGTTCAt                                                                               >  2:837954/1‑150 (MQ=255)                                                                   gAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAAAACCCGC‑AGAATTGCCGCG‑AAAAGGGTGTGGGGGAAAAACATCGACGTGGGACACCTGCGGGTTCAtt                                                                              >  2:33920/1‑149 (MQ=255)                                                                      aaGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGC‑AGAATTGCCGCG‑GAACGGGTGGGGGGGAAAAACAACGACGGGGGACAACCGCGGGTTCATTaa                                                                            >  1:215248/1‑148 (MQ=255)                                                                       aGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGtt                                                                        <  2:220376/151‑1 (MQ=255)                                                                             ttCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCt                                                                   <  1:626918/150‑1 (MQ=255)                                                                               ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACGGC‑CGAAATGCCGCG‑GAACAGGTGTGGGGGAAAAGCAAAGAAGCGGGACACCCGCGGGTTTATTAACGTTCTCCcgt                                                                 >  1:433890/1‑150 (MQ=255)                                                                                 gAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGC‑AGAAATGCCGCG‑GAAAGGGTGTGGGGGGAAAACATCGACGGGGGGGACCTCCGGGTTTTTTAAAGTTCTTCtgtg                                                                >  1:718789/1‑149 (MQ=255)                                                                                         tAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGG‑‑AGAAATGCCGCG‑GAAAGGGTGGGGGGGAAAAAAAAAGAAGGGGGAAACCCGCGGGTTTAATAAAGGTCTCCTGTGAGTGGGGaa                                                       >  2:334037/1‑150 (MQ=255)                                                                                         tAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑AGAATTTCCGCG‑GAACGGGTGGGGGGGAAAAACATCGACGGGGTACACCCGCGGGGTTATTAACGGTCTCCTTTGACTggggga                                                       >  2:791759/1‑146 (MQ=255)                                                                                              cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑AGAATTGGCGCG‑GAACGGGTGTGAGGGAAAAAAAACGACGGGGGACACCTGCGGGTTTATTAAAGTTCTCCTGTGACTGGGGAACTa                                                    >  2:834282/1‑148 (MQ=255)                                                                                                cgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATa                                               <  2:480945/151‑1 (MQ=255)                                                                                                 gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGGAAAAACACCCGG‑AGACATGCCGCG‑GAACGGGGGTGAGGTAAAAAGATCCACGCGGGATACCTGCGGGTTTGTTAACGTTCTCCtgtgtacgggggactat                                                   >  1:798014/1‑133 (MQ=255)                                                                                                    cGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCt                                           <  2:401271/151‑1 (MQ=255)                                                                                                     gATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTg                                          <  2:691299/151‑1 (MQ=255)                                                                                                       tAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCCGC‑AGAAATTCCGCG‑GAACGGGTGTGGGGGGAAAAAAACGACGGGGGACACCCGCGGGTTTTTTAAAGTTTTCCTGTTACTTGGGAAATATCATaaac                                            >  2:270514/1‑144 (MQ=255)                                                                                                       tAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGC‑AGGAATTTCGCG‑AAAAGGGTGTGGGGGAAAAACAAACACGCGGGACACCCGCGGGTTCATTAAAGATTTCCCGTGGCTGGGGGAATAACATAgac                                            >  1:752396/1‑145 (MQ=255)                                                                                                        aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                       <  1:480364/151‑1 (MQ=255)                                                                                                         aGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGG‑AGGATTTTCGCGCAAACGG‑TGGGGGGTGAAAAAATCGAAGGGGGACACCTGCGGGTTCCTTAACATTCTCCCCTGAATGGAGGAATATCATa                                               >  1:950665/1‑142 (MQ=255)                                                                                                           cTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAACACCCGG‑GGAAATTCCGCC‑AAACGGGTGTGGGGGCAAAACCTCGCCGTGGAAGACCCGCGGGTTTTCTAACGCTTTCCTGTGTCCGGGGAACAATTATAACCTGccaa                                      >  2:767935/1‑146 (MQ=255)                                                                                                             gCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTgg                                   <  1:611427/150‑1 (MQ=255)                                                                                                                  tGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGC‑CGAATTGCCGCG‑CAAAGGGTGTGGGGGAAAAACATCGACGTGGGACACCTCCGGGTTTCTTAACGGTCTCCTGTGACTGGGGGACTATCATAACCCGCCAGTgggc                                 >  1:911729/1‑145 (MQ=255)                                                                                                                     gcTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGc                          <  1:799052/151‑1 (MQ=255)                                                                                                                         cGCTAAGATTATCGGC‑GGGGGAAACAACACCGGC‑AGAAATGCCGCG‑AAAAGGATGTGGGGGAAAAGCAACGACCGGGGACACCCGCGGGTTTATTAAAGCTCTCCTGTGGCTGGGGAAATATCATAACCCGCAAGTGGGCGGggag                           >  1:899703/1‑146 (MQ=255)                                                                                                                                  tATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATcc              <  2:731589/150‑1 (MQ=255)                                                                                                                                     cggc‑ggGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGcc          <  1:784896/151‑1 (MQ=255)                                                                                                                                     cggc‑ggGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGc           <  2:223699/150‑1 (MQ=255)                                                                                                                                            ggAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg    <  2:718928/151‑1 (MQ=255)                                                                                                                                            ggAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg    <  1:49641/151‑1 (MQ=255)                                                                                                                                              aaaaCAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaa  <  1:286071/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                         AAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGC‑AGAATTGCCGCG‑AAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAA  >  NC_000913/1535645‑1535935

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