Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,294,011 C→G 17.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_0029474,294,0110CG17.9% 57.9 / 8.3 28D13H (GAC→CAC) PP_3763precorrin‑6A synthase
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base G (3/2);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

CAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGC  >  NC_002947/4293862‑4294151
                                                                                                                                                     |                                                                                                                                            
cAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTc                                                                                                                                              >  1:161954/1‑150 (MQ=255)
 aGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCa                                                                                                                                             <  2:37624/150‑1 (MQ=255)
     gATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCAccta                                                                                                                                         <  2:304655/150‑3 (MQ=255)
                             gcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGt                                                                                                                 >  2:475409/1‑150 (MQ=255)
                                  aGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTc                                                                                                            <  2:161954/150‑1 (MQ=255)
                                           ttGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCTCGGGTCACCCGCGCCGATACCGATCAACAACAg                                                                                                                   >  2:225848/1‑134 (MQ=255)
                                           ttGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAg                                                                                                                   <  1:225848/134‑1 (MQ=255)
                                               gcAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCAtt                                                                                                         >  1:261083/1‑140 (MQ=255)
                                               gcAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCAtt                                                                                                         <  2:261083/140‑1 (MQ=255)
                                               gcAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTgccctgccct                                                                                               >  2:94178/1‑150 (MQ=255)
                                                     gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGACCCGCGccag                                                                                                                                   <  2:422431/108‑3 (MQ=255)
                                                     gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGACCCGCGccag                                                                                                                                   >  1:422431/1‑106 (MQ=255)
                                                              cATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGc                                                                                           >  2:431852/1‑139 (MQ=255)
                                                              cATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGc                                                                                           <  1:431852/139‑1 (MQ=255)
                                                                    tGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGGACCCGCGACCAtt                                                                          <  1:59310/150‑1 (MQ=255)
                                                                                    gAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCgg                                                          <  2:432218/150‑1 (MQ=255)
                                                                                         cGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTAttt                                                     >  1:579752/1‑150 (MQ=255)
                                                                                          gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccct                                                                                                    >  2:349056/1‑102 (MQ=255)
                                                                                          gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccct                                                                                                    <  1:349056/102‑1 (MQ=255)
                                                                                                     tGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTc                                                                                                            >  2:77969/1‑83 (MQ=255)
                                                                                                     tGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTc                                                                                                            <  1:77969/83‑1 (MQ=255)
                                                                                                          cTGGGCAGGGCAGCGACCGCTTGATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTa                                    <  2:40091/150‑1 (MQ=255)
                                                                                                           tGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTag                                   <  1:94178/150‑1 (MQ=255)
                                                                                                                     aTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTagaaga                               <  1:117243/144‑1 (MQ=255)
                                                                                                                     aTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTagaaga                               >  2:117243/1‑144 (MQ=255)
                                                                                                                            gCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTa                  <  1:475409/150‑1 (MQ=255)
                                                                                                                                            ggCGCGGGTGACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCcgc  >  2:365580/1‑150 (MQ=255)
                                                                                                                                             gcgcGGGTGACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAccc                                                                                   >  2:188976/1‑68 (MQ=255)
                                                                                                                                             gcgcGGGTGACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAccc                                                                                   <  1:188976/68‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
CAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGC  >  NC_002947/4293862‑4294151

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: