Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,628,082	C→G	9.5%	intergenic (+11/‑129)	PP_3198 → / → PP_3199	2Fe‑2S ferredoxin/monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,628,082	0	C	G	9.5%	124.9 / 4.3	42	intergenic (+11/‑129)	PP_3198/PP_3199	2Fe‑2S ferredoxin/monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (20/18);  new base G (2/2);  total (22/20)
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.86e-01

GCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTG  >  NC_002947/3627944‑3628222                                                                                                                                           |                                                                                                                                              gcgcgcgcTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAgg                                                                                                                                     <  8:170740/149‑1 (MQ=255)  cgcgcgcTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGACGCCTGGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGt                                                                                                                                    >  5:250399/1‑149 (MQ=255)                 cGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcgacgc                                                                                                                     >  7:99625/1‑149 (MQ=255)                  gAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATATGTCGCATCCcgcgacgc                                                                                                                     >  6:265632/1‑148 (MQ=255)                                 gTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCg                                                                                                     >  8:102704/1‑149 (MQ=255)                                   tCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGgc                                                                                                   >  7:160133/1‑149 (MQ=255)                                     aCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATcgcgccgcgc                                                                                                                                           >  3:176499/1‑107 (MQ=255)                                     aCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATcgcgccgcgc                                                                                                                                           <  4:176499/107‑1 (MQ=255)                                      cGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcg                                                                                                                           <  3:148689/122‑1 (MQ=255)                                      cGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcg                                                                                                                           >  4:148689/1‑122 (MQ=255)                                               ccGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCTATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                        <  7:102704/148‑1 (MQ=255)                                               ccGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                        <  8:191912/148‑1 (MQ=255)                                               ccGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAg                                                                                       <  8:160133/149‑1 (MQ=255)                                                 gAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGAt                                                                                                                                        >  6:318773/1‑98 (MQ=255)                                                 gAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGAt                                                                                                                                        <  5:318773/98‑1 (MQ=255)                                                  aGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATATGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCtg                                                                                    <  5:265632/149‑1 (MQ=255)                                                    gCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCTATAGGTCGCATCCCGCGACGCCCTTGGTCATGGACCGGCGCAGGCGGGAAGCtgag                                                                                  >  7:328179/1‑147 (MQ=255)                                                                     cTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                        <  2:325008/126‑1 (MQ=255)                                                                          tCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgccgc                                                            >  6:216839/1‑149 (MQ=255)                                                                              aaGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCg                                                        >  6:47827/1‑149 (MQ=255)                                                                              aaGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCg                                                        >  5:90807/1‑149 (MQ=255)                                                                                 gcaGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CAt                                                    <  6:250399/147‑1 (MQ=255)                                                                                            ctggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGTCATG‑CGGGgcggc                                         <  6:90807/149‑3 (MQ=255)                                                                                               gctgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAg                                                                       <  5:236883/117‑1 (MQ=255)                                                                                               gctgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAg                                                                       >  6:236883/1‑117 (MQ=255)                                                                                                ctgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGa                                                                                         >  5:230147/1‑98 (MQ=255)                                                                                                ctgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGa                                                                                         <  6:230147/98‑1 (MQ=255)                                                                                                 tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCacaac                                     >  2:272307/1‑148 (MQ=255)                                                                                                   ctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCacaaca                                    <  5:47827/147‑1 (MQ=255)                                                                                                   ctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCAcaacaa                                   <  5:216839/148‑1 (MQ=255)                                                                                                      ttGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGaa                                >  3:222400/1‑148 (MQ=255)                                                                                                                agcagcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGTga                      <  8:270978/148‑1 (MQ=255)                                                                                                                     cTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGGACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAGAACGTGAGACATg                >  5:183412/1‑149 (MQ=255)                                                                                                                     cTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGTGAGACATg                >  4:102476/1‑149 (MQ=255)                                                                                                                        gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAg                                  <  4:178441/128‑1 (MQ=255)                                                                                                                        gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAg                                  >  3:178441/1‑128 (MQ=255)                                                                                                                            gCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGt                        >  4:89856/1‑134 (MQ=255)                                                                                                                            gCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGt                        <  3:89856/134‑1 (MQ=255)                                                                                                                              cgACCTATCGCGCGGCGCGATAGGTCACATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCg                                                                                            >  6:233597/2‑65 (MQ=37)                                                                                                                              cgACCTATCGCGCGGCGCGATAGGTCACATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCg                                                                                            <  5:233597/64‑1 (MQ=37)                                                                                                                                aCCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACt     >  5:202928/1‑149 (MQ=255)                                                                                                                                aCCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACt     <  1:272307/149‑1 (MQ=255)                                                                                                                                  cTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgt   >  1:273633/1‑149 (MQ=255)                                                                                                                                     tCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgtg  <  2:336745/147‑1 (MQ=255)                                                                                                                                           |                                                                                                                                              GCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCG‑CATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTG  >  NC_002947/3627944‑3628222

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