Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,264,338:1	+T	12.8%	intergenic (‑44/+40)	zinU ← / ← PP_4640	P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,264,338	1	.	T	12.8%	161.6 / 20.5	47	intergenic (‑44/+40)	zinU/PP_4640	P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (20/21);  new base T (4/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

CGATCACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGT  >  NC_002947/5264189‑5264486                                                                                                                                                       |                                                                                                                                                     cgctCACGGCGAGCTTAAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGcc                                                                                                                                                       <  1:1059646/147‑1 (MQ=255)  taTCACGGCGAGCTTCAGGCCAGGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑g                                                                                                                                                     <  1:521590/149‑1 (MQ=255)  gATCACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑g                                                                                                                                                     <  2:99485/150‑1 (MQ=255)  gATCACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑g                                                                                                                                                     >  1:663066/1‑150 (MQ=255)  gATCACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑g                                                                                                                                                     <  1:227125/150‑1 (MQ=255)    tCACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGGAGCCTGTGAGGACGGTTACGGGAATGGGCCTCTGCAAGCAAAATACTCCTGAATCGGAAAATGACTGTTGGAGCC‑ggg                                                                                                                                                   >  2:521590/1‑149 (MQ=255)       cGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTAc                                                                                                                                                <  1:491205/150‑1 (MQ=255)         gCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTAccc                                                                                                                                              >  2:1119417/1‑150 (MQ=255)              cTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCt                                                                                                                                         <  2:663066/150‑1 (MQ=255)               ttCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTc                                                                                                                                        >  2:439461/1‑150 (MQ=255)                     ccATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAg                                                                                                                                  <  1:765167/150‑1 (MQ=255)                     ccATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGAGGCCACAg                                                                                                                                  >  1:60290/1‑150 (MQ=255)                      cATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCGTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAgg                                                                                                                                 >  2:1174534/1‑150 (MQ=255)                                   gAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCg                                                                                                                    <  1:439461/150‑1 (MQ=255)                                      cATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAg                                                                                                                 <  2:2863/150‑1 (MQ=255)                                       aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAAGACTCCTGAATCTGAAAATGACTCTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAgg                                                                                                                <  2:640365/150‑1 (MQ=255)                                          tACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTc                                                                                                             >  2:1267096/1‑150 (MQ=255)                                                  cAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCAc                                                                                                     >  1:140796/1‑150 (MQ=255)                                                     ggtggtCTTGCCGGCACCGAGGAAGCCTGTAAATCTACACACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTg                                                                                                  <  2:1181912/150‑1 (MQ=255)                                                          tCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTGCTGAATCTGCAAGTGACTGTGGGAGCC‑GGTTACCCAGGCCCCCACAGGTTTACACGAACCGCAGGCCGAGCCGCACCTGCGccc                                                                                             >  2:48433/1‑150 (MQ=255)                                                            ttGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCgg                                                                                           >  2:618881/1‑150 (MQ=255)                                                                 ggCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTg                                                                                      >  2:796842/1‑150 (MQ=255)                                                                    aCCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCa                                                                                   <  1:1267096/150‑1 (MQ=255)                                                                            aGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACg                                                                           <  2:1330853/150‑1 (MQ=255)                                                                             gCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCCACGAGGTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGt                                                                          >  1:375038/1‑150 (MQ=255)                                                                                   gAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTcc                                                                    <  2:688630/150‑1 (MQ=255)                                                                                      gACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACC‑GGCTCCcaca                                                                                                                                   >  2:964784/1‑84 (MQ=25)                                                                                      gACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACC‑GGCTCCcaca                                                                                                                                   <  1:964784/84‑1 (MQ=25)                                                                                           ttACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTc                                                            <  1:796842/150‑1 (MQ=255)                                                                                           ttACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGAACTCGCGTCCGTAACGTGCTTCCTGACGttt                                                            >  1:714704/1‑149 (MQ=255)                                                                                              cGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGttagga                                                         <  1:1119417/150‑5 (MQ=255)                                                                                               gggAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCgg                                                        >  2:1335337/1‑150 (MQ=255)                                                                                                       gCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GTGTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTTCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAAttc                                                >  2:200530/1‑148 (MQ=255)                                                                                                              cacagaaaaTACGCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAg                                         <  1:48433/145‑1 (MQ=255)                                                                                                              cACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAg                                         <  1:1174534/150‑1 (MQ=255)                                                                                                              cACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAg                                         <  1:823991/150‑1 (MQ=255)                                                                                                                cGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAg                                                                                                        >  1:610177/1‑85 (MQ=255)                                                                                                                cGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAg                                                                                                        <  2:610177/85‑1 (MQ=255)                                                                                                                   aaaaTACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAg                                                                                                                                  <  2:1380525/56‑1 (MQ=255)                                                                                                                   aaaaTACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAg                                                                                                                                  >  1:1380525/1‑56 (MQ=255)                                                                                                                      aTACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCAc                                 <  2:60290/150‑1 (MQ=255)                                                                                                                       tactacGGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCAcc                                >  1:1058984/1‑150 (MQ=255)                                                                                                                        gcTCCTGAATCGGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCg                               <  2:454302/149‑1 (MQ=255)                                                                                                                              gAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGt                         >  2:869762/1‑150 (MQ=255)                                                                                                                                 tCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCCCATCACCGCCTTGTCCg                      >  2:227101/1‑150 (MQ=255)                                                                                                                                    gAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGAt                   <  1:618881/150‑1 (MQ=255)                                                                                                                                     tgtaaaccTGTGGGAGCCTGGGTAACC‑GGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATg                  >  2:842008/8‑150 (MQ=38)                                                                                                                                            cTGTGGGAGCCTGGGTAACC‑GGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                 <  2:700613/64‑1 (MQ=18)                                                                                                                                            cTGTGGGAGCCTGGGTAACC‑GGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTgg                                                                                                 >  1:700613/1‑64 (MQ=18)                                                                                                                                            cTGTGGGAGCCTGGGTAACC‑GGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTc           >  1:821490/1‑150 (MQ=255)                                                                                                                                             tgtgGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCt          >  1:96869/1‑150 (MQ=255)                                                                                                                                                 ggAGCC‑GGTTACCCAGGCTCCCACAGGCATGTACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt      >  1:37787/1‑150 (MQ=255)                                                                                                                                                   aGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAtgt    <  2:140796/150‑1 (MQ=255)                                                                                                                                                     cc‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATgtgt  >  1:472887/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                     CGATCACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGT  >  NC_002947/5264189‑5264486

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