Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,574,051	0	G	T	46.2%	80.6 / 28.2	52	V384V (GTC→GTA)	pqqL	putative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (5/23);  new base T (24/0);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 3.00e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573910‑1574201                                                                                                                                              |                                                                                                                                                       cGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCg                                                                                                                                               <  2:907349/151‑1 (MQ=255)   aaaGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGa                                                                                                                                              <  2:319370/150‑1 (MQ=255)   aaaGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTAaatc                                                                                                                                                   >  1:231202/1‑145 (MQ=255)       ttGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTATATCATCGAggttt                                                                                                                                         >  2:232497/1‑147 (MQ=255)       ttGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGCAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCTTCATCGAGtttt                                                                                                                                         >  1:441269/1‑149 (MQ=255)            atGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTca                                                                                                                                      <  2:90551/149‑1 (MQ=255)                   ttCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGTTAGATTTTATAATATTGAAGTTTTCAGCCGaaaa                                                                                                                              >  2:589004/1‑150 (MQ=255)                    tCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGaaaa                                                                                                                              <  2:983639/149‑1 (MQ=255)                      cGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTTAGCAGaaaaa                                                                                                                             >  1:450941/1‑148 (MQ=255)                       ggCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATCAACGAGTTCTTTCGCAGaaaaaac                                                                                                                           >  2:267941/1‑147 (MQ=255)                         cGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAGTTTGCATCATCTAGGTTTTCATCAaaaaaaaccag                                                                                                                        >  2:1014164/1‑145 (MQ=255)                           aCAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTAAATCATCGGGTTCTTTTACAGaaaaaaaaagctg                                                                                                                     >  1:18829/1‑143 (MQ=255)                            cAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGGTCTTCTGCAGAAAAAACaatct                                                                                                                      >  1:1060303/1‑145 (MQ=255)                               gAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTACATCATCGAGGTTTTTTGCCGAAAAAACATGCTGAt                                                                                                                   >  1:913112/1‑149 (MQ=255)                                     cGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTAACTCATCGAGATCTTTAGCAGAAAAAACaaagcgaacaataggt                                                                                                           >  2:707354/1‑136 (MQ=255)                                       tattatTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATCAACGAGTTTTTTAGCAGAAAAAACACTCTGATCAAAGGt                                                                                                            >  2:1166539/1‑148 (MQ=255)                                           attTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATTATCGAGGTCTTTAGCAGAAAAAACATGCTTATTaa                                                                                                                >  1:41196/1‑140 (MQ=255)                                            tttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTCATCACAAAAAACATGCTGATCAATGGTTGCCAAc                                                                                                     >  2:535976/1‑150 (MQ=255)                                                aTGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTCTTTATCAGAAAAAACATGCTGATCAAAGGTTTCCaa                                                                                                      >  2:117296/1‑145 (MQ=255)                                                    gCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCAtt                                                                                            <  2:971035/151‑1 (MQ=255)                                                       ggATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAAt                                                                                         <  1:262587/151‑1 (MQ=255)                                                           gCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCAt                                                                                     <  1:589170/151‑1 (MQ=255)                                                           gCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGCGGTTTTTCGCAGAAAAAAAATGCTGATCAATGGGTGCCCAATCTGCCCtt                                                                                            >  2:796729/1‑144 (MQ=255)                                                             cAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTc                                                                                   <  1:335137/151‑1 (MQ=255)                                                              aGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCAACGAGTTTTTTAGCAGAAAAAACAAGGTGTTTAATGGGTGGCAACACTTCCATTAATGCAtt                                                                                    >  2:824632/1‑149 (MQ=255)                                                               ggCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTATATCATCGAGTTTTTTAGCCGAAAAAACCTGCTGATTAATTGTTTTCAACCCTTCCATTAATAGATTctc                                                                                 >  2:616410/1‑149 (MQ=255)                                                                gCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCgca                                                                                <  2:733415/151‑1 (MQ=255)                                                                 cGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTACATCATCGGAGGCTTTAGCAGAAAAAACATGCTGGTCAAAGGTTGCCAACTCTGCCATTAAAGCAtt                                                                                    >  1:737056/1‑146 (MQ=255)                                                                    cTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATTATCGAGTTTTTCAGCAGAAAAAACAAGCTGATCAATGGGTGCCCACCCTGCCCTTTATTCATTCgc                                                                                 >  2:360589/1‑146 (MQ=255)                                                                                cGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCAtattat                                                                <  1:796980/151‑1 (MQ=255)                                                                                    aaTCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGtc                                                            <  2:286013/151‑1 (MQ=255)                                                                                      tCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGTTCTTCAGCAGAAAAACCAAGCTGATCAATAGTTGCCCACTCTTCCATTAATGCATTTGCAGCATCCTGCATATTATCGt                                                             >  1:94027/1‑148 (MQ=255)                                                                                          gTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATTAACGAGGTTTTTAGCCGAAAAAACAATCTGATCAAAGGTTGCCCACTCTTCCCTTAATGCATTCGCAACATCCTGCATATTTTCGt                                                             >  2:1086341/1‑144 (MQ=255)                                                                                          gTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTTTTCAGCAGAAAAAACATTCTGATCAATGGTTGCCAACCCTTCCCTTTATGCATTCGCCGCATCCTGCAtatttttgt                                                             >  1:518000/1‑139 (MQ=255)                                                                                          gTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGc                                                      <  2:570002/151‑1 (MQ=255)                                                                                           ttCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGAGGTCTTTAGCAGAAAAAACAAGCTGATCAATGGTTGCCAACTCTGCCCTTAATGCATTCGCAGCATCCTGCCTATTTTCGTCTCCTg                                                       >  1:409608/1‑149 (MQ=255)                                                                                             cAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCAtt                                                   <  1:232544/151‑1 (MQ=255)                                                                                             cAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCAt                                                    <  1:820466/150‑1 (MQ=255)                                                                                                  tGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTaca                                               <  2:957488/150‑1 (MQ=255)                                                                                                    cTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACg                                             <  1:473331/150‑1 (MQ=255)                                                                                                      gATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACgg                                            <  2:18825/149‑1 (MQ=255)                                                                                                      gATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACgg                                            <  2:450812/149‑1 (MQ=255)                                                                                                             cGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGc                                   <  1:808531/151‑1 (MQ=255)                                                                                                                      cAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATaa                          <  1:707572/151‑1 (MQ=255)                                                                                                                       aGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAAt                         <  2:1059982/151‑1 (MQ=255)                                                                                                                                  gCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAtt              <  1:634884/151‑1 (MQ=255)                                                                                                                                    gggTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAtttt            >  1:544922/1‑151 (MQ=255)                                                                                                                                     ggTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTa           <  1:411279/151‑1 (MQ=255)                                                                                                                                       tAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAAc         >  2:524746/1‑151 (MQ=255)                                                                                                                                        aGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACg        >  1:706340/1‑151 (MQ=255)                                                                                                                                          aTTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCt      >  1:524012/1‑151 (MQ=255)                                                                                                                                           tttGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTg     >  1:341665/1‑151 (MQ=255)                                                                                                                                              gACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcgc  <  1:1087722/151‑1 (MQ=255)                                                                                                                                              gACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcg   <  2:429581/150‑1 (MQ=255)                                                                                                                                              gACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcg   <  1:947333/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                       CGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATCAATGGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573910‑1574201

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