Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,842,925	0	G	A	41.7%	28.4 / 28.8	49	K185K (AAG→AAA)	gdhA	glutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base G (1/27);  new base A (20/0);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.25e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGAAGGGATGCG  >  NC_000913/1842785‑1843071                                                                                                                                              |                                                                                                                                                            aGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑TGAAAAAAATCTCcaac                                                                                                                                                    >  1:1029574/1‑150 (MQ=255)      aGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGA‑GGAAAAAACACTCccacaa                                                                                                                                                  >  1:1116038/1‑147 (MQ=255)       ggTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCg                                                                                                                                             <  1:871249/151‑1 (MQ=255)       ggTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCg                                                                                                                                             <  1:858286/151‑1 (MQ=255)       ggTGAAATGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAAGACCg                                                                                                                                             <  2:365877/151‑1 (MQ=255)                aTGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGGAAAAAAACCCCCACAATATCGCCCG‑G‑GTC‑t                                                                                                                                 >  1:193079/1‑151 (MQ=255)                                cgcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGGAAAAAATCTCCCACAAAAACCCCCG‑C‑GTC‑ttttcccggaaagg                                                                                                                    >  1:265593/1‑136 (MQ=255)                                           gAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTggcg                                                                                                   <  2:547222/151‑1 (MQ=255)                                            aaCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑TGGGAAAAGTCTCCCACAAAAACCCCCG‑C‑GGC‑TTT‑A‑CCCGG‑TAGGGGCGTTCAttttgtgg                                                                                                  >  2:590304/1‑146 (MQ=255)                                             aCTGTATCGCCCCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAGGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTggcg                                                                                                   <  1:915051/149‑1 (MQ=255)                                              cTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGA‑TGGAAAAAACCTCCCACAAAACCGCCCG‑C‑GTC‑TTC‑C‑CCCGG‑GAGGGGCCTTCATTTtgcg                                                                                                   >  1:8523/1‑148 (MQ=255)                                                       cACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTc                                                                                      <  1:1066230/151‑1 (MQ=255)                                                        aCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑TGGAAAAAATCTCCCACAAAAACCCCCG‑C‑GTCTTTC‑C‑CCGGT‑AGGGGCCTTTATTTTggggggggctttttt                                                                                        >  1:1031383/1‑140 (MQ=255)                                                        aCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑GGGAAAAAAACTCCAACAAAAACGCCCG‑C‑GCC‑CTCTC‑CCGGG‑AGGGGCCTTT‑ATTTTGCGGCAG‑TCTTAttt                                                                                      >  2:1013110/1‑149 (MQ=255)                                                         ccTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGGAAAAACTCTCCCACCATACCACCTC‑C‑GTC‑CTT‑A‑CCAGTAAAGGGCCTTTCttttttcgcacgcctatt                                                                                         >  1:927680/1‑134 (MQ=255)                                                         ccTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGc                                                                                    <  2:406341/151‑1 (MQ=255)                                                          cTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAACACTTCAACAAAAACGCCCG‑CGGTC‑TTT‑C‑CCGGG‑CGGGCCCTTTC‑TTTTGCGGCGG‑TCTTATTCTcc                                                                                   >  1:871909/1‑151 (MQ=255)                                                           tGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGG‑GGGAAAAAACCTCCTACAAAAACCCCCG‑C‑GTT‑CTC‑T‑CCCGG‑AAGGGGCCTTCATTTTGCGGGGG‑TCTTATTcttcc                                                                                  >  2:401754/1‑147 (MQ=255)                                                           tGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGA‑TGAAAAAAAACTCCCACCAAAACGCCCG‑C‑GTC‑TTT‑A‑CCGGG‑AAGGGGCCTTCATTTGGGGGGAG‑GCTTTTTcttcc                                                                                  >  2:445352/1‑147 (MQ=255)                                                               cgcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAAAACTCCCACAAAACCGCCCG‑C‑GTC‑TTC‑T‑CCCGG‑AAGGGGCTTTCATTTTGGGGGGG‑GCTTATTTGccccg                                                                                >  1:979933/1‑147 (MQ=255)                                                                gcgGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTAAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGaa                                                                              <  2:663347/150‑1 (MQ=255)                                                                  gaaacccgccgTTCCGGCAGTTGATATCGGGTTTGTTGGTCTGGAACTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCt                                                                           <  2:964640/146‑1 (MQ=255)                                                                      aCCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑c‑c                                                                      <  1:848953/150‑1 (MQ=255)                                                                       ccGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑Cgg                                                                    <  1:165599/151‑1 (MQ=255)                                                                       ccGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑Cg                                                                     <  1:652644/150‑1 (MQ=255)                                                                        cGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAGGA‑GGAAAAAAAACTCCAACAAAAACCCCCG‑C‑GTC‑TTC‑A‑CCCGG‑AAGGGGCCTTTATTTTGGGGCGG‑GGTTTTTTTCCCCGGAGCTA‑C‑Cgg                                                                    >  1:999873/1‑150 (MQ=255)                                                                             ttCCGGCAGGTGATATCGAGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACg                                                               <  1:77160/150‑1 (MQ=255)                                                                               ccGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGTCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGt                                                             <  2:545702/150‑1 (MQ=255)                                                                               ccGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACgg                                                              <  2:259232/149‑1 (MQ=255)                                                                                  gCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTgg                                                         <  2:981550/151‑1 (MQ=255)                                                                                   cAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGAGTAAAAAAAATCACC‑ACAAAAACACCGC‑C‑GGC‑TTT‑T‑CCAGT‑CAAGGGCTTTCTTTTTTGGGCAG‑GCTTTTTTGCCCGGGAGCTA‑CGCGGCt                                                                  >  2:1020501/1‑142 (MQ=255)                                                                                      gTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTa                                                     <  2:372951/151‑1 (MQ=255)                                                                                      gTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTa                                                     <  2:1116008/151‑1 (MQ=255)                                                                                        gATATCGGGGTTGGTGTTCTTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAACCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTAtt                                                   <  2:921434/151‑1 (MQ=255)                                                                                        gATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGA‑TGAAAAAAATCTCCAACAATAACGCCCG‑C‑GGC‑TTC‑C‑CCCGG‑AAGGGGCTTTTATTTTGCGGCCG‑GCTTATTTGCCCCGAGGCTA‑C‑CGGCCACGGGCTGGGTtttt                                                   >  1:723860/1‑151 (MQ=255)                                                                                         atatCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAAACGCCCCACAAAACAGGCTGCC‑GAC‑TT‑‑A‑CCCGA‑CAGGGGCTTTCCTTTTGCGGCAGATGTCTTTTCCCCGGACG‑TA‑C‑AGGCTACGGCCTGGGTTAttt                                                  >  2:751794/1‑151 (MQ=255)                                                                                           atCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTc                                                 <  2:154747/150‑1 (MQ=255)                                                                                            tCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTca                                                <  1:211578/150‑1 (MQ=255)                                                                                               gggTTGGTGGTCGTGAAGTCGGCTTTATGGGGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGa                                            <  2:906551/151‑1 (MQ=255)                                                                                                  ttGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGG‑TGAAAAAAATCTCCCACAAAACCGCCCG‑C‑GTC‑TTC‑C‑CCCGG‑AAGGGGCTTTTTTTTTGCGGCGG‑TCTTTTTTTCCCCGGAGCCA‑C‑CGGCTACGGGCTGGGTTTTTTTTCAGa                                            >  1:538300/1‑148 (MQ=255)                                                                                                       ggtCGTGAAGTCGGCTTTATGGCGGGGGGGA‑TGAAAAAACACTCCCAAAAAAAACCCCC‑C‑GTC‑TTT‑C‑CCGGG‑AAGGGCCCTTTCTTTTGCGGCAG‑GCTTTTTTTCCCGGGAGGGA‑C‑CCGCCCCGGGCTGGGttttttttcaaaaacaa                                       >  2:408364/1‑137 (MQ=255)                                                                                                          cGTGAAGTCGGCTTTATGGCGGGGATGG‑GGAAAAAAAATTCCAAAAAAAACGCCCG‑C‑GGC‑TTC‑ATCGGGT‑AGGGGGCTTT‑TTTTTGCGGCGG‑TCTTTTTTTCCCGGAAGGTACC‑CGGCT‑CGGGCTGGGTTTTTTTTCAGAAGCaacact                                   >  2:597745/1‑145 (MQ=255)                                                                                                          cGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTa                                  <  1:895292/150‑1 (MQ=255)                                                                                                              aaGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCACCAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACg                             <  2:708321/151‑1 (MQ=255)                                                                                                                  cGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCa                          <  2:297303/150‑1 (MQ=255)                                                                                                                         aTGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATgg                  <  1:1013136/151‑1 (MQ=255)                                                                                                                          tGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATggg                 >  1:886585/1‑151 (MQ=255)                                                                                                                            gCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGtt               <  2:1098383/151‑1 (MQ=255)                                                                                                                                          aaaaGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGAAGGGATgcg  <  1:531668/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                            AGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGA‑TGAAAAAGCTCTCCAACAATACCGCCTG‑C‑GTC‑TTC‑A‑CCGGT‑AAGGGCCTTTCATTTGGCGGCAG‑TCTTATTCGCCCGGAAGCTA‑C‑CGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATGGGTTTTGAAGGGATGCG  >  NC_000913/1842785‑1843071

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