Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,842,904:1	+G	9.8%	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,904	1	.	G	9.8%	140.8 / 11.8	41	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 . (22/15);  new base G (2/2);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATAC  >  NC_002947/4842758‑4843047                                                                                                                                                    |                                                                                                                                                ccGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                <  3:2005/148‑1 (MQ=255)           aacaacTGCAGGGCAACGACAAGTTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCATGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCggg                                                                                                                                     <  6:297893/149‑1 (MQ=38)            acaacTGCAGGGCAACGACAAGGTGCACTTGCTGGTGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGt                                                                                                                                    <  6:253323/149‑1 (MQ=255)                       gCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTAc                                                                                                                                >  4:283190/1‑142 (MQ=255)                       gCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTAc                                                                                                                                <  3:283190/142‑1 (MQ=255)                         aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAACCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGaa                                                                                                                         >  7:59378/1‑147 (MQ=255)                              cacggggCACTTGCGGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTACGGCAGGCGTGGCCTCTTCGCGGGTAACCCCGCTCCTACACAGCATCCT‑TGTCGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                  <  2:188723/143‑1 (MQ=37)                               aaGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                 <  4:226731/149‑1 (MQ=255)                                     cACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTCCACAGGAGCCT‑TGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTgg                                                                                                           <  6:276211/149‑1 (MQ=255)                                                  cTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTAACCCCCTCCTACCCAGGCTGCG‑TGTTGGCGCGGGGTTTCCCGCGAAGCGGCCGGGTCTGGGCCAGCAGTtctc                                                                                              >  1:118004/1‑149 (MQ=255)                                                           aCAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTGCTACACAGGATCCG‑TGTTGGTGCGGGGTTTCCCGCGAAGAGGCCGGTTATGGGACAGCAATTGTAGCAGGCTAg                                                                                     >  7:157638/1‑149 (MQ=25)                                                               ccTGTCGCACCAGGCTGGGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                >  2:32225/1‑86 (MQ=255)                                                               ccTGTCGCACCAGGCTGGGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                <  1:32225/86‑1 (MQ=255)                                                                   tCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACAGTATCCT‑TGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTgg                                                                             <  8:67698/149‑1 (MQ=255)                                                                     gCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGa                                                                            >  8:54073/1‑148 (MQ=255)                                                                     gCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGAt                                                                           <  7:54073/149‑1 (MQ=255)                                                                      cACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGCAGCGGCCGGTACTGGCCCAGCCTTTCTCGCAGGCAACTGACCgggct                                                                           >  3:33841/1‑146 (MQ=37)                                                                               gAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGacac                                                                 >  5:295809/1‑149 (MQ=255)                                                                                    tgCCAAGGGTGCTTAAGGCTGGCGTGGCCTCTGCGGGGTTAAACCCGCTCCAACACCGGTTCCT‑TGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGc                                                             <  3:294251/147‑1 (MQ=37)                                                                                          ggCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGt                                                                                    >  1:28174/1‑119 (MQ=255)                                                                                          ggCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGt                                                                                    <  2:28174/119‑1 (MQ=255)                                                                                          ggCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAaca                                                      >  1:251943/1‑149 (MQ=255)                                                                                          ggCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAATTCTCGCAGGCAAGTGACCTGGATCAATTGAACACCGCCATAaca                                                      >  8:163579/1‑149 (MQ=255)                                                                                          ggCTCAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTAGCACAGCATTTCTCGCAGGCAAGTGa                                                                                  <  1:248671/121‑1 (MQ=255)                                                                                             tGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGATGCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGt                                                   >  7:19652/1‑149 (MQ=255)                                                                                                 tAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAGTTCTCGCAGGGAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCa                                               >  6:225392/1‑149 (MQ=255)                                                                                                   aggctggcGTGGCCTCTTCGCGGGTAAACCCGCTCGTACACTGGATCCT‑TGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCAt                                              <  8:59378/148‑1 (MQ=38)                                                                                                           ctggcctCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCAc                                         <  7:70948/145‑1 (MQ=255)                                                                                                           ctggcctCTTCGCGGGTAAACCCGCTCCTACACAGGATACT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCAc                                         >  8:70948/1‑145 (MQ=255)                                                                                                                    tCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTtc                                                                                                >  1:64752/1‑81 (MQ=255)                                                                                                                    tCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTtc                                                                                                <  2:64752/81‑1 (MQ=255)                                                                                                                        gggTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTa                         >  2:234863/1‑148 (MQ=255)                                                                                                                          gTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAg                      <  7:163579/149‑1 (MQ=255)                                                                                                                            aaaCCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGt                    >  6:264576/1‑149 (MQ=255)                                                                                                                               cccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCg                  >  6:276326/1‑148 (MQ=255)                                                                                                                                ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCg                  >  8:225901/1‑147 (MQ=255)                                                                                                                                  gCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGt              >  3:19692/1‑149 (MQ=255)                                                                                                                                  gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGt              >  3:36391/1‑149 (MQ=255)                                                                                                                                  gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGt              >  8:148370/1‑149 (MQ=255)                                                                                                                                   cTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCg                  <  2:45880/144‑1 (MQ=255)                                                                                                                                   cTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCAGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCg                  >  1:45880/1‑144 (MQ=255)                                                                                                                                       ctacaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGc            <  4:36391/142‑1 (MQ=255)                                                                                                                                          acaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCTTTTCTCGCAGGCAAGTGACCTGAATCAATTGACACACGCCATAACACGTTTCATACCACGAACCAAACGCGACTAAAGGTCGAAGTTCGAAcc       >  5:16885/1‑148 (MQ=39)                                                                                                                                            aGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa     >  2:242609/1‑148 (MQ=255)                                                                                                                                               tcctt‑tGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATAc  <  7:189919/144‑1 (MQ=255)                                                                                                                                                   t‑tGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAc                                            >  5:176490/1‑102 (MQ=255)                                                                                                                                                   t‑tGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAc                                            <  6:176490/102‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                CCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATAC  >  NC_002947/4842758‑4843047

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