Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	46.9%	62.2 / 13.0	49	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (20/6);  new base C (0/23);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.63e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCGT  >  NC_000913/907460‑907752                                                                                                                                                |                                                                                                                                                       cgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCAc                                                                                                                                                <  2:293056/151‑1 (MQ=255)  ggTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCAc                                                                                                                                                <  1:183210/150‑1 (MQ=255)   gTATCATCCGCCAGAAGAGCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGCAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCg                                                                                                                                              <  1:90679/151‑1 (MQ=255)                          tttCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTCCGGAAGTAATCACATTGAAAAACCACCCCGAAACCCAGCGGGTGGGGAGGGGGAGGGGGATCACCCAGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGtt                                                                                                                         <  2:662901/149‑1 (MQ=255)                           ttCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCt                                                                                                                      >  2:849276/1‑151 (MQ=255)                            tCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                                     >  1:866739/1‑151 (MQ=255)                                cccagaagtCAACCCCCTGGCGGTGGTGGGCGGTATGGAAGAAATCACATTGAAAAACCACGCCAAAGTCCGGGGGGTGGGGGGGGGAATGGGTATAAACCAGCGG‑TTCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATa                                                                                                                  <  2:229673/142‑1 (MQ=255)                                    tGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                              >  2:95941/1‑150 (MQ=255)                                         cagcGCCCGGGCGGGGGGGGGCGGTAGGAAAGAAATCAAATGGAAAACCCACCCCAAAACCCGGCGGGGGGGGAGGGGGTTGGGAATCACCCGGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCa                                                                                                          <  1:502726/146‑1 (MQ=255)                                         cAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                        >  1:89320/1‑151 (MQ=255)                                          aaCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACct                                                                                                       >  2:472382/1‑151 (MQ=255)                                          aaCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACct                                                                                                       >  1:492468/1‑151 (MQ=255)                                                  ggCGGTGGGGGGCGCTCCGGAATTAATCACATGAAGAACCCACCCCAAAATCCGGGGGGTGGGGGGGGGGAGGGGGACCAACCAGCGG‑CTTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                               <  2:25204/151‑1 (MQ=255)                                                      gtggggGGGGGTAGGGAAGAAACCACATTGAGAACCCCCACCAAAATCCGGGGGGGGGGGGGGGGGAGGGGGATCACCCAGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                           <  2:591794/151‑1 (MQ=255)                                                           ggggggTACGGAAGTAATCACATGGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCAGGGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                                       <  2:345944/150‑1 (MQ=255)                                                            ggCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCg                                                                                      >  2:3304/1‑150 (MQ=255)                                                                 tACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATg                                                                                >  1:326880/1‑151 (MQ=255)                                                                 tACGGAAGTAATCACATTGAGAAACCACACCGAAACCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGAt                                                                                 >  2:136813/1‑150 (MQ=255)                                                                    ggAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                             >  1:680575/1‑151 (MQ=255)                                                                    ggAAGAAATCACATTGAAAAACCACCCCAAAACCCGGGGGGGGGGGGGGGGGGGGGGGATAAACCAGCGG‑CGTTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGAt                                                                                 <  2:597997/147‑1 (MQ=255)                                                                    gaaaGTAACCCCATGGAAAAACCCCCCCAAAACCCGGGGGGGGGGGAGGGGGAGGGGGATCAACCGGCGC‑CGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                             <  2:878524/149‑1 (MQ=255)                                                                      aaaaagtCACATGGAGAACCCACACTGAAATCGGGGGGGGGGGGAGGGGGGGGGGGATCACCCAGCGG‑CTTTTCCGCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgt                                                                            <  2:3742/144‑1 (MQ=255)                                                                      aaGTAATCACATAGAAAAACCACCCCAAAACCCGGGGGGGGGGGAGGGGGAGGGGGATCAACCAGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                           <  2:812626/151‑1 (MQ=255)                                                                      aaGTAACCACATGAAAAAACCACACCAAAACCCGGGGGGGGGGGGGGGGGAGGGGAACCAACCAGGGG‑CCCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                                              <  1:807363/148‑1 (MQ=255)                                                                          atccacaTTAAGAAACCCCCCCAAAACCCGGGGGGGGGGGGGGGGGGGGGGGAACACCAAGCGG‑CGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGc                                                                               <  2:252900/140‑1 (MQ=255)                                                                                  ggagaAACCACACCAAAATCCGGCGGGGGGGGGGGGGGAGGGGGATCACCCAGCGG‑CGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                               <  1:107651/150‑1 (MQ=255)                                                                                   gagaAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                               >  1:280354/1‑150 (MQ=255)                                                                                    aaaacccACCCCAAAACCCGGCGGGTGGGGGGGGGGATGGGTACCACCCAGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                             <  2:882672/146‑1 (MQ=255)                                                                                       acccACCCCGAACTCCGGGGGGTGGGGGGGGGGAGGGGGACCAACCAGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                                          <  2:365683/149‑1 (MQ=255)                                                                                         ccACCCCAAAACCCAGGGGGTGGGGAGGGGGATGGGGATAACCCAGCGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                        <  2:415636/151‑1 (MQ=255)                                                                                            ccccGAAACCCGGGGGGGGGGGAGGGGGAGGGGGATCACCCGGCGC‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTa                                                       <  1:393162/147‑1 (MQ=255)                                                                                              ccgaaggccgggggggggggggggggAGGGGAAGCAACCAGGGG‑TGTTTCCCGTTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaa                                                   <  2:737084/139‑1 (MQ=255)                                                                                              ccaaaacccgggggggggggAGGGGGAGGGGGAACACCCAGGGG‑CGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaaca                                                    <  2:382770/138‑1 (MQ=255)                                                                                                 aaaTCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                                                >  2:350145/1‑151 (MQ=255)                                                                                                  accccAGGGGGGGGGGGGGGGGAGGGGGATCAACCAGCGG‑CGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                               <  2:357154/148‑1 (MQ=255)                                                                                                    tCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAc                                             >  2:571219/1‑151 (MQ=255)                                                                                                     ccAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                            >  2:272011/1‑151 (MQ=255)                                                                                                          ggggggggAGGGGGAGGGGACCCAACCAGCGG‑CGCTCCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATa                                       <  1:750372/151‑1 (MQ=255)                                                                                                                      ggatggggatCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                           >  1:95820/1‑151 (MQ=255)                                                                                                                       gatggggatCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                          >  1:845995/1‑151 (MQ=255)                                                                                                                          gggaatCAACCGGCGC‑CCCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                          <  2:649838/148‑1 (MQ=255)                                                                                                                          gggaaccacccAGGGG‑CCCTTCCCGCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCCGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTc                        <  1:295327/141‑1 (MQ=255)                                                                                                                           ggggctcaCCCGCGGCCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                           <  2:508841/140‑1 (MQ=255)                                                                                                                              aacAACCGGCGG‑CGCTTACCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAg                     <  1:677815/147‑1 (MQ=255)                                                                                                                                 aaCCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa                 >  2:283433/1‑150 (MQ=255)                                                                                                                                 aaCCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa                 >  1:620819/1‑150 (MQ=255)                                                                                                                                      gcgg‑cgCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATgg             <  2:287755/149‑1 (MQ=255)                                                                                                                                             cTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtggtg     >  2:464971/1‑151 (MQ=255)                                                                                                                                              ttACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGc    >  1:135397/1‑151 (MQ=255)                                                                                                                                                aCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCGt  >  1:95654/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                       CGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAGCGG‑CGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCGT  >  NC_000913/907460‑907752

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