Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,842,914	0	A	G	45.5%	40.3 / 17.3	46	M182V (ATG→GTG)	gdhA	glutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base A (0/24);  new base G (20/0);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 5.68e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATG  >  NC_000913/1842764‑1843054                                                                                                                                                       |                                                                                                                                                  gATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGga                                                                                                                                                   <  2:378628/151‑1 (MQ=255)      cGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGgatgat                                                                                                                                               <  1:86124/150‑1 (MQ=255)        aTCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGgggat‑g                                                                                                                                             >  1:101942/1‑149 (MQ=255)         tCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGgggat‑gga                                                                                                                                           >  1:374280/1‑148 (MQ=255)           cGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑Gaaaaa                                                                                                                                        <  2:664488/151‑1 (MQ=255)                     aaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGG‑GAAAAAACACT‑ccca                                                                                                                              >  2:542417/1‑148 (MQ=255)                     aaGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑Ccaac                                                                                                                             <  2:399537/151‑1 (MQ=255)                      aGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGAGGTGGT‑aaaaaa                                                                                                                                        >  1:718373/1‑140 (MQ=255)                      aGCGAAAGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTTGGCGCGGATACCGACGTTCCGGGAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑Ccaac                                                                                                                             <  2:371022/150‑1 (MQ=255)                       gCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGT‑Gaaaaaaaactcc                                                                                                                                 >  2:140305/1‑139 (MQ=255)                       gCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑Ccaacaa                                                                                                                           <  1:418256/151‑1 (MQ=255)                        cGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAAt                                                                                                                          <  2:99247/151‑1 (MQ=255)                          aaGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGAT‑GGAAAAAATCT‑Cccacaaaa                                                                                                                         >  1:633441/1‑148 (MQ=255)                              tGAAGTGATGCGTTTTTGCCAGACGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGT‑GGAAAAAATCT‑Cccacaaaaacccc                                                                                                                    >  2:750453/1‑144 (MQ=255)                               gAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCt                                                                                                                   <  2:265275/151‑1 (MQ=255)                                   tGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGCGTT‑GGAAAAGGTTT‑TCCACAATACAGc                                                                                                                     >  2:102118/1‑145 (MQ=255)                                       gCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGCTGAT‑GCAAAAGGACT‑CCCACCAAAACGCCTGccgctt                                                                                                            >  1:202951/1‑145 (MQ=255)                                        cGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGT‑GGAAAAGGACG‑CCCACAAAAACCCCCGCGTc                                                                                                              >  2:721626/1‑147 (MQ=255)                                         gTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGT‑GGAAAAACACTCCCAACAA‑ACCCCCTGCGTCTTccc                                                                                                         >  2:720235/1‑149 (MQ=255)                                                ccAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGTGGAAAAAATCCT‑CCACCAA‑AACGCCTGCGTGTTCA‑c                                                                                                        >  1:286056/1‑144 (MQ=255)                                                  aGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGT‑GAAAAAAATCT‑CCCACAAAACCGCCCCCCTCTTTA‑ccccgaa                                                                                                  >  2:556697/1‑143 (MQ=255)                                                    gcgcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGcc                                                                                           <  2:397010/151‑1 (MQ=255)                                                     cgcTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGAT‑GAAAAAAAACT‑CCCACAAAAACGCCCGCGTCTTCA‑cccggaaagg                                                                                               >  1:26011/1‑140 (MQ=255)                                                            gaCTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGT‑GAAAAAACACT‑CCCACAAAAACCCCTGCGGCTTCACCC‑GGAAA‑GGGCCttccttt                                                                                    >  1:520814/1‑146 (MQ=255)                                                                gAACTGTATCGCCCCCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTggc                                                                                <  1:238441/150‑1 (MQ=255)                                                                  aCTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTggcggc                                                                             <  2:175996/151‑1 (MQ=255)                                                                     gTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGT‑TGAAAAAATAT‑CCTCCCATACCGCCTGCGTCTTCT‑TC‑GGCAA‑AGGGCATTCtttttgtggggg                                                                           >  1:102857/1‑142 (MQ=255)                                                                      tATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAg                                                                           <  1:781854/149‑1 (MQ=255)                                                                               cTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGGT‑GGAAAAAAACT‑CCCACAAAAACGCCCGCGGCTTTC‑CC‑CGGAA‑GGGGCTTTTAtttttggggaggcttat                                                                     >  2:55555/1‑132 (MQ=255)                                                                                 gggCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCg                                                              <  2:701401/151‑1 (MQ=255)                                                                                     gcggtcaCCGACGTCCCGGCAGGTGATTTCGGGGTTGGTGGTCGGGAAGTCGGCTTTATGGCGGGGTTGAT‑GAAAAAGCTCT‑CAAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAg                                                          <  2:40299/145‑1 (MQ=255)                                                                                          tACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGT‑GAAAAAAATCT‑CCCACAAAAACGCCTGCGGCTTTA‑CC‑GGGAA‑GGGGCCTTTATTTGGCGGCGGGCTTTTTTGCCCGGGAGCTAcc                                                     >  1:223881/1‑151 (MQ=255)                                                                                              gACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGGGAT‑GAAAAAACGCC‑CCCACCAAAACCCCTCCCTCTTTT‑CCGGGAAA‑GGGCCTTTTTTTTGG‑GGGAGGCTTTTTTGCCCCGGAGCTccc                                                     >  2:469583/1‑147 (MQ=255)                                                                                                cGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTAc                                               <  1:254122/151‑1 (MQ=255)                                                                                                      ggCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTg                                         <  1:743612/151‑1 (MQ=255)                                                                                                      ggCAGGTGATACCGGGGTTGGTGGTCGTGAAGTCGGCTTTTTGGCGGGGATGAG‑AAAAAAGCTCC‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTg                                         <  2:342535/151‑1 (MQ=255)                                                                                                       gCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGGT‑GAAAAAAGTCT‑CCCACAATACCCCCCGCGGCTTTA‑CC‑CGGAA‑GGGGCTTTCTTTTTGGGGGAGGCTTTTTTTCCCCGGAGCCAcccggcccgggctgg                                        >  1:95262/1‑138 (MQ=255)                                                                                                                   ggggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTcaca                             <  1:408559/150‑1 (MQ=255)                                                                                                                   ggggTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAg                            <  2:81226/151‑1 (MQ=255)                                                                                                                       ttGGTGGTCGTGAAGTCGGCTTTATGGCGGGGGTGAT‑GAAAAAAATCT‑CCCACAAAACCGCCCGCGTCTTCA‑CC‑CGGAA‑GGGCCTTTTCTTTTGCGGGAGGCTTATTTTCCCGGGAGCTACCCGCTACGGGCTGGGTTTTTTTACAGaa                          >  1:759072/1‑149 (MQ=255)                                                                                                                         ggtggtCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCa                       <  2:101941/150‑1 (MQ=255)                                                                                                                            ggtCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGc                   <  2:216849/151‑1 (MQ=255)                                                                                                                              tCGTGAAGTCGGCTTTATGGCGGGGGTGAT‑GAAAAAACACT‑CCCACAAAAACGCCCGCGGCTTTA‑CC‑GGGAA‑GGGGCTTTCATTTTGGGGGAGGCTTATTTTCCCCGGAGCCACCCGCTCCGGGCTGGTTTTTTTTACAGAAACAAAGc                   >  2:752975/1‑149 (MQ=255)                                                                                                                               cGTGAAGTCGGCTTTATGGCGGGGGTGAT‑GAAAAAACTCT‑CCCACAAAAACCCCCCCGTCTTTA‑CC‑GGGAAGGGGCCTTTTTTTTGG‑GGGAGTCTTTTTTTCCCCGGAGGCACCCGCCACGGGCTGGGTTTTTTTTCAGAAGCAAAGCt                  >  1:693633/1‑149 (MQ=255)                                                                                                                                gTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTaaa               <  1:648227/151‑1 (MQ=255)                                                                                                                                   aatccGGCTTTATGGCGGGGATGAT‑GAAAAGGCGTT‑CCAACAACACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACg            <  2:304387/147‑1 (MQ=255)                                                                                                                                   aaGTCGGCTTTATGGCGGGGATGAT‑GAAAAATCTCT‑CCAACAATGCCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAAc             <  2:23695/150‑1 (MQ=255)                                                                                                                                      tCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGc           <  1:129997/149‑1 (MQ=255)                                                                                                                                           tttggggCGGGGATGAT‑CAAAAAGCCCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACgg      <  2:428918/144‑1 (MQ=255)                                                                                                                                               tGGCGGGGATGAT‑GAAAAAGCACT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATg  <  2:136738/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  GATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGAT‑GAAAAAGCTCT‑CCAACAATACCGCCTGCGTCTTCA‑CC‑GGTAA‑GGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGGTCTGGTTTATTTCACAGAAGCAATGCTAAAACGCCACGGTATG  >  NC_000913/1842764‑1843054

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