Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,429,984	0	T	G	43.4%	73.0 / 12.0	53	L41V (TTG→GTG)	cycA	D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base T (11/19);  new base G (23/0);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 8.36e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTC  >  NC_000913/4429839‑4430133                                                                                                                                                   |                                                                                                                                                        cagacagGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGttt                                                                                                                                                    <  2:732658/151‑1 (MQ=255) cagacagGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGttt                                                                                                                                                    >  2:562553/1‑151 (MQ=255)  agacagGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTa                                                                                                                                                   <  2:556001/151‑1 (MQ=255)    acagGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCCTTTGGACGGGGTTGttt                                                                                                                                                    >  1:394100/1‑148 (MQ=255)                agaaAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGGGTTTATGGGGGCTGGCa                                                                                                                                       >  1:677573/1‑149 (MQ=255)                  aaaaaaaCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTTTGGGGTCTGGCaaaa                                                                                                                                    >  2:261855/1‑150 (MQ=255)                        aCATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGG‑GCCATTGGTAAGGGGTTGTTTATGGGGTCTGGCAAAAAGa                                                                                                                                 >  1:65223/1‑147 (MQ=255)                         cATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTTTTGGGTCTGGCAAAACGA‑tt                                                                                                                              >  1:357612/1‑148 (MQ=255)                           tGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTTTTGGGTCTGGCAAAACGA‑TTAGcc                                                                                                                          >  2:233038/1‑150 (MQ=255)                            ggTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGTGGCCTTTGGTA‑GGGGGTGTTTTTGGGGTCTGGCAAAAAGG‑TTAGCCtt                                                                                                                        >  1:643640/1‑151 (MQ=255)                              tAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGGACGGGGGTGTTTTTGGGGTCTGGCAAAACGA‑TTAGCCtt                                                                                                                        >  2:681875/1‑149 (MQ=255)                                gATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCCTTGGTACGGGGTTGTTTATTGGGTCTGGCAAAACGA‑TTAGCCTTGcc                                                                                                                     >  2:131308/1‑150 (MQ=255)                                 tacAGGTAAAAGTCGTTGCCGAGGATAAGGCTCCGGCTGAACAGACGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GACTGTGCGAAGGGGGGGTTTCTGGGGTGTGGAAAAACGG‑TTAGCCtctca                                                                                                                     >  1:510709/3‑145 (MQ=255)                                    aGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGcc                                                                                                                <  1:280209/151‑1 (MQ=255)                                    aGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGA‑TTAGCCTTGCCggg                                                                                                                  >  2:230825/1‑149 (MQ=255)                                     ggTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCg                                                                                                               <  1:74716/151‑1 (MQ=255)                                       tAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTGACAAGCCGACATATTCATCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGcc                                                                                                                <  1:435775/148‑1 (MQ=255)                                         aaaGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTATGGGGGCTGGCAAAAAGA‑TTAGCCTTGCCGGGGCGTCg                                                                                                            >  2:306868/1‑150 (MQ=255)                                          aaGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGAt                                                                                                          <  1:712467/151‑1 (MQ=255)                                              cGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCCTTGGTACGGGGGTGTTTTTGGGGGCTGGGAAAAAGA‑TTAGCCTTTCCGGGCCGTCGattatt                                                                                                      >  1:322990/1‑151 (MQ=255)                                                ttGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCCTTGGTAAGGGGGTGTTTTTCGGGTCTGGCACAACGC‑TTTGCCTTGCCGGGCCGCCGCTctcttt                                                                                                    >  1:264158/1‑146 (MQ=255)                                                      atgatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTTGGTCGGGGTTGTTTCTGGGGGTTGGCCAAACGG‑TTAGGCTTGCCGGGCCGGCGATCATTTTCGttt                                                                                               >  2:259753/1‑150 (MQ=255)                                                      atgatCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTta                                                                                              <  1:214732/151‑1 (MQ=255)                                                         atCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCCTTGGTACGGGGGTGGTTATGGGGTCTGGCAAAACGG‑TTAGCCTTGCCGGGCCGTCGGTCATTTTCGttt                                                                                               >  1:236832/1‑147 (MQ=255)                                                                   ggCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCAATGGTACGGGGTTTTTTTTGGGGTCTGGCCAAACGAGTTGCCCTTGCC‑GGCCGTCGGTCCTTTTCTGTTTTTT‑GATCCtt                                                                                   >  2:651605/1‑148 (MQ=255)                                                                   ggCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGACGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGt                                                                                >  2:772812/1‑151 (MQ=255)                                                                        aaCAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTAt                                                                           <  2:861558/151‑1 (MQ=255)                                                                         aCAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATTATTGGtttttt                                                                           >  1:330637/1‑148 (MQ=255)                                                                          cAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGc                                                                         <  2:494434/151‑1 (MQ=255)                                                                               gCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCttt                                                                    <  2:370773/151‑1 (MQ=255)                                                                                cTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTTGTTTATAT‑GATCATTGGTTTTTTGCTCtttt                                                                   >  2:324911/1‑151 (MQ=255)                                                                                cTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGGACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCtttt                                                                   >  2:85504/1‑151 (MQ=255)                                                                                cTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCACTCGTACGGGGGTGTTTTTCGCGGCTGGCACAACGA‑TTTGGCTTCCCGGGGCCACGATCCTTTTCGTTTATATGGTTCTTTGGTTTT‑TGGTCtttt                                                                   >  1:91001/1‑151 (MQ=255)                                                                                   cGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTTTGGGGGCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCttttt                                                                  >  2:494473/1‑149 (MQ=255)                                                                                    ggCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGt                                                               <  2:770032/151‑1 (MQ=255)                                                                                          aTCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTTGTACGGGGGTGTTTTTGGGGTCTGGCAAAACGA‑TTAGCCTTTCCGGGCCGTCGATCATTTTCGTTTTTAT‑GATCATTTGTTTTATGCTCTTTTTTGTGAt                                                            >  1:638395/1‑148 (MQ=255)                                                                                              cacaAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCggg                                                       <  2:104322/149‑1 (MQ=255)                                                                                              cacaAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCg                                                         <  1:256610/147‑1 (MQ=255)                                                                                              cacaAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCa                                                     <  1:369756/151‑1 (MQ=255)                                                                                                caAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCCTTGGTACGGGGGTGTTTTTGGGGTCTGGGAAAACGA‑TTAGCCTTTCCGGGCCGGCGATCATTTTTGTTTATAT‑GGTCCTTGGTTTTTTGCTCTTTTTTGTGATGCGGGcga                                                    >  2:239883/1‑148 (MQ=255)                                                                                                 aaaCCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCCTTTGTACGGGGGTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGGTTATTTTCGTTTATAT‑GATCATTTGTTTTTTGCTCTTTTTTGTGGTGCGGGCaa                                                    >  1:290578/1‑149 (MQ=255)                                                                                                   aCCGACATATTCAGCTGCTTTCCATTGGCGGT‑GCCATTGGTCCGGGGTTGTTTATGGGATCTGGCTAAACGA‑TTAGCCTTGCCGGGCCGTCGCTAATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCCTGATGCgggaa                                                     <  1:702979/146‑3 (MQ=255)                                                                                                    ccGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGGCGTCGGTCATTTTCGTTTATAT‑GATAATTGGTTTTATGCTCTTTTTCGTGATGCGGGCCATgggg                                               >  1:752894/1‑151 (MQ=255)                                                                                                     cGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGGTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGGAATggggg                                              >  2:50466/1‑151 (MQ=255)                                                                                                      gACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGa                                             <  1:602163/151‑1 (MQ=255)                                                                                                                    aTTGCCATTGGCGGG‑GCCATTGGGACGGGGGTGTTTATGGGGTCTGGCAAAACGC‑TTAGCCTTGCCGGGGCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTTGTGATGCGGGCAATGGGGGAATTGCTGtttt                                  >  1:684497/1‑148 (MQ=255)                                                                                                                      tGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGaa                              <  2:661485/150‑1 (MQ=255)                                                                                                                       gCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATc                            <  1:324874/151‑1 (MQ=255)                                                                                                                             ggCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGaa                       <  1:206920/150‑1 (MQ=255)                                                                                                                                 gT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACa                   >  1:375968/1‑150 (MQ=255)                                                                                                                                         tGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCtttt           >  2:642956/1‑151 (MQ=255)                                                                                                                                           gTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCtttt           >  2:789439/1‑149 (MQ=255)                                                                                                                                              cGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTg       >  2:336533/1‑150 (MQ=255)                                                                                                                                                  gttgttTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTc  <  1:815844/151‑1 (MQ=255)                                                                                                                                                   ttgttTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTc  <  1:450985/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                        CAGACAGGTACAGGAAGAAAAAAACATGGTAGATCAGGTAAAAGTCGTTGCCGATGATCAGGCTCCGGCTGAACAGTCGCTACGGCGCAATCTCACAAACCGACATATTCAGCTTATTGCCATTGGCGGT‑GCCATTGGTACGGGGTTGTTTATGGGGTCTGGCAAAACGA‑TTAGCCTTGCCGGGCCGTCGATCATTTTCGTTTATAT‑GATCATTGGTTTTATGCTCTTTTTCGTGATGCGGGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTC  >  NC_000913/4429839‑4430133

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