Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,293,967	C→G	6.9%	D13H (GAC→CAC)	PP_3763 ←	precorrin‑6A synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,293,967	0	C	G	6.9%	177.4 / 3.8	58	D13H (GAC→CAC)	PP_3763	precorrin‑6A synthase
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base G (2/2);  total (30/28)
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.49e-01

GCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGC  >  NC_002947/4293821‑4294107                                                                                                                                                    |                                                                                                                                             gCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTCGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCAc                                                                                                                                            >  1:116634/1‑149 (MQ=255)   gATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCAccc                                                                                                                                          <  8:282135/149‑1 (MQ=255)            ccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATa                                                                                                                                 >  5:354174/1‑149 (MQ=255)            ccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATa                                                                                                                                 >  7:311031/1‑149 (MQ=255)             cTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATAc                                                                                                                                >  8:316062/1‑149 (MQ=255)                ggggCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCg                                                                                                                              >  2:168038/1‑148 (MQ=255)                       tAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATcaacaa                                                                                                                      <  1:279711/149‑1 (MQ=255)                         gcgcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCaacaac                                                                                                                     <  3:20360/148‑1 (MQ=255)                         gcgcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCaacaac                                                                                                                     <  8:311031/148‑1 (MQ=255)                          cgcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAg                                                                                                                   <  8:43287/149‑1 (MQ=255)                            cTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGt                                                                                                                 >  6:117157/1‑149 (MQ=255)                            cTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGt                                                                                                                 >  2:50336/1‑149 (MQ=255)                               cAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCtt                                                                                                              >  1:333448/1‑149 (MQ=255)                                      gCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtg                                                                                                        <  6:354174/148‑1 (MQ=255)                                       ccTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgcc                                                                                                      <  2:116634/149‑1 (MQ=255)                                             gcAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccctgccc                                                                                                >  8:307173/1‑149 (MQ=255)                                                       cAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCCTCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGc                                                                                           <  5:342315/144‑1 (MQ=255)                                                       cAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGCCCGCTTCAT‑AGGTGACCTGTCGCGGGTCCCCCGCGCCGAGACCGATCAACAACCGGTCTTTCAGTGCCCGGCCCTTTGCgcgct                                                                                      >  6:342315/1‑145 (MQ=255)                                                       cAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCCGCCCTTTGCGGGaa                                                                                      <  7:307173/149‑1 (MQ=255)                                                         gTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAcc                                                                                    >  3:143339/1‑149 (MQ=255)                                                          ttCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAcc                                                                                    <  4:143339/148‑1 (MQ=255)                                                             aTCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCcgcg                                                                                >  6:105838/1‑149 (MQ=255)                                                                tttGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCCCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGAc                                                                              >  8:292580/1‑148 (MQ=255)                                                                  tGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGAc                                                                              <  1:50336/146‑1 (MQ=255)                                                                          cccTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATAc                                                                                                                                <  2:131572/88‑1 (MQ=255)                                                                          cccTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATAc                                                                                                                                >  1:131572/1‑88 (MQ=255)                                                                             ttGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCAGCCGCGCCGATACCGAGCAACAACAGGGCTTTCCTTGCCCTGCCCTTTGCGGGCACCCGCGACCGTTGgcgggcggg                                                                >  8:329427/1‑148 (MQ=255)                                                                                tCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTAc                                                             >  7:330871/1‑149 (MQ=255)                                                                                    gCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGc                                                                                           <  7:260632/115‑1 (MQ=255)                                                                                    gCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGc                                                                                           >  8:260632/1‑115 (MQ=255)                                                                                        gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTGACCCgcgc                                                                                                                                      >  7:307050/1‑68 (MQ=255)                                                                                        gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTGACCCgcgc                                                                                                                                      <  8:307050/68‑1 (MQ=255)                                                                                           gaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGc                                                  <  5:117157/149‑1 (MQ=255)                                                                                              cacGCTGGCTCTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCaaa                                               <  1:168038/149‑1 (MQ=255)                                                                                                   tGGCACTGCGCAGCGCATCGACCGCTTCATGAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAagcag                                           >  1:214167/1‑149 (MQ=255)                                                                                                       aCTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAgctgct                                      >  4:144410/1‑149 (MQ=255)                                                                                                        cTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTAGTTGGCAACGCAgctgct                                      >  6:356309/1‑148 (MQ=255)                                                                                                           cgcagcgcaTCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTaga                                  <  5:105838/149‑1 (MQ=255)                                                                                                           cgcagcgcaTCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTa                                    <  7:292580/147‑1 (MQ=255)                                                                                                                cgcaTCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATc                             <  2:333448/149‑1 (MQ=255)                                                                                                                 gcaTCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCg                            <  1:207485/149‑1 (MQ=255)                                                                                                                 gcaTCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCg                            <  7:230891/149‑1 (MQ=255)                                                                                                                  caTCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATc                             >  1:42657/1‑147 (MQ=255)                                                                                                                    tCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCgag                          >  6:52796/1‑148 (MQ=255)                                                                                                                    tCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCgag                          >  3:227218/1‑148 (MQ=255)                                                                                                                      gACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCgag                          >  6:201404/1‑146 (MQ=255)                                                                                                                      gACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGaa                        <  5:201404/148‑1 (MQ=255)                                                                                                                      gACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGaa                        <  7:316062/148‑1 (MQ=255)                                                                                                                       aCCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGt                      >  3:198160/1‑149 (MQ=255)                                                                                                                         cGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTg                    >  5:268624/1‑149 (MQ=255)                                                                                                                          gCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGt                   <  2:42657/149‑1 (MQ=255)                                                                                                                           cTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTg                    >  5:86811/1‑147 (MQ=255)                                                                                                                           cTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTa                  <  2:214167/149‑1 (MQ=255)                                                                                                                                   ggTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGAc            <  8:330871/148‑1 (MQ=255)                                                                                                                                      gtatcggCGCGGGTGACCCGCGCCAATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCga                           <  5:101987/125‑1 (MQ=255)                                                                                                                                      gtatcggCGCGGGTGACCCGCGCCAATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCga                           >  6:101987/6‑130 (MQ=255)                                                                                                                                            gcgcGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCcgc  >  2:132736/1‑149 (MQ=255)                                                                                                                                                  gTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTa                  <  4:97207/127‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             GCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAT‑AGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGC  >  NC_002947/4293821‑4294107

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