Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,491,144	G→C	7.4%	S267T (AGT→ACT)	cysNC →	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,491,144	0	G	C	7.4%	166.5 / 4.9	54	S267T (AGT→ACT)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (22/28);  new base C (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

TCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAG  >  NC_002947/1490996‑1491283                                                                                                                                                     |                                                                                                                                             tCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGt                                                                                                                                             <  2:198983/150‑1 (MQ=255)   gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTgg                                                                                                                                           >  1:689323/1‑150 (MQ=255)    ccGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGc                                                                                                                                          >  1:1145469/1‑150 (MQ=255)       tGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tg                                                                                                                                      <  1:1411995/150‑1 (MQ=255)       tGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tg                                                                                                                                      <  1:1089124/150‑1 (MQ=255)         gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGGCgatggt                                                                                                                                    <  2:1376719/150‑6 (MQ=255)             gCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGcac                                                                                                                                >  1:1417593/1‑150 (MQ=255)                 gCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAg                                                                                                                            <  2:1145469/150‑1 (MQ=255)                   ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGCCCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTg                                                                                                                                            >  1:531215/1‑133 (MQ=255)                   ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTg                                                                                                                                            <  2:531215/133‑1 (MQ=255)                   ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGt                                                                                                                          <  2:723326/150‑1 (MQ=255)                   ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGt                                                                                                                          <  2:689323/150‑1 (MQ=255)                   ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGGCGATGGTGc                                                                                                                                  >  2:499031/1‑142 (MQ=255)                   ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCA‑CGCCACTGGCGATGGTGc                                                                                                                                  <  1:499031/142‑1 (MQ=255)                    cGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTg                                                                                                                         <  1:1090082/150‑1 (MQ=255)                     gACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGa                                                                                                                        <  2:1089155/150‑1 (MQ=255)                      aCGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGGGCACAAAGGTGAc                                                                                                                       <  2:1417593/150‑1 (MQ=255)                                  aTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tggt                                                                                                                                    <  2:1335222/125‑1 (MQ=255)                                  aTCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCg‑tggt                                                                                                                                    >  1:1335222/1‑125 (MQ=255)                                     cTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGAc                                                                                                                       >  2:132892/1‑135 (MQ=255)                                     cTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGAc                                                                                                                       <  1:132892/135‑1 (MQ=255)                                      tCGAAACCGTCGAGGTTCCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                       >  2:720013/1‑150 (MQ=255)                                            ccGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                 <  1:433083/150‑1 (MQ=255)                                            ccGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCg                                                                                                  >  1:815217/1‑149 (MQ=255)                                            ccGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCg                                                                                                  <  2:815217/149‑1 (MQ=255)                                                      ttCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGagcag                                                                                       <  1:237189/150‑1 (MQ=255)                                                         ggCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCACACCCTGAACTTCAGCGGCTTCTCGGTCAACATCGCCAGTGGTG‑TGGTGAAAAAAGGTGACGAGATTTTCGTGCTGCCGTCGGGCAAGAGCAACCg                                                                                    >  2:352414/1‑150 (MQ=255)                                                          gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                   <  1:1129967/150‑1 (MQ=255)                                                          gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                   >  2:1275987/1‑150 (MQ=255)                                                                           cGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCAc                                                                  <  1:720013/150‑1 (MQ=255)                                                                                   ttttCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                 <  1:786278/110‑1 (MQ=255)                                                                                   gTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                 >  2:786278/1‑111 (MQ=255)                                                                                                gTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaaaa                                             >  1:1339218/1‑150 (MQ=255)                                                                                                        aCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATAGACAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGcc                                     >  1:175752/1‑150 (MQ=255)                                                                                                         ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGccc                                    >  2:846125/1‑150 (MQ=255)                                                                                                               cAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                   <  2:1171397/97‑1 (MQ=255)                                                                                                               cAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                   >  1:1171397/1‑97 (MQ=255)                                                                                                                   cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGtg                                                                                                           >  1:556027/1‑69 (MQ=255)                                                                                                                   cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGtg                                                                                                           <  2:556027/69‑1 (MQ=255)                                                                                                                    tGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGt                                                        <  2:1284324/119‑1 (MQ=255)                                                                                                                    tGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGt                                                        >  1:1284324/1‑119 (MQ=255)                                                                                                                       gcTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGa                      <  2:1339218/149‑1 (MQ=255)                                                                                                                         ttCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGCTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGAcc                    >  1:890552/1‑150 (MQ=255)                                                                                                                             gcgGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                 <  1:854223/69‑1 (MQ=255)                                                                                                                             gcgGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                 >  2:854223/1‑69 (MQ=255)                                                                                                                             gcgGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaa                                               <  2:907391/119‑1 (MQ=255)                                                                                                                             gcgGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaa                                               >  1:907391/1‑119 (MQ=255)                                                                                                                                     gcgGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGaa        <  1:846125/150‑1 (MQ=255)                                                                                                                                        ggggccgtCGGCAGGGGCG‑TGGTGGACAAAGGGGACGAGATTGTCGTGCTGCCGGCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAAGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGAc     <  2:175752/143‑1 (MQ=255)                                                                                                                                            accacGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACgag  >  1:83668/5‑150 (MQ=255)                                                                                                                                              aTCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                        <  2:228972/109‑1 (MQ=255)                                                                                                                                              aTCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                        >  1:228972/1‑109 (MQ=255)                                                                                                                                               tctccAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGa         >  1:1425469/4‑139 (MQ=255)                                                                                                                                               tCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGa         <  2:1425469/139‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             TCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCG‑TGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAG  >  NC_002947/1490996‑1491283

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