Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,079,468	0	A	G	44.9%	65.3 / 17.3	50	D55G (GAT→GGT)	putP	proline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (1/26);  new base G (22/0);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.63e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACA  >  NC_000913/1079320‑1079614                                                                                                                                                     |                                                                                                                                                            ccGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGat‑a                                                                                                                                                          <  2:1058319/151‑1 (MQ=255)            gACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCggctg                                                                                                                                               <  2:204751/150‑1 (MQ=255)               aTTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGT‑A‑TGAg                                                                                                                                                     >  1:1036331/1‑141 (MQ=255)                ttttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGT‑C‑TGatggggcgggcgt                                                                                                                                          >  2:920029/1‑139 (MQ=255)                ttttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGTGGGGGGT‑G‑TGAGCGggtg‑gctg                                                                                                                                          >  2:537628/1‑150 (MQ=255)                 tttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGT‑A‑TTAGGgggtg‑g                                                                                                                                             >  1:64193/1‑146 (MQ=255)                   tgtgTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGG‑A‑TGAGCgggtg‑ggggt                                                                                                                                         >  1:85728/1‑144 (MQ=255)                    gtgtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AAt                                                                                                                                   <  1:340124/151‑1 (MQ=255)                        cTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGA‑A‑TGAGGGGCTG‑GGTG‑TT‑TATg‑ggg                                                                                                                              >  1:294783/1‑151 (MQ=255)                          atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATg‑ggg                                                                                                                              <  2:123332/149‑1 (MQ=255)                          atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGtt                                                                                                                            <  1:150983/151‑1 (MQ=255)                          atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGt                                                                                                                             <  2:128763/150‑1 (MQ=255)                                   gCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGgcg‑c                                                                                                                  <  2:811054/151‑1 (MQ=255)                                     aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGG‑A‑TAAGCGGGTGCGTGG‑TT‑AGTG‑GGGTTGTCggcg                                                                                                                     >  2:162572/1‑146 (MQ=255)                                        atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGT‑A‑TGAGCGGGTG‑GGGG‑TT‑TATT‑GGGGTGCCggggggtggtt                                                                                                              >  1:329774/1‑143 (MQ=255)                                            tgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTGGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTtttc                                                                                                        <  1:967334/150‑1 (MQ=255)                                             gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTtttctt                                                                                                      <  2:935237/151‑1 (MQ=255)                                               ttgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGAGCGTCGGGA‑A‑TGAGAGGCGG‑GCTG‑GT‑GATG‑GGGGTGCCGGGCG‑G‑GG‑Gttttttttt                                                                                                    >  1:424992/1‑150 (MQ=255)                                                  ggTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGT‑A‑TGAGCGGGTG‑GGTG‑TT‑AATG‑GGGGTGCCGGGGG‑G‑TG‑Ttttttttt                                                                                                     >  1:985668/1‑147 (MQ=255)                                                    tttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGG‑A‑TGGGGGGGTG‑GGCG‑TT‑TTTG‑GGGGTGTCGGGGG‑G‑TG‑Tttttttttttcggg                                                                                               >  1:698751/1‑145 (MQ=255)                                                      tATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCCGGT‑ATTGAG‑GGGTG‑GCCG‑TT‑TATT‑GTGGTGTCGGGCG‑CGGG‑TTTTTTTTTCCtgggt                                                                                              >  2:634776/1‑149 (MQ=255)                                                       aTCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGAttt                                                                                            <  1:968542/151‑1 (MQ=255)                                                       aTCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGT‑A‑TTAGCGGGTG‑GGTG‑TT‑TATT‑GGGGTGCCGGGGG‑C‑TG‑TTTTTTTTTTCCGGGttt                                                                                            >  1:265446/1‑151 (MQ=255)                                                               gCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGT‑A‑TTAGGGGCTG‑GGTG‑GT‑TATG‑GGGGTTCCGGGGG‑G‑GGTTTTTTTTTTCCGGGTTTTcgcaa                                                                                       >  1:386413/1‑145 (MQ=255)                                                                cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCt                                                                                   <  2:902648/151‑1 (MQ=255)                                                                  aTCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTgg                                                                                 <  2:1013734/151‑1 (MQ=255)                                                                    cAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGCCGGGG‑A‑TTAGGCGCTG‑GCTG‑GT‑ATTGTGGGGGGCCCGCGG‑G‑GG‑TTTTTTTTTCCGGGATTTCGGaaactggt                                                                                 >  2:95777/1‑145 (MQ=255)                                                                      aCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCg                                                                             <  1:724956/151‑1 (MQ=255)                                                                      aCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCg                                                                             <  2:724465/151‑1 (MQ=255)                                                                      aCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGG‑A‑TTAGGGGCCG‑GGGG‑TT‑TATG‑GGGGTGCGCCGCGCC‑GG‑TTTTTTTTTCTGGGTTTTCCCaaacgtg                                                                                  >  2:560615/1‑143 (MQ=255)                                                                         aaaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCa                                                                          <  2:849368/151‑1 (MQ=255)                                                                           aaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCAtt                                                                        <  2:30133/151‑1 (MQ=255)                                                                           aaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGGCGCT‑A‑TGGGCGGGGG‑GGGG‑TTGAATG‑GGGTGCGCGGGGC‑C‑GG‑TTTTTTTTTCCGGGGTTTCCGAAGGC‑GGGTCGCCCtt                                                                        >  1:1015480/1‑151 (MQ=255)                                                                              cTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGT‑A‑TGAGCGGGTG‑GCGGTTT‑TTTG‑GGGTGGCCGGGGG‑T‑GG‑TTTTTTTTTCGGGGTTTTCC‑CAAGCTGGGTCGCCCtt                                                                        >  1:833565/1‑148 (MQ=255)                                                                                 tGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCt                                                                   <  2:424990/150‑1 (MQ=255)                                                                                  gacgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTGT‑TGGGGGGGGG‑GTGG‑TT‑ATTG‑GGGTGGCGGGGGG‑T‑GT‑TTTTTTTTTTCGGGGTTTCC‑AAAGCTGGGGCGGCATTGGCCt                                                                   >  2:939326/1‑149 (MQ=255)                                                                                                gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGcgcg                                                      <  1:1047308/148‑1 (MQ=255)                                                                                                gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGGGGGT‑A‑TTAGGGGCTG‑GGGG‑TT‑TATG‑GGGGTGGCGGGCG‑C‑CG‑GTTTTTTTTTCGGGGTTTTCGGAAGCTGGGTGGCCCTTGGCCTGTGCTTTGGGGCGTgg                                                   >  2:94387/1‑151 (MQ=255)                                                                                                gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGAGGTGGGG‑A‑TGAGGGGCGG‑GGTG‑TT‑TATT‑GGGGTGCCGGGCG‑G‑TG‑TTTTTTTTTCTCGGGTTTTCGAAAACGGGATCGGCCTTTGCCCGCCAtttagcg                                                        >  1:315420/1‑141 (MQ=255)                                                                                                          aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTgg                                          <  1:627716/150‑1 (MQ=255)                                                                                                          aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGG‑A‑TGAGGGGGGG‑GGTG‑TT‑TATT‑GGGGTGGCGGGGG‑C‑GG‑TTTTTTTTTTCGGGGTTTTCGAAAAATGGGTCGCCCATTGCCCGACACTAGGGGCGGGGGTTTAATgg                                          >  1:530800/1‑150 (MQ=255)                                                                                                            tCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGT‑A‑TTAGGGGCCG‑GGTG‑GT‑TATG‑GGGTTTGCGGGGG‑C‑GG‑GTTTTTTTTTCGGGGTTTTCCAAAAATGGGTTGCCCTTTGCCCGACCTTAAGGGCGGGGGtttaatt                                           >  2:610601/1‑142 (MQ=255)                                                                                                               tGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAg                                       <  1:1067363/148‑1 (MQ=255)                                                                                                                  gCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggtg                                 <  1:33976/151‑1 (MQ=255)                                                                                                                  gCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggtg                                 <  1:696482/151‑1 (MQ=255)                                                                                                                        cGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCgg                            <  2:431699/150‑1 (MQ=255)                                                                                                                                  ttATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgt                  <  1:228413/150‑1 (MQ=255)                                                                                                                                        gcggGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATa           <  2:944591/151‑1 (MQ=255)                                                                                                                                             tGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGa       >  1:50966/1‑150 (MQ=255)                                                                                                                                                 tCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGACGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGCAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGaa      <  2:624128/147‑1 (MQ=255)                                                                                                                                                  cGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATaca  <  2:1047325/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                            CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGAT‑A‑TGAGCGGCTG‑GCTG‑TT‑AATG‑GGGTTGCCGGGCG‑C‑TG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACA  >  NC_000913/1079320‑1079614

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: ‑