Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,991	0	G	T	21.6%	142.8 / 10.4	52	L207F (TTG→TTT)	pdeL	DNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base G (15/25);  new base T (10/1);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 2.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTTC  >  NC_000913/332853‑333127                                                                                                                                           |                                                                                                                                         tctGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTgg                                                                                                                                         >  1:983783/1‑140 (MQ=255)  ctGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGGTTTTTTGc                                                                                                                                        >  2:867044/1‑140 (MQ=255)   tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCa                                                                                                                                       >  1:656584/1‑140 (MQ=255)   tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCa                                                                                                                                       >  1:277643/1‑140 (MQ=255)    ggCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAg                                                                                                                                      <  2:752813/140‑1 (MQ=255)     gCGGAGTCATCCGGTCTTATTGTCATAATTACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGCTTGTTTTTTGGCAGc                                                                                                                                     >  2:770627/1‑140 (MQ=255)          gTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCggggt                                                                                                                                >  2:1130813/1‑138 (MQ=255)            cATCCGGTCTTATTGTAATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATGTCCATATTGGCATCATCGTCTCGGCGGGTTGTTTTTTGGCAGCGGGAttt                                                                                                                              <  2:520976/140‑1 (MQ=255)                    cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                      <  2:277643/140‑1 (MQ=255)                            cATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGa                                                                                                              <  2:983783/140‑1 (MQ=255)                             aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGGGTCTGaa                                                                                                             >  1:481004/1‑140 (MQ=255)                             aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGaa                                                                                                             >  2:909054/1‑140 (MQ=255)                                    cccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGTGGGATTTGAAAAAGAGTGTCTGAACCTGGtt                                                                                                      <  2:512152/140‑1 (MQ=255)                                    cccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGtt                                                                                                      <  1:19214/140‑1 (MQ=255)                                     ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCTGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta                                                                                                     >  2:227858/1‑140 (MQ=255)                                     ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAaga                                                                                                                      >  2:623408/1‑123 (MQ=255)                                     ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                      <  1:623408/123‑1 (MQ=255)                                     ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta                                                                                                     <  2:950504/140‑1 (MQ=255)                                     ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTta                                                                                                     <  2:621089/140‑1 (MQ=255)                                     ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGGTTTTTTGGAGCGGGGTTTGAAAAAAAGAGTCTGAAACTGGttt                                                                                                     >  2:32998/1‑139 (MQ=255)                                        ccAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGTGTCTGAACCTGGGtaata                                                                                                  >  2:903951/1‑140 (MQ=255)                                               atgaAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAAtatgt                                                                                           >  1:1069356/1‑136 (MQ=255)                                                      aGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTGGCAGCGGGATTTGAAAAAAAGGGGCTGGAACCGGTTAATAAATTAGGTAAAGAt                                                                                    >  2:913999/1‑140 (MQ=255)                                                       gACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATa                                                                                   <  1:134762/140‑1 (MQ=255)                                                         cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTAAAAAGAGTGTCTGAAACTGGTTAATAAATTAGGTAACGATaaa                                                                                 >  2:134762/1‑140 (MQ=255)                                                          tGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTT‑GCAGCGGGGTTTGAAAAAGAGGGTCTGAACCTGGGTAATAAATTAGGTAACGATaaaaa                                                                               >  2:295497/1‑139 (MQ=255)                                                            cGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGactgtctc                                                                              >  1:960478/1‑133 (MQ=255)                                                            cGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATc                                                                              >  1:540500/1‑140 (MQ=255)                                                                tatTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGc                                                                          <  1:1130813/140‑1 (MQ=255)                                                                  ttCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTg                                                                        >  1:627601/1‑140 (MQ=255)                                                                      gATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAAAGGGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtc                                                                    >  1:1069238/1‑140 (MQ=255)                                                                       aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGGTTTGAAAAAGAGAGTCTGAACCTGGTTAATAAATTAGGGAAAGATAAAATCAAACTGGTtct                                                                   >  1:1090842/1‑140 (MQ=255)                                                                       aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtct                                                                   <  1:475461/140‑1 (MQ=255)                                                                       aTGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGGTTTTAAAAAAAGGGTCTGAACCTGGTTAATAAATTAGGGAAAGATAAAATCAAGCTGGTtct                                                                   >  1:252491/1‑140 (MQ=255)                                                                             gTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCt                                                             >  1:1081114/1‑140 (MQ=255)                                                                               aaaaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTaa                                                           <  1:685807/140‑1 (MQ=255)                                                                                 aaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACg                                                         <  1:713830/140‑1 (MQ=255)                                                                                      ttGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGTCTGAACCTGGGTAATAAATTAAGGAACGATAAAAACAAACTGGTTCTCGAGCTAACGGAACg                                                    >  2:686307/1‑140 (MQ=255)                                                                                      ttGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACg                                                    <  1:903951/140‑1 (MQ=255)                                                                                      ttGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACg                                                    <  2:540500/140‑1 (MQ=255)                                                                                       tgctgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGGGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGt                                                   >  1:1097911/1‑140 (MQ=255)                                                                                            ccGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGGGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAAccc                                              >  1:990152/1‑140 (MQ=255)                                                                                                       ccATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTaa                                   <  2:68958/140‑1 (MQ=255)                                                                                                       ccATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTaa                                   <  2:327922/140‑1 (MQ=255)                                                                                                         atatTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGTGTCTGAACCTGGGTAATAAATTAGGGAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAAACCTATTCCGGGAAcc                                 >  2:536187/1‑139 (MQ=255)                                                                                                           atTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGGTTTGAAAAAGAGGGTCTGAACCTGGTTAATAAATTAGGGAAAGATAAAATCAAACTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGGAAAGcc                               >  1:366747/1‑140 (MQ=255)                                                                                                              ggCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGa                            >  1:401398/1‑140 (MQ=255)                                                                                                              ggCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAAAGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGa                            >  2:1084206/1‑140 (MQ=255)                                                                                                                cATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTTATAAATTAGGGAAAGATAAAATCAAGCTGGTTTTCGAGCTAACGGAACGTAACCCTATTTCGGTAAAGCCAGaaa                          >  2:324431/1‑139 (MQ=255)                                                                                                                cATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGGGGCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAg                          >  1:969347/1‑140 (MQ=255)                                                                                                                 aTCAACGTCTCGGCGGGTTG‑TTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGcc                        >  1:112714/1‑140 (MQ=255)                                                                                                                   cAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCa                       <  1:686307/140‑1 (MQ=255)                                                                                                                      cGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCagag                    <  2:677139/140‑1 (MQ=255)                                                                                                                      cGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCagag                    <  2:252491/140‑1 (MQ=255)                                                                                                                                ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa          <  1:909054/140‑1 (MQ=255)                                                                                                                                ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa          <  1:32998/140‑1 (MQ=255)                                                                                                                                ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa          <  1:324431/140‑1 (MQ=255)                                                                                                                                ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa          <  1:227858/140‑1 (MQ=255)                                                                                                                                ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa          <  2:1069238/140‑1 (MQ=255)                                                                                                                                        tttGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTtc  >  1:598152/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                         TCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTTC  >  NC_000913/332853‑333127

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