Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,978,081 A→C 32.6% A251A (GCT→GCG PP_3506 ← magnesium chelatase subunit ChII

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,978,0810AC32.6% 20.2 / 2.1 24A251A (GCT→GCGPP_3506magnesium chelatase subunit ChII
Reads supporting (aligned to +/- strand):  ref base A (5/9);  new base C (4/4);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGC  >  NC_002947/3977954‑3978214
                                                                                                                               |                                                                                                                                     
ggggCCCCCACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGcgcgcagcgcgcag                                                                                                                               <  1:318527/134‑1 (MQ=255)
          gggcggGCTGGCCTGGGGCGTCATACGGCGGCTTTGGCGCAGGGCTAATTCGCCCACCGCCTCGAAATCCACTTCTTCGATAGCTGCCGCCCCCCGCCAGGCACTGTGGGCGCGCGCAGGGCGCAGCCACACCAgg                                                                                                                     <  2:275988/135‑1 (MQ=255)
          gggcggGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAgg                                                                                                                     <  2:317288/135‑1 (MQ=255)
                                aTACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGc                                                                                               <  1:241588/136‑1 (MQ=255)
                                aTACGGCGGCGAGGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGc                                                                                               <  2:134486/136‑1 (MQ=255)
                                 tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCCCCCCCCCCCCCACCAGGTCGGCCCGCCCGCCGGCCACGcc                                                                                              >  2:464183/1‑136 (MQ=255)
                                 tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCAcccc                                                                                              >  2:223833/1‑136 (MQ=255)
                                 tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGcc                                                                                              >  1:164046/1‑136 (MQ=255)
                                   cggcggCGAGGCCCCAGGGCCAATCCGCCCACCGCCCCCACACCCACTTCTCCCATCCCGCCCCCCCCCCCCCCGCCACTGGGGGCGCGCGCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCcgg                                                                                            <  2:507629/136‑1 (MQ=255)
                                           aTGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCACCAGGTCGGCGCGCACGCCGTCCACGCCGGCGGCGTag                                                                                    >  1:146581/1‑136 (MQ=255)
                                             ggCGCAGGGCTAATTCGGCCCCGCCCCCGACATCCACTTCTTCGAACGCTCCAGCCCCCCGCCTGGCACTGTGGGCGCGCGCAGCGCGAAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTagca                                                                                  <  2:144021/136‑1 (MQ=255)
                                                      cGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGt                                                                         <  2:129511/136‑1 (MQ=255)
                                                       gAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCACCCCCACCAAGTCCGCGCGCACGCCGTCCACCCCGGCCGCGGAGCAGCGCGCGGTa                                                                        >  1:139685/1‑136 (MQ=255)
                                                        aaTTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCAGGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTaa                                                                       <  1:322296/136‑1 (MQ=255)
                                                           tCGGCCACCGCCTCGACATCCACTTCTTCGATCGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATcc                                                                    <  2:208959/136‑1 (MQ=255)
                                                           gcgcccACCGCCTCCACACCCACTTCTGCGATCGAGCCCGCCCCACGCCCGCCACTGGGGGCGGGCGCTGGGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATcc                                                                    <  1:300235/132‑1 (MQ=255)
                                                            cGGCCCCCGCGTCGACATCCACTTCTTCTATCGCTGCCGCCCCCCGCCCGGCACTGTGGGCGCGCGCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCa                                                                   <  1:442596/136‑1 (MQ=255)
                                                                cACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCa                                                               <  1:140114/136‑1 (MQ=255)
                                                                 aCCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCCGCCCCACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCACGGTAATCCAGGCaa                                                              >  2:279792/1‑136 (MQ=255)
                                                                                   tcttcGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCCCCCCGCCCCACCACGGCGGCGCGCAGGCCGTCCACCCCCGCCGCGTAGCCGCGCACGGTAATACACGCAAGCCCCAGGTCATCCAggg                                            >  2:531045/1‑136 (MQ=255)
                                                                                      tcGATAGCTGCCGCCCCACGCCTGGCACAGGGGGCGCGCGCTGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCaa                                         <  2:537320/136‑1 (MQ=255)
                                                                                      cctaTAGCGGCGCCCCCCCCCCCGGCACGGGGGGGGCGCCCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCaa                                         <  1:216438/133‑1 (MQ=255)
                                                                                                    cccACGCCAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCCCCCGGTCGGCGCGCCGGCCGGCCCCGCCGGCGGCGGAGGCGGGCGCGGGTATACAGGCCAGCGCCGGGGGCTCCCGGGCCAAGTCGGGCCAGGcc                           >  2:338278/1‑136 (MQ=255)
                                                                                                        cGCCCGGCACCGGGGGCGCGCGCCGCGCGCCGCCACACCAGGGCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGgc                       <  1:391119/136‑1 (MQ=255)
                                                                                                          ccAGGCACAGTGGGCGCGGGCCCCGCGCCGCCCCACCCGGGGGGCGCGCCGGGCGGCCACCCCCGCGGCGTTGCCGGGCGGGGGTATACCGGCCAGCGCCGGGTGATCCCGGGCCATGTGGGCCCAGGCCAGgggc                     >  2:406999/1‑136 (MQ=255)
                                                                                                                        cgcgcgccgcgcgcCGCCACCCCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCGCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAgcg       >  2:482746/1‑136 (MQ=255)
                                                                                                                             gcagcgcgcagcCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTcgc  >  1:520136/1‑136 (MQ=255)
                                                                                                                               |                                                                                                                                     
GTGGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGC  >  NC_002947/3977954‑3978214

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: