Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,842,949:1	+G	22.0%	intergenic (+48/‑249)	ccoP‑II → / → PP_4259	cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,842,949	1	.	G	22.0%	148.9 / 54.7	59	intergenic (+48/‑249)	ccoP‑II/PP_4259	cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (23/23);  new base G (8/5);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACT  >  NC_002947/4842799‑4843093                                                                                                                                                        |                                                                                                                                                 gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCGTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCt                                                                                                                                                   <  2:221547/150‑1 (MQ=255)   tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                 <  1:1297100/150‑1 (MQ=255)     ccGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGt                                                                                                                                               <  2:171751/150‑1 (MQ=255)      cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTa                                                                                                                                              <  2:255676/150‑1 (MQ=255)       ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCCTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAg                                                                                                                                             <  2:2056887/150‑1 (MQ=255)             ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTTTTCCACCAGACGAAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                       <  1:1135845/150‑1 (MQ=37)             ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg                                                                                                                                       <  2:323026/150‑1 (MQ=255)             ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                       <  2:1454964/150‑1 (MQ=255)              gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCggg                                                                                                                                      <  2:2123682/150‑1 (MQ=255)              gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCggg                                                                                                                                      <  1:907685/150‑1 (MQ=255)                     tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTcccc                                                                                                                               <  2:1328685/150‑1 (MQ=255)                     tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTAccc                                                                                                                               <  2:2175000/150‑1 (MQ=255)                          ggCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAACGATGCCCTGGGTAGGCCTCGCTGGTTGGTGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGaa                                                                                                                          <  1:2676429/150‑1 (MQ=25)                           gCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGCTCCT‑TGTAGGAGCGGGTTTACCCGCGAag                                                                                                                         >  1:118141/1‑150 (MQ=255)                                      gTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGaa                                                                                                                          <  2:95719/138‑1 (MQ=255)                                      gTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGaa                                                                                                                          >  1:95719/1‑138 (MQ=255)                                            ttGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGccg                                                                                                                                       >  1:1946876/1‑117 (MQ=255)                                            ttGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGccg                                                                                                                                       <  2:1946876/119‑3 (MQ=255)                                            ttGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGc                                                                                                                                         >  1:1502188/1‑117 (MQ=255)                                            ttGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGc                                                                                                                                         <  2:1502188/117‑1 (MQ=255)                                                tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAg                                                                                                                                             <  1:1201292/109‑1 (MQ=255)                                                tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAg                                                                                                                                             >  2:1201292/1‑109 (MQ=255)                                                tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTAc                                                                                                                                 >  1:322155/1‑121 (MQ=255)                                                tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTAc                                                                                                                                 <  2:322155/121‑1 (MQ=255)                                                       ttACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGc                                                                                                                                         >  2:2536404/1‑106 (MQ=25)                                                       ttACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGc                                                                                                                                         <  1:2536404/106‑1 (MQ=25)                                                              tACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGt                                                                                                                                     >  1:3038799/1‑103 (MQ=255)                                                              tACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGt                                                                                                                                     <  2:3038799/103‑1 (MQ=255)                                                                                                                       ttCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                       >  2:2577980/1‑108 (MQ=255)                                                                                                                       ttCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                       <  1:2577980/108‑1 (MQ=255)                                                                                                                             ggTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAg                       >  2:1549225/1‑150 (MQ=255)                                                                                                                                  aCCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTg                                                                                                             >  1:2598040/1‑59 (MQ=255)                                                                                                                                  aCCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTg                                                                                                             <  2:2598040/59‑1 (MQ=255)                                                                                                                                  aCCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGc                  >  2:1752934/1‑150 (MQ=255)                                                                                                                                    ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATaa                                                         >  2:1548677/1‑109 (MQ=255)                                                                                                                                    ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATaa                                                         <  1:1548677/109‑1 (MQ=255)                                                                                                                                    ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                >  2:2935252/1‑150 (MQ=255)                                                                                                                                    ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                >  2:2544275/1‑150 (MQ=255)                                                                                                                                    ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                >  1:693581/1‑150 (MQ=255)                                                                                                                                    ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAcc                                            >  2:1197766/1‑122 (MQ=255)                                                                                                                                    ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAcc                                            <  1:1197766/122‑1 (MQ=255)                                                                                                                                    ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                >  1:2549829/1‑150 (MQ=255)                                                                                                                                    ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                >  1:2911687/1‑150 (MQ=255)                                                                                                                                      gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg              >  2:1451451/1‑150 (MQ=255)                                                                                                                                      gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg              >  1:2457296/1‑150 (MQ=255)                                                                                                                                          cTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCg                                                               >  2:1803389/1‑97 (MQ=255)                                                                                                                                          cTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCg                                                               <  1:1803389/97‑1 (MQ=255)                                                                                                                                           tACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAttt                                                                                                  <  1:1098026/61‑1 (MQ=255)                                                                                                                                           tACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAttt                                                                                                  >  2:1098026/1‑61 (MQ=255)                                                                                                                                            acacAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGttt                                                  >  1:2610742/1‑108 (MQ=255)                                                                                                                                            acacAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGttt                                                  <  2:2610742/108‑1 (MQ=255)                                                                                                                                            acacAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAcc        >  2:100049/1‑150 (MQ=255)                                                                                                                                             cacaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAccc       >  2:2094794/1‑150 (MQ=255)                                                                                                                                             acaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGAAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAg                       >  1:539075/4‑135 (MQ=255)                                                                                                                                             acaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGAAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAg                       <  2:539075/132‑1 (MQ=255)                                                                                                                                               aaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa      >  2:1348416/2‑150 (MQ=255)                                                                                                                                                aGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                       >  2:2661272/1‑83 (MQ=255)                                                                                                                                                aGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                       <  1:2661272/83‑1 (MQ=255)                                                                                                                                                aGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAAcg                                <  2:2915490/122‑1 (MQ=255)                                                                                                                                                aGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAAcg                                >  1:2915490/1‑122 (MQ=255)                                                                                                                                                  gATCCA‑TGCAGGACCGGGTTTACCCGCGAAGAGGCCGGTACTAGCACAGCAGTTCTCGCAGGCAGGTGACCTGGCTCCATTGACACCCGCCATACCACGTTTCCTACCACGAACCCAACCCCACCAACGGTCGAAGCGCGAAccacgact  >  2:172181/1‑144 (MQ=255)                                                                                                                                                        |                                                                                                                                                 GATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACT  >  NC_002947/4842799‑4843093

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