Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,718	0	G	C	42.3%	53.5 / 28.3	55	P292R (CCG→CGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (10/20);  new base C (22/0);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 2.64e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTG  >  NC_000913/2984573‑2984868                                                                                                                                                  |                                                                                                                                                         attgacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑Ca                                                                                                                                                    <  1:327296/151‑1 (MQ=255)              cGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAATTCCCATCaa                                                                                                                                         >  1:508041/1‑149 (MQ=255)                aCGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCG‑CCATTCCCATc                                                                                                                                           >  2:821121/1‑145 (MQ=255)                    tCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAAATCCCATCAATGcc                                                                                                                                     >  1:636726/1‑147 (MQ=255)                      aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCt                                                                                                                               <  1:44872/151‑1 (MQ=255)                        cccACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCATCAATGCCCGGTCt                                                                                                                               >  1:393153/1‑149 (MQ=255)                          cACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCAACACCCGCG‑GCAATCCCATCAAATCCGGGTCt                                                                                                                               >  1:133442/1‑147 (MQ=255)                          cACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACAAACGCG‑CCATTCCCATAAATGCCCGATAT‑aca                                                                                                                           >  1:22412/1‑150 (MQ=255)                               aaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCG‑CAAATCCCCTCAATGCCGGATCT‑ACACCaa                                                                                                                       >  1:758337/1‑149 (MQ=255)                                aGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTCCTGCG‑CAATACGCCTCAATCCGAGATCT‑ACACCAACGc                                                                                                                    <  2:859313/151‑1 (MQ=255)                                  cATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTTCGGCG‑CAATTCGCAGAAATCCCGGATCT‑ACACCAACGCt                                                                                                                   <  2:107323/150‑1 (MQ=255)                                    tCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCCTCAATGCCGGATCT‑ACACCAACGCTgg                                                                                                                 >  1:240777/1‑150 (MQ=255)                                         ccAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATATTTTGTAA‑CGATACGCAGCCCTGCGAGATCT‑ACACCAACGCTGGCAAATg                                                                                                           <  2:70587/151‑1 (MQ=255)                                           aaCCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAGCGCTGGCAAATgcg                                                                                                         <  1:812382/151‑1 (MQ=255)                                             ccGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGcgcg                                                                                                       <  1:593703/151‑1 (MQ=255)                                              cGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACAACCGCG‑CCATTCCCATCAATGCCGGGTCT‑ACACCAACCCTGGCAAAAgcg                                                                                                         >  1:12853/1‑148 (MQ=255)                                                 aaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATg                                                                                                   <  2:55119/151‑1 (MQ=255)                                                  aGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CCATTCCCATCAATGCCCGATCT‑ACACCCACCCCGGCAAATGcgcgcaagc                                                                                                  >  1:101553/1‑146 (MQ=255)                                                  aGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGc                                                                                                  <  1:828451/151‑1 (MQ=255)                                                      ccTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCACCGGCGCCAATTCCC‑CCAACGCCGGAACC‑ACACCAACCCCCGCCAACGCGGGAAAGCgg                                                                                                >  2:476815/1‑149 (MQ=255)                                                           gcagcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTCTTTAA‑CGATTCCGAGCACTCCGAGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCa                                                                                         <  2:16118/151‑1 (MQ=255)                                                             agcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACg                                                                                       <  1:625742/151‑1 (MQ=255)                                                             agcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAAACCGCC‑CAAATCCCCTCAATCACGGGTCC‑ACAACCACGCCGGCAACAGGGCGAATGCGGGGCAg                                                                                           >  2:320738/1‑147 (MQ=255)                                                             agcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCGGCG‑CAATTTCCCTCAAATCCCGATCT‑ACACCCAACCCGGCAAAAGCGCGAAAGCggggc                                                                                             >  1:550722/1‑143 (MQ=255)                                                              gcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACgg                                                                                      <  1:170736/151‑1 (MQ=255)                                                                aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTTCCATCAATGCCGGATCT‑ACAACAACGCCGGCAAAAGCGCGAAAGCGGGCCCGCACgggg                                                                                    >  2:923581/1‑150 (MQ=255)                                                                aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGCCGGAAAAACCCGGC‑CAATTTCCATCAATGCCCGATCT‑ACACCAAAACTGGCGAAAACGCGGAGGCGGGTCCGCCCgggg                                                                                    >  1:450780/1‑150 (MQ=255)                                                                       tCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CCATTTCCCTCAATGCCGGCTCT‑TCACCAACCCCGGGAAAAGCGCCGATGCGGGCCAGCCCGGGGaaa                                                                                 >  1:3134/1‑147 (MQ=255)                                                                       tCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGATCCCATCCACCAGCG‑CCATTCCCACCCATGCCGGGTCT‑TCACCAACGCCGCGCAACGGGCGAAAGCGCGCCCGCACgggggaaat                                                                               >  2:34526/1‑143 (MQ=255)                                                                           tAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAAATCCCCTCCATGCCCCATCT‑ACACCCACCCTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAATTCAAcac                                                                         >  2:202989/1‑149 (MQ=255)                                                                            aGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCG‑CAAATTCCCCCCATGCCGGACCT‑ACACCCACCCCGGCAAAAGCCCGAAAGCGGGGCAGCCCGGGGAAAATCAAAGcc                                                                        >  2:555806/1‑150 (MQ=255)                                                                                tCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCTCATACACCCGCG‑CCAATCCCCCCAATTCCCGAGCT‑ACACCCACGCTGGCAAATGCGCGGATGCGGGCCCGCCCGGGGAAATTTAAAGCCCGt                                                                     >  2:796729/1‑150 (MQ=255)                                                                                  gATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCAACAACCGGG‑CAATTCCCATCAATGCCGGGTCT‑ACACCCACCCTGGCAAATACGCGAAAGCGGGGCAGCACCGGTAAAATTAACGccc                                                                       >  1:567612/1‑144 (MQ=255)                                                                                   aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCtt                                                                 <  2:697926/151‑1 (MQ=255)                                                                                       gATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCATCAATACCGGATCT‑ACACCAACGCCGGCAAAAGCGCGAATGCGGGCCAGCACCGGGAAATTCAACGCCCGTGCttttg                                                              >  1:462669/1‑148 (MQ=255)                                                                                         tCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGATCGCATAAACCCCCG‑CAATTTCCACCAATGGCGGGTCT‑AAAACAACGCTGGGAAATTCGACAATGACGGCCTCCCCGGGGAAATTCAACACTcccgcctttgcct                                                           >  1:623078/1‑139 (MQ=255)                                                                                          cGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTc                                                          <  1:197768/151‑1 (MQ=255)                                                                                           gACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCg                                                         <  1:181519/151‑1 (MQ=255)                                                                                            aCTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAATTTCCATCAATGCCGGATCT‑ACACCAACCCCGGCAAAAGCGCGCATGCGGGGCCGCACGGGGAAAATCAACCCTCCGGCttttg                                                              >  2:1039654/1‑143 (MQ=255)                                                                                                 aGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACAACGGCG‑CAAATCCCCTCAATGCCCGCGCT‑ACACCAACGCCGGCCAAAGCGCGAATGCGGGCCAGCACGGGGAAAATCAACGCCCGTGCTTTGGCTTCGc                                                        >  2:196465/1‑146 (MQ=255)                                                                                                   ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTGGCATACACCCGCG‑CACTTCCCCCCCATACCGGGTCTCACCCAAACCCT‑GCAAATGCGCCGCTGCCGGCCAGCACGGGAAGAATCAAACCTCGTGCTTTGGCTTctctactcc                                                  >  2:899324/1‑145 (MQ=255)                                                                                                      aCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTcatcatcat                                              <  1:624046/151‑1 (MQ=255)                                                                                                      aCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTcatcatcat                                              <  2:644444/151‑1 (MQ=255)                                                                                                       cTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCCTCAATGCCGGATCT‑ACACCAACGCTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAATTCAACGCTCGTGCTTTTGCTTCGCTCATCATCAtt                                             >  2:119471/1‑151 (MQ=255)                                                                                                           cAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCg                                         <  1:768084/151‑1 (MQ=255)                                                                                                              ccTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAATTCCCCTCAATGCCGGATCT‑ACACCAACGCTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAATTCAACGCTCGTGCTTTTGCTTCGCTCATCATCCTTACCGc                                        >  2:989695/1‑149 (MQ=255)                                                                                                                         cAGGCAACGGCGGGTCGCATTCACCGGCG‑CAATTCCCATCAATCCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATc                            <  2:983789/150‑1 (MQ=255)                                                                                                                               aCGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATg                     <  1:56340/151‑1 (MQ=255)                                                                                                                                  gcggGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgat                  <  1:368104/151‑1 (MQ=255)                                                                                                                                      gTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGatgatg              >  2:596797/1‑151 (MQ=255)                                                                                                                                       tCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGc             >  2:519763/1‑151 (MQ=255)                                                                                                                                        cGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGc             <  2:393046/150‑1 (MQ=255)                                                                                                                                           aTACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt          >  1:94352/1‑150 (MQ=255)                                                                                                                                           aTACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt          >  1:531490/1‑150 (MQ=255)                                                                                                                                           aTACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt          >  2:242046/1‑150 (MQ=255)                                                                                                                                                  ggCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTg  <  1:777436/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                         ATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCTG  >  NC_000913/2984573‑2984868

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