Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,499,756	G→A	10.5%	A178A (GCG→GCA)	cobK →	precorrin‑6x reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,499,756	0	G	A	10.5%	102.9 / 2.7	38	A178A (GCG→GCA)	cobK	precorrin‑6x reductase
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base A (2/2);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

TCACGCTAGGACGCGAGCCGCTGCAGCACCTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCA  >  NC_002947/5499608‑5499900                                                                                                                                                     |                                                                                                                                                  tCACGCTAGGACGCGAGCCGCTGCAGCACCTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCg                                                                                                                                                   >  6:332226/1‑149 (MQ=255)   aCGCTAGGACGCGAGCCGCTGCAGCACCTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCgt‑t                                                                                                                                                >  8:225674/1‑149 (MQ=255)     gCTAGGACGCGAGCCGCTGCAGCACCTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCgt‑tgt                                                                                                                                              >  2:105910/1‑149 (MQ=255)          gACGCGAGCCGCTGCAGCACCTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAg                                                                                                                                         <  4:42420/149‑1 (MQ=255)                 gCCGCTGCAGCACCTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCgccgcc                                                                                                                                  <  3:291207/149‑1 (MQ=255)                         agcaCCTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGAt                                                                                                                          <  2:214614/149‑1 (MQ=255)                             ccTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCgt‑tgt                                                                                                                                              >  3:291147/1‑125 (MQ=255)                             ccTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCgt‑tgt                                                                                                                                              <  4:291147/125‑1 (MQ=255)                                    gagaTTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAg                                                                                                               >  2:109858/1‑149 (MQ=255)                                    gagaTTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAg                                                                                                               >  3:39263/1‑149 (MQ=255)                                     agaTTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGc                                                                                                              <  7:225674/149‑1 (MQ=255)                                                gAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCAGCCACATCGATGTGCTGATCAGCAAGAACAGCgg                                                                                                   <  7:255403/149‑1 (MQ=255)                                                           gACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGa                                                                                        <  1:167633/149‑1 (MQ=255)                                                             ccTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGTCTGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGAc                                                                                       <  1:105910/149‑1 (MQ=255)                                                              cTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTCGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TCTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCg                                                                                     >  8:71876/1‑149 (MQ=255)                                                                  cgtgcgCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGcc                                                                                 <  5:332226/149‑1 (MQ=255)                                                                    tgcgCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGc                                                                                                                      <  2:143872/110‑1 (MQ=255)                                                                       gcTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCAT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGc                                                                            >  1:199160/1‑149 (MQ=255)                                                                        cTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCt                                                                           <  4:130972/149‑1 (MQ=255)                                                                           gAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTgg                                                                         <  2:132017/148‑1 (MQ=255)                                                                              gCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCgt‑tgt                                                                                                                                              >  2:227794/1‑76 (MQ=255)                                                                              gCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCgt‑tgt                                                                                                                                              <  1:227794/76‑1 (MQ=255)                                                                              gCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGt                                                                     <  1:109858/149‑1 (MQ=255)                                                                                     cTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAgc                                                                                                                                        >  5:288565/1‑75 (MQ=255)                                                                                     cTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAgc                                                                                                                                        <  6:288565/75‑1 (MQ=255)                                                                                           aCGAGCGCTGTGAAGGGATGGGCGCACGTGGGCCGGTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCtt                                                        <  6:21208/149‑1 (MQ=255)                                                                                              aGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGc                                                     <  6:122254/149‑1 (MQ=255)                                                                                                     tgAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCgg                                              >  4:272736/1‑149 (MQ=255)                                                                                                      gAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGt                                             >  7:47148/1‑149 (MQ=255)                                                                                                                gcgcACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCt                                    >  5:149591/1‑148 (MQ=255)                                                                                                                gcgcACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCt                                    >  2:102070/1‑148 (MQ=255)                                                                                                                  gcAAGTGGGCCGTTACATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGa                                  >  5:170359/1‑148 (MQ=255)                                                                                                                      gTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATAGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGaac                                                                             >  8:9454/1‑100 (MQ=255)                                                                                                                       tGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTc                            <  6:170359/149‑1 (MQ=255)                                                                                                                          gCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGa                         >  3:127128/1‑149 (MQ=255)                                                                                                                                ccATATCGAGGATGAACGGGCAT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGc                   >  4:297654/1‑149 (MQ=255)                                                                                                                                 cATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGcc                  >  8:356282/1‑149 (MQ=255)                                                                                                                                    atCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGc                   <  7:154842/145‑1 (MQ=255)                                                                                                                                    atCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGc                   >  8:154842/1‑145 (MQ=255)                                                                                                                                       gAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGt            >  6:70185/1‑149 (MQ=255)                                                                                                                                       gAGGATGAACGGGCAT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAgg             <  3:297654/148‑1 (MQ=255)                                                                                                                                            tGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTCCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGa         >  7:45183/1‑147 (MQ=255)                                                                                                                                              aaCGGGCAT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAgg     <  2:199160/149‑1 (MQ=255)                                                                                                                                                 gggCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCa  >  4:308684/1‑149 (MQ=255)                                                                                                                                                 gggCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCa  >  8:26126/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TCACGCTAGGACGCGAGCCGCTGCAGCACCTGGACGAGATTCCGCCGGAGCAGTTCTGGACCTTGCGTGCGCTGGAGGCGTGCCCTGGCAACGAGCGCTGTGAAGTGATTGGCGCACGTGGGCCGTTCCATATCGAGGATGAACGGGCGT‑TGTTCCAGCGCCGCCACATCGATGTGCTGATCAGCAAGAACAGCGGCAGTGTGGCGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCA  >  NC_002947/5499608‑5499900

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: ‑