Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,293,967 C→G 7.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,9670CG 7.0% 183.0 / 3.9 57D13H (GAC→CAC) PP_3763precorrin‑6A synthase
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base G (2/2);  total (27/30)
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.09e-01

GCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCA  >  NC_002947/4293821‑4294113
                                                                                                                                                    |                                                                                                                                                  
gCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGAc                                                                                                                                                  >  5:299036/1‑149 (MQ=255)
       ccaccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGACCCGCGcc                                                                                                                                           <  6:299036/149‑1 (MQ=255)
       ccaccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGcc                                                                                                                                           >  3:243527/1‑149 (MQ=255)
         accCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGa                                                                                                                                         >  5:222937/1‑149 (MQ=255)
         accCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCA‑CCGCGCCGa                                                                                                                                         >  8:106963/1‑148 (MQ=255)
          cccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGa                                                                                                                                         >  2:220784/1‑148 (MQ=255)
                          gcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAgg                                                                                                                        <  6:213594/149‑1 (MQ=255)
                              cAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGcc                                                                                                                                           >  1:163929/1‑126 (MQ=255)
                              cAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGcc                                                                                                                                           <  2:163929/126‑1 (MQ=255)
                               aGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCCCCCGCGCCGAGACCGCTCAACAACAGGTCgtt                                                                                                                   >  8:229536/1‑149 (MQ=255)
                               aGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCttt                                                                                                                   <  5:331382/149‑1 (MQ=255)
                                gAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGTGAAATGCGCAGCGCTTCGACGGGTGCTTTGGTGACATGGCGCGGGGCACCCGCGCGGATCCCGAGCAAAAACGGGTGtt                                                                                                                    >  1:279711/1‑149 (MQ=255)
                                   tGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCAtt                                                                                                               <  4:243527/149‑1 (MQ=255)
                                       cTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGACCCgcgc                                                                                                                                            >  7:265283/1‑116 (MQ=255)
                                       cTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGACCCgcgc                                                                                                                                            <  8:265283/116‑1 (MQ=255)
                                        ttGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccct                                                                                                          >  3:253672/1‑149 (MQ=255)
                                                  gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCg                                                                                                <  5:256593/149‑1 (MQ=255)
                                                   cACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCgg                                                                                               <  5:113070/149‑1 (MQ=255)
                                                              cTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGAc                                                                                    >  1:126661/1‑149 (MQ=255)
                                                               tttGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGa                                                                                     <  7:319794/147‑1 (MQ=255)
                                                               tttGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGa                                                                                     >  8:319794/1‑147 (MQ=255)
                                                                   tcgtcgCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTg                                                                               >  3:180180/1‑149 (MQ=255)
                                                                     gtcgCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCA‑CCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCg                                                                            <  7:106963/149‑1 (MQ=255)
                                                                          ccTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCg                                                                        >  3:202118/1‑149 (MQ=255)
                                                                              gTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccctgcc                                                                                                       <  4:250407/114‑1 (MQ=255)
                                                                              gTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccctgcc                                                                                                       >  3:250407/1‑114 (MQ=255)
                                                                                 gAGCACGAAGATCACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCTCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCg                                                                 <  6:195069/149‑1 (MQ=255)
                                                                                 gAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCg                                                                 <  2:19790/149‑1 (MQ=255)
                                                                                   gCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCgg                                                                >  5:83442/1‑148 (MQ=255)
                                                                                       gaagaaCACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTAttt                                                           >  2:236584/1‑149 (MQ=255)
                                                                                        aagaaCACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTg                                                          <  2:126661/149‑1 (MQ=255)
                                                                                         agaaCACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTgg                                                         <  1:220784/149‑1 (MQ=255)
                                                                                         agaaCACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTgg                                                         <  2:279711/149‑1 (MQ=255)
                                                                                                gCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAag                                                    >  6:305463/1‑147 (MQ=255)
                                                                                                gCTG‑‑GCACTGCGCAGCACATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAagca                                                  <  1:186089/149‑1 (MQ=255)
                                                                                                   g‑‑gCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAgct                                               >  2:33898/1‑149 (MQ=255)
                                                                                                      gCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAgct                                               <  1:33898/148‑1 (MQ=255)
                                                                                                            cgcagcgcaTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTaga                                        <  4:180180/149‑1 (MQ=255)
                                                                                                            cgcagcgcaTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTag                                         <  7:103737/148‑1 (MQ=255)
                                                                                                            cgcagcgcaTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTag                                         >  6:128937/1‑148 (MQ=255)
                                                                                                                gcgcaTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTag                                         >  8:125288/1‑144 (MQ=255)
                                                                                                                   caTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATc                                   >  6:344690/1‑147 (MQ=255)
                                                                                                                     tCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCgaga                               <  6:222937/149‑1 (MQ=255)
                                                                                                                        aCCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGt                            <  1:236584/149‑1 (MQ=255)
                                                                                                                        aCCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGt                            <  5:305463/149‑1 (MQ=255)
                                                                                                                         ccGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGtt                           <  4:253672/149‑1 (MQ=255)
                                                                                                                            cTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGACAAGTTg                          >  3:85509/1‑147 (MQ=255)
                                                                                                                             ttCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTg                                                          <  6:127830/114‑1 (MQ=255)
                                                                                                                             ttCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTg                                                          >  5:127830/1‑114 (MQ=255)
                                                                                                                               cATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAg                       <  5:128937/147‑1 (MQ=255)
                                                                                                                               cATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCt                     >  1:203101/1‑149 (MQ=255)
                                                                                                                                aTAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCt                     >  6:287458/1‑148 (MQ=255)
                                                                                                                                    gTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACaa                <  7:125288/149‑1 (MQ=255)
                                                                                                                                       aCCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGa             <  5:287458/149‑1 (MQ=255)
                                                                                                                                            ccgcgcgGCACCCGCGCCGATACCGAGCAACAACAGGGCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGAccc          <  4:202118/146‑1 (MQ=255)
                                                                                                                                                gggTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGtct    <  6:83442/149‑1 (MQ=255)
                                                                                                                                                  gTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCa  >  8:315177/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTG‑‑GCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCA  >  NC_002947/4293821‑4294113

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