Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	41.5%	81.3 / 15.7	53	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (26/5);  new base C (0/22);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 2.07e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTG  >  NC_000913/907454‑907749                                                                                                                                                         |                                                                                                                                                    gTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTAcc                                                                                                                                                   <  1:397754/151‑1 (MQ=255)     cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTc                                                                                                                                               <  2:650247/151‑1 (MQ=255)     cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTc                                                                                                                                               >  1:638383/1‑151 (MQ=255)     cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTc                                                                                                                                               >  1:145528/1‑151 (MQ=255)            tcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGcc                                                                                                                                         >  1:876896/1‑150 (MQ=255)                          ttggtttttcgcccctTGAATCACCGCCCGGGCGGTGGGGGGGGGTAGGGAGGTAATCAAATTGAAAACCCCCCCCA‑AAACCC‑A‑GGGGGTGGGGGGGGGGAGGGGGATCACCCGGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCaaa                                                                                                                           <  1:315144/138‑1 (MQ=255)                              tgtTTCCGCCCCTGAATCACCGCCCTGGGGGTGGGGGGGGGTACGGAAGTAATCACATTGAGAAACCACACCA‑AAACCC‑A‑GGGGGGGGGGGGGGGGAGGGGGATCAACCAGCGG‑CGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGt                                                                                                                       <  1:831162/150‑1 (MQ=255)                                   cAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                                  >  1:375002/1‑150 (MQ=255)                                      ccccTGAACCACCCCCCGGGGGGGGGGGGGGGGTAGGAAGGTACGCACATGGAGAACCCACCCCA‑AACCCC‑A‑GCGGGGGGGGAGGGGGAGGGTGATCAGCCGCCGG‑CCCTTCCCACCCCCCGCCATAAAATGTTCAAACTGTTGCTGAt                                                                                                                <  2:279881/149‑1 (MQ=255)                                        ccgaaaccAACGCCCTGGGGGTGGTGGGGGGTGCGGAAGTAATCCCAGTGAGAACCTACCCCG‑AACCCC‑A‑GCGGGTGGGGGGGGGGATGGGAATCAACCGGCGG‑CGCTTACCCCACACCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                                            <  2:236625/144‑1 (MQ=255)                                        ccTAAATCACCGCCGTGGGGGTGGTGGGGGGTAGGGAAGTAATCACATTGAAAACCCACCCCA‑AAACCC‑G‑GGGGGTGGGGGGGGGGAGGGGGATAAACCAGCGG‑CGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                                            <  1:634474/151‑1 (MQ=255)                                            aaacAACGCCCGGGGGGTGGGGGGGGGTACGGAGGTAATCACATTGAGAAACCACCCCA‑AAACCC‑G‑GGGGGGGGGGGGGGGGATGGGGATCACCCGGGGG‑CTCTTCCCGCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGc                                                                                                        <  2:944170/148‑1 (MQ=255)                                            aaTCACCCCCCGGGGGGTGGGGGGCGGTACGGAATTAATCACATTGAAAACCCCCCCCA‑AAACCC‑G‑GGGGGTGGGGGGGGGGAGGGGGATCAACCAGGGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAg                                                                                                         <  1:492116/150‑1 (MQ=255)                                              tCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAc                                                                                                      >  1:151985/1‑151 (MQ=255)                                                accccccTGGGGGGGGGGGGGGGTACGGAAGTAATCACATTGAAAACCCACCCCA‑AAACCC‑A‑GGGGGGGGGGGGGGGGAGGGGAACCACCCAGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCa                                                                                                       <  1:296112/144‑1 (MQ=255)                                                aaCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACct                                                                                                    >  1:783340/1‑151 (MQ=255)                                                  cGCCGTGGCGGTGGGGGGAGGGACGGAGCTAATCATTGTGAAAATCCCCCCGG‑AACCCG‑C‑GCGGGTGGGGAGGGGGATGGGTAGC‑ACCAACGGACGCTTCCCCCACCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGc                                                                                                        <  2:708785/145‑1 (MQ=255)                                                  cGCCCGGGGGGGGGTGGGGGGCACGAAAGTAATCACATTGACAACCCACACCA‑AATCCC‑C‑GGGGGTGGGGAGGGGGAGGGGGAACAGCCAGCGC‑CGCCAACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctc                                                                                                   <  2:478420/150‑1 (MQ=255)                                                     ccTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTgg                                                                                                >  2:291003/1‑150 (MQ=255)                                                        gggggtggtGGGCGGTACGAAATTAACAAAATGGAGAACCCCCCCCG‑AAACCC‑A‑GCGGGGGGGGGGGGGGAGGGGGACCACCCAGCGC‑CGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                            <  2:786838/148‑1 (MQ=255)                                                          cGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggc                                                                                          >  1:442112/1‑151 (MQ=255)                                                          cGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGGTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcgg                                                                                           >  2:192087/1‑150 (MQ=255)                                                           ggtggtGGGCGGTAGGAAAGAAATCACATTGAGAACCCACCCCA‑AACCCC‑G‑GGGGGGGGGGGGGGGGAGGGGAATCAACCAGCGG‑CGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                         <  1:182613/151‑1 (MQ=255)                                                            gggggggggggTACGGAAGTAATCAAATTGAAAAACCCCCCCA‑AAACCC‑G‑GGGGGGGGGGGGGGGGATGGGAATCAACCGGCGG‑CGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                        <  1:143846/146‑1 (MQ=255)                                                              ggtggGCGGTACGGAAGTAATCACATGGAAAAACCCCCCCA‑AACCCC‑G‑GGGGGGGGGGGGGGGGATGGGGATCAACCAGCGG‑CGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                         <  1:104113/148‑1 (MQ=255)                                                              ggggggggtaaCGGAAGTAATGACAGAGAAAACCCAC‑CCA‑GAATCC‑AGGCGGTTGGGGGGGGGGAGGGGGACGAGCCAGCGC‑CGCTACCCCCCCGCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTc                                                                                              <  2:739819/133‑1 (MQ=255)                                                                  ggCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCg                                                                                   >  2:659195/1‑150 (MQ=255)                                                                     ggTAGGGAAGTAATCACATGGAGAACCCACACCA‑AAACCC‑G‑GGGGGTGGGGGGGGGGAGGGGGATCAACCGGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAg                                                                                <  1:327400/150‑1 (MQ=255)                                                                     ggTACGGAAGAAACCATATGGAAAACCCCCCCCA‑AACCCC‑G‑GGGGGGGGGGGGGGGGAGGGGAATCAACCGGGGG‑CTCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAg                                                                                <  1:305878/150‑1 (MQ=255)                                                                         cGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                                           >  1:358737/1‑151 (MQ=255)                                                                          ggAACTAATCATATTAAAAACCCACCCCA‑AACCCC‑G‑GGGGGTGGGGGGGGGGAGGGGAACCAACCAGGGG‑CCTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                                           <  1:695511/150‑1 (MQ=255)                                                                           gAATTATTCACATTGAGAACCCAC‑CCG‑AACCCCAA‑GGGGGGGGGGAGGGGTATGGGAATAACCCGGGGC‑CGCTTCCCACTCACGGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGag                                                                              <  2:111970/146‑2 (MQ=255)                                                                           gAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                          >  2:663745/1‑150 (MQ=255)                                                                            aattAATCAAATGGAGAACCCACACCA‑AAACCC‑G‑GGGGGGGGGGGGGGGGAGGGGAATCACCCAGCGG‑CGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                          <  1:121655/146‑1 (MQ=255)                                                                                aaTCACATTAAAAACCCACCCCG‑AAACCC‑G‑GGGGGGGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGtt                                                                      <  1:758199/149‑1 (MQ=255)                                                                                aaTCACA‑TGAAAAACCCCCCCGAAAACCC‑G‑GGGGGGGGGGGGGGGGGTGGGAATCACCCACGGG‑CCCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCg                                                                    <  2:373357/151‑1 (MQ=255)                                                                                    acaTTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCCTCGAAGATGCGGTGTTCGCTgg                                                                >  2:128356/1‑151 (MQ=255)                                                                                     caTTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGc                                                               >  2:16538/1‑151 (MQ=255)                                                                                      aTTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                                              >  1:299113/1‑151 (MQ=255)                                                                                          aaacaccaccCCG‑AACCCC‑G‑GGGGGTGGGGGGGGGGTTGGGGATCACCGAGCGG‑CGCTCCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                          <  2:187172/147‑1 (MQ=255)                                                                                          aaaacccACACCA‑AAACCC‑A‑GGGGGGGGGGAGGGGGAGGGGGATCAACCGGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                                             <  1:814447/143‑1 (MQ=255)                                                                                               ccACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                     >  1:140953/1‑151 (MQ=255)                                                                                                            cc‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                         >  1:704863/1‑151 (MQ=255)                                                                                                                 gcgggcggGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTACCAACGCAAcc                                         >  2:311014/1‑148 (MQ=255)                                                                                                                   gggTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                     >  1:444043/1‑150 (MQ=255)                                                                                                                           cgggggATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc                             >  2:154752/2‑150 (MQ=255)                                                                                                                                    gggatAACCCAGCGG‑CGCTTCCCATTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTTGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCa                   <  1:774296/151‑1 (MQ=255)                                                                                                                                      gatCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgc                 >  1:442527/1‑151 (MQ=255)                                                                                                                                              aGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg         >  1:873594/1‑151 (MQ=255)                                                                                                                                               gcgg‑cgCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGc        >  1:125418/1‑151 (MQ=255)                                                                                                                                                  g‑cgcTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtg     >  1:646141/1‑151 (MQ=255)                                                                                                                                                    cgcTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtgg    >  1:924544/1‑151 (MQ=255)                                                                                                                                                      cTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtggtg  >  2:700993/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                    GTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCG‑AAATCC‑A‑GCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTG  >  NC_000913/907454‑907749

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