Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,293,967 C→G 13.0% 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_0029474,293,9670CG13.0% 127.2 / 10.2 46D13H (GAC→CAC) PP_3763precorrin‑6A synthase
Reads supporting (aligned to +/- strand):  ref base C (17/23);  new base G (3/3);  total (20/26)
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.85e-01

GGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCAT  >  NC_002947/4293820‑4294114
                                                                                                                                                   |                                                                                                                                                   
ggCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCa                                                                                                                                                    <  5:268718/149‑1 (MQ=255)
 cgatcctccaccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCAc                                                                                                                                                   <  1:183631/142‑1 (MQ=255)
        cgaccCTCGTGGCGGTAGCGCTGCAGAATTGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGGCACCCGCGcc                                                                                                                                            <  5:209636/147‑1 (MQ=255)
         caccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGACCCgcgc                                                                                                                                             >  3:212729/1‑147 (MQ=255)
          accCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACCCGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGa                                                                                                                                          <  7:151385/149‑1 (MQ=255)
          accCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGACCCGCGccag                                                                                                                                          <  4:212729/149‑3 (MQ=255)
          accCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCCGGCGCGGGTGACCCGCGccag                                                                                                                                          <  2:171047/149‑3 (MQ=255)
                    cGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATcaa                                                                                                                                <  7:210599/149‑1 (MQ=255)
                         gcgcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCaacaac                                                                                                                            <  2:194622/148‑1 (MQ=255)
                                  aaTGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCa                                                                                                                  >  6:264457/1‑149 (MQ=255)
                                        tttGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATcaac                                                                                                                               <  8:155893/129‑1 (MQ=255)
                                                  cgcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGc                                                                                                  <  2:251123/149‑1 (MQ=255)
                                                   gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCGTCATAGGTGCCCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCg                                                                                                 >  6:179725/1‑149 (MQ=255)
                                                       cAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGTAGAAGACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGaa                                                                                             <  4:220288/149‑1 (MQ=255)
                                                        aGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAc                                                                                            >  8:146141/1‑149 (MQ=255)
                                                         gTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAcc                                                                                           >  4:235404/1‑149 (MQ=255)
                                                                  tGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGAc                                                                                     <  5:169544/146‑1 (MQ=255)
                                                                      gtcgCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccc                                                                                                            >  4:237791/1‑119 (MQ=255)
                                                                      gtcgCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccc                                                                                                            <  3:237791/119‑1 (MQ=255)
                                                                       tcgCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGaa                                                                                             <  7:25251/133‑1 (MQ=255)
                                                                       tcgCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGaa                                                                                             >  8:25251/1‑133 (MQ=255)
                                                                          cccTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGc                                                                          <  3:235404/149‑1 (MQ=255)
                                                                           ccTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCg                                                                         >  6:270286/1‑149 (MQ=255)
                                                                                  gAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCg                                                                  <  1:177354/149‑1 (MQ=255)
                                                                                  gAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCg                                                                  >  1:231478/1‑149 (MQ=255)
                                                                                  gAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCg                                                                  <  7:146141/149‑1 (MQ=255)
                                                                                       cGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTAtt                                                             >  4:134685/1‑149 (MQ=255)
                                                                                        gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTAcc                                                                   >  7:75689/1‑142 (MQ=255)
                                                                                        gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTAcc                                                                   <  8:75689/142‑1 (MQ=255)
                                                                                                                      gACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGaa                               <  5:264457/148‑1 (MQ=255)
                                                                                                                      gACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAg                              <  3:134685/149‑1 (MQ=255)
                                                                                                                       aCCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGt                             >  7:41881/1‑149 (MQ=255)
                                                                                                                       aCCGCTTCATAGGTGACCTGGCGCGGGTCACCAGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCGCTTTGCGGGAAACCGCGACCATTTTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTCCTTAGAAGATagag                                 >  2:73417/1‑145 (MQ=255)
                                                                                                                         cGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTg                           <  5:270286/149‑1 (MQ=255)
                                                                                                                         cGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTg                           >  5:185440/1‑149 (MQ=255)
                                                                                                                              cATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAg                        <  2:231478/147‑1 (MQ=255)
                                                                                                                               aTAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCt                      >  7:110818/1‑148 (MQ=255)
                                                                                                                                 aGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAACTGCTTAGAAGATCGAGaa                               <  3:1503/137‑1 (MQ=255)
                                                                                                                                 aGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAACTGCTTAGAAGATCGAGaa                               >  4:1503/1‑137 (MQ=255)
                                                                                                                                  ggTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGAc                   >  2:25351/1‑148 (MQ=255)
                                                                                                                                       ccTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGAc             <  2:83707/149‑1 (MQ=255)
                                                                                                                                        cTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGAcc            >  4:194425/1‑149 (MQ=255)
                                                                                                                                          ggCGCGGGTGACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCgg                                                                 >  1:4450/1‑94 (MQ=255)
                                                                                                                                          ggCGCGGGTGACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCgg                                                                 <  2:4450/94‑1 (MQ=255)
                                                                                                                                          ggCGCGGGTGACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCcg          >  7:5558/1‑149 (MQ=255)
                                                                                                                                                  tCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCAt  <  6:185440/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCAT  >  NC_002947/4293820‑4294114

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