Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,981,108 A→C 23.1% 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,981,1080AC23.1% 23.4 / 2.2 29A251A (GCT→GCG) PP_3506magnesium chelatase subunit ChII
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base C (5/2);  total (16/13)
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

GCTGTGGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCG  >  NC_002947/3980978‑3981240
                                                                                                                                  |                                                                                                                                    
gcTGTGGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCgcgggccccgcg                                                                                                                                 >  1:71800/1‑136 (MQ=255)
     ggCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGAACAGTGGGCGcgcgcagcgagcagac                                                                                                                            >  1:232898/1‑134 (MQ=255)
                    cTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCCCGCAGCCACCCCAGGTCGgcgg                                                                                                             >  2:492534/1‑135 (MQ=255)
                      ggCGTGGGGCGTCATACGGGGGGGAGGGCGCTGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCTATAGCTGCAGCCCCACGCCTGGCACTGTGGGGGGGCGCTGCGCGCAGCCACACCAGGTCGgcgcgc                                                                                                           <  1:494798/136‑1 (MQ=255)
                      ggCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCCCCCGGGCGgcgggc                                                                                                           >  1:220141/1‑136 (MQ=255)
                      ggCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCCCGCAGCCCCCCCCGGTCGgcgccc                                                                                                           >  1:305804/1‑134 (MQ=255)
                      ggCCTGGGGCGGCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTGCTTCGATAGCTGCCGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGgcgcgc                                                                                                           <  1:330931/136‑1 (MQ=255)
                                   aTACGGCGGCGATGGCGCAGGGCGAATTCGGCCCCCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGc                                                                                              <  1:137663/136‑1 (MQ=255)
                                   aTACGGCGGCGAGGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTCCCCCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGc                                                                                              <  1:246190/136‑1 (MQ=255)
                                    tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCCCCCGCCGCCACACCACGGCGGCGCCCCGGCCCTCCAcccc                                                                                             >  1:491452/1‑136 (MQ=255)
                                              aTGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTag                                                                                   >  1:241450/1‑136 (MQ=255)
                                              aTGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCAGTCCACGCCGGCGGCGTag                                                                                   >  1:241859/1‑136 (MQ=255)
                                                ggCGCTGGGCGATTGCGCCCACCGCCCGCACATCCACTTCTTCTATCGCGGCCCCCCCCCGCCTGCAACTGGGGGCGCGCGCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTagca                                                                                 <  2:373379/136‑1 (MQ=255)
                                                ggCGCCGGGCTAATGCGCCCACCCCCTCGACATCCACTTCTTCGATCGCTGCCGCCCCCCGCCCGGCACTGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTagca                                                                                 <  2:332686/136‑1 (MQ=255)
                                                     tgggCGAATTCGGCCACCGCCGCGACATCCACTTCTTCGATCGCTGCAGCCCCACGCCAGGCACAGGGGGGGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCt                                                                            <  2:22628/135‑1 (MQ=255)
                                                     aGGGCGAATGCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCt                                                                            <  2:334595/136‑1 (MQ=255)
                                                                gccccccGCCTCTAAATCCACTTCTTCGATAGCTGCAGCCCCCCGCCCGGCACAGTGGGCGCGCGAAGCGCGCAGCCACACCAGGTCGGCGCTCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAg                                                                 <  1:190123/131‑1 (MQ=255)
                                                                    cccgccgCTACATCCACTTCTTCGATAGCTCCACCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCaa                                                             <  2:500419/135‑1 (MQ=255)
                                                                     ccGCATCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAg                                                            >  1:455584/1‑136 (MQ=255)
                                                                              aaTCCACTTCTTCGAAAGCTCCACCCCCACGCCCGCCACTGGGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTc                                                   <  2:455584/135‑1 (MQ=255)
                                                                                    cttcttCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAg                                             >  2:384072/1‑136 (MQ=255)
                                                                                    cttcttCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAg                                             >  2:381618/1‑136 (MQ=255)
                                                                                                     gCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCCGCGGCGTAGCAGCGCTCCGTAATCCAGGCAAGCGCCAGGTCATCCAGGGCAATGTAGGCCaacg                            >  2:258944/1‑134 (MQ=255)
                                                                                                           ccccAGCCACCGGGGGCGCGCGCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGgc                      <  2:241450/134‑1 (MQ=255)
                                                                                                           cGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCAGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGgc                      <  2:241859/136‑1 (MQ=255)
                                                                                                            gCCAGGCACAGTGGGCGCGCGCCGCGCGCAGCCCCACCAGGTAGGGGCGCCGGCCCTGCAAGCCCGCGGCGGAGCCGGGGTCGGGTAAACAGGCAAGCGCCAGGTTCTGGAGGGCGATTTGGGCCCAAGCCAGgcc                     >  2:318798/1‑135 (MQ=255)
                                                                                                              cAGGCACAGTGGGCGCGCGCCGCGCGCAGCCACCCCAGGTCGGCGCGCAGGCCGTCCACCCCCGCGGCGGAGCCGCGCGCGGTAATACCGGCAAGCGCCTGGTCCTCCAGGGCAAAGTCGGCCAAGGCCAGgcgcc                   >  1:41070/1‑135 (MQ=255)
                                                                                                                   acaGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCcgcg              >  2:245790/1‑136 (MQ=255)
                                                                                                                        gggCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCa         >  1:100440/1‑136 (MQ=255)
                                                                                                                          gcgcgcgcAGCGCGCCGCCCCCCCAGGTCGGCGCGCCGGCCGTCCACGCCGGCGGCGTAGCCGCGCGCGGTAATCCCGGCAAGCGCCGGGGGATCCAGGGGAATGTCGGCCAAGGCCGGGCGCGCCGCGGGGCCgc       >  2:106602/1‑136 (MQ=255)
                                                                                                                          gcgcgcgcAGCGCGCAGCCCCACCAGGTCGGCGCGCAGGCCGGCCACGCCGGCGGGGTTGCAGCGCTCGGTTATCCAGGCCAGCGCCTGGTGATCCAGGGCCATGTCGGCCAAGGCCGGGCGCGCCGCGTGGCCgc       >  2:419289/1‑136 (MQ=255)
                                                                                                                             ggcgcagcgcgcagcCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCt    <  1:392097/135‑1 (MQ=255)
                                                                                                                              gcgcagcgcgcagcCACACCAGGTCGGCGCGCCGGCCGTCCACGCCGGCGGCGTAGCCGCGCGCGGTAATCCAGGCAAGCGCCGGGGCATCCAGGGGAATGTCGGCCCAGGCCAGGCGCGCCGCGTGGCCGCGCTc   >  2:238093/1‑136 (MQ=255)
                                                                                                                               cgcagcgcgcagcCACACCAGGTTGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCCGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTcg  >  2:485589/1‑136 (MQ=255)
                                                                                                                                  |                                                                                                                                    
GCTGTGGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCG  >  NC_002947/3980978‑3981240

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