Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,628,161	G→C	8.6%	intergenic (+90/‑50)	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,161	0	G	C	8.6%	136.9 / 2.4	46	intergenic (+90/‑50)	PP_3198/PP_3199	2Fe‑2S ferredoxin/monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base C (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGCCTTGCGCCGCAAGGGCATCGCCGTCGATCTGG  >  NC_002947/3628028‑3628307                                                                                                                                      |                                                                                                                                                        cgggcggCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACg                                                                                                                                          <  7:215792/149‑1 (MQ=255)         tggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGca                                                                                                                                   >  6:70743/1‑148 (MQ=255)         tggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGCCCGCCCCGCATGACgggg                                                                                                                                       <  2:153288/144‑1 (MQ=255)         tggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGCCCGCCCCGCATGACgggg                                                                                                                                       <  2:158786/144‑1 (MQ=255)         tggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGCCCGCCCCGCATGACgggg                                                                                                                                       >  1:158786/1‑144 (MQ=255)         tggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGCCCGCCCCGCATGACgggg                                                                                                                                       >  1:153288/1‑144 (MQ=255)             tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCacaac                                                                                                                               >  4:314156/1‑148 (MQ=255)              gctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCACGCGGGAAGCTCTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCTCacaaca                                                                                                                              >  8:179097/1‑148 (MQ=255)                  ttGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGaa                                                                                                                          <  1:317003/148‑1 (MQ=255)                  ttGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGaa                                                                                                                          <  5:108686/148‑1 (MQ=255)                    gCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACa                                                                                                                        >  6:185507/1‑148 (MQ=255)                     ccGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACaaa                                                                                                                      <  5:185507/149‑1 (MQ=255)                           gAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGt                                                                                                                  >  2:320795/1‑147 (MQ=255)                            agcagcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                                                                                                                <  1:320795/148‑1 (MQ=255)                            agcagcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                                                                                                                <  5:70743/148‑1 (MQ=255)                                gcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACAt                                                                                                           >  1:325104/1‑149 (MQ=255)                                     aTCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCAt                                                                                                      <  7:57217/149‑1 (MQ=255)                                             ccTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtg                                                                                              >  1:323646/1‑149 (MQ=255)                                                   gcgccgcgcGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaac                                                                                         <  3:314156/148‑1 (MQ=255)                                                       cgcgcgATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑gt                                                                                >  6:95290/1‑149 (MQ=255)                                                         cgcgATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑gtgt                                                                              <  2:325104/149‑1 (MQ=255)                                                            gATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCacaac                                                                                                                               >  4:27075/1‑101 (MQ=255)                                                            gATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCacaac                                                                                                                               <  3:27075/101‑1 (MQ=255)                                                                 gTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATc                                                                        >  4:187560/1‑146 (MQ=255)                                                                 gTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATc                                                                        <  3:187560/146‑1 (MQ=255)                                                                      aTCCCGCGACGCGTGTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGggcgg                                                                >  7:305199/1‑149 (MQ=255)                                                                      aTCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGggcgg                                                                <  6:240342/149‑1 (MQ=255)                                                                      aTCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGggcgg                                                                >  6:369150/1‑149 (MQ=255)                                                                         cctggacgcgtTTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGgcggcg                                                              <  8:305199/144‑1 (MQ=255)                                                                                  gaTTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCgg                                                    <  5:95290/149‑1 (MQ=255)                                                                                  gaTTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCAGCATGACGGGGCGCACAACAAGAACAAAACGAGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCgg                                                    >  3:77200/1‑149 (MQ=255)                                                                                          aTGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGc                                                      >  1:1347/1‑139 (MQ=255)                                                                                          aTGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGc                                                      <  2:1347/139‑1 (MQ=255)                                                                                            ggACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTgcc                                          >  7:68139/1‑149 (MQ=255)                                                                                              aCCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTgccgc                                        >  7:261049/1‑149 (MQ=255)                                                                                                cGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCa                                      >  7:19922/1‑149 (MQ=255)                                                                                                    gcAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGc                                  <  5:369150/149‑1 (MQ=255)                                                                                                         cGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGCCTTgc                             <  8:68139/149‑1 (MQ=255)                                                                                                          gggAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGCCTTgcg                            >  5:179335/1‑149 (MQ=255)                                                                                                                   tgCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGCCTTGCGCCGCAAggg                   <  6:179335/149‑1 (MQ=255)                                                                                                                          gCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGCCTTGCGCCGCAAGGGCATCGcc            <  8:261049/149‑1 (MQ=255)                                                                                                                              aGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACAGGGGGGTATGGATGTGAGCGGCAAAATTGGCGGGCTGTGAGCCCACATCGAATGGCGCCGCAGGGGCATCGc             >  1:115915/1‑149 (MQ=255)                                                                                                                              aGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGCCTTGCGCCGCAAGGGCATCGCCGTCg        >  7:299944/1‑149 (MQ=255)                                                                                                                              aGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCAGCCATCGCCTTGCGCCGCAAGGGCATCGCCGTc         >  2:240947/1‑148 (MQ=255)                                                                                                                                gAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGCCTTGCGCCGCAAGGGCATCGCCGTCGAt      <  8:299944/149‑1 (MQ=255)                                                                                                                                    cGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCTAGGGCGTCGCAATTTGAGTAAAATATTCAGCCATCGCCTTGCGCCGCAAGGGCATCGCCGTCGATCTgg  >  8:220729/1‑149 (MQ=255)                                                                                                                                      |                                                                                                                                                        CGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAAC‑‑‑‑GTGT‑TGATCGTGGGCGGCGGAATTGGCGGGCTGTGTGCCGCCATCGCCTTGCGCCGCAAGGGCATCGCCGTCGATCTGG  >  NC_002947/3628028‑3628307

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