Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,976,607 C→T 8.7% A87A (GCG→GCA)  PP_3506 ← magnesium chelatase subunit ChII

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,976,6070CT 8.7% 125.8 / 3.3 46A87A (GCG→GCA) PP_3506magnesium chelatase subunit ChII
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base T (2/2);  total (25/21)
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.96e-01

GGTGCGAGATACCGTCACGCTCGATGCGGTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAGGGTGCTCTTGGCCATGCCACGCGG  >  NC_002947/3976459‑3976753
                                                                                                                                                    |                                                                                                                                                  
ggTGCGAGATACCGTCACGCTCGATGCGGTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTgc                                                                                                                                                    >  6:173645/1‑149 (MQ=255)
                aCGCTCGATGCGGTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTtgccctgccc                                                                                                                                      >  2:224370/1‑147 (MQ=255)
                  gCTCGATGCGGTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAgc                                                                                                                                  <  1:224370/149‑1 (MQ=255)
                      gATGGGGTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGAGCGTCCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACCCCGCCGTCGGCCTGGGCCAGCTCGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGcgcg                                                                                                                               <  7:281387/148‑1 (MQ=255)
                      gATGCGGTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCgg                                                                                                                              >  7:283928/1‑149 (MQ=255)
                       aTGCGGTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGc                                                                                                                             <  5:173645/149‑1 (MQ=255)
                            gTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGTGCCTTGCCCTGCCCCAGCGCGGCATc                                                                                                                          >  3:83375/1‑147 (MQ=255)
                             tttGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGTGCCTTGCCCTGCCCCAGCGCGGCATCCAg                                                                                                                       <  4:83375/149‑1 (MQ=255)
                             tttGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAg                                                                                                                       <  1:55252/149‑1 (MQ=255)
                                   ccGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAg                                                                                                                 >  6:287124/1‑149 (MQ=255)
                                           aGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGTGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAg                                                                                                                 >  4:61305/1‑141 (MQ=255)
                                           aGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGTGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAg                                                                                                                 <  3:61305/141‑1 (MQ=255)
                                            gCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGAc                                                                                                        >  1:324095/1‑149 (MQ=255)
                                                cGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGc                                                                                                                             >  2:249317/1‑124 (MQ=255)
                                                          cAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCAc                                                                                          >  4:335503/1‑149 (MQ=255)
                                                                 aCCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATcc                                                                                                                         >  2:182791/1‑111 (MQ=255)
                                                                 aCCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATcc                                                                                                                         <  1:182791/111‑1 (MQ=255)
                                                                   cAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGccc                                                                                  >  8:105826/1‑148 (MQ=255)
                                                                     aGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAAt                                                                               >  5:103605/1‑149 (MQ=255)
                                                                            gggCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGcc                                                                                                           >  1:344276/1‑114 (MQ=255)
                                                                            gggCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGcc                                                                                                           <  2:344276/114‑1 (MQ=255)
                                                                                    ggTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAAc                                                                <  2:324095/149‑1 (MQ=255)
                                                                                     gTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGc                                                                            >  8:94645/1‑136 (MQ=255)
                                                                                     gTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGc                                                                            <  7:94645/136‑1 (MQ=255)
                                                                                      ttGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACgg                                                              <  7:105826/149‑1 (MQ=255)
                                                                                        gACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTg                                                                                        <  7:112600/121‑1 (MQ=255)
                                                                                        gACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTg                                                                                        >  8:112600/1‑121 (MQ=255)
                                                                                            tCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTc                                                        >  2:347727/1‑149 (MQ=255)
                                                                                              aTCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGAc                                                                    <  4:213729/135‑1 (MQ=255)
                                                                                              aTCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGAc                                                                    >  3:213729/1‑135 (MQ=255)
                                                                                               tCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGcc                                                     <  6:103605/149‑1 (MQ=255)
                                                                                                      aCAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGt                                                                       >  7:134923/1‑124 (MQ=255)
                                                                                                      aCAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGt                                                                       <  8:134923/124‑1 (MQ=255)
                                                                                                        aGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTc                                            >  5:855/1‑149 (MQ=255)
                                                                                                         gCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGc                                                            >  4:27481/1‑132 (MQ=255)
                                                                                                         gCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGc                                                            <  3:27481/132‑1 (MQ=255)
                                                                                                           acacCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTc                                                        >  6:157110/1‑134 (MQ=255)
                                                                                                           acacCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTc                                                        <  5:157110/134‑1 (MQ=255)
                                                                                                                        ccTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGc                             >  6:316269/1‑148 (MQ=255)
                                                                                                                            ggCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAggg                        <  3:335503/149‑1 (MQ=255)
                                                                                                                            ggCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAgg                         <  6:855/148‑1 (MQ=255)
                                                                                                                               cAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAGGGTGc                     >  4:122591/1‑149 (MQ=255)
                                                                                                                                    cGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAGGGTGCTCTTg                <  3:110893/149‑1 (MQ=255)
                                                                                                                                      cccGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAGGGTGCTCTTGGc              <  8:283928/149‑1 (MQ=255)
                                                                                                                                             aaaaCTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAGGGTGCTCTTGGCCATGCCa       >  4:122415/1‑149 (MQ=255)
                                                                                                                                                 cTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAGGGTGCTCTTGGCCATGCCAcgcg   >  3:285575/1‑149 (MQ=255)
                                                                                                                                                  tGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAGGGTGCTCTTGGCCATGCCACGCgg  >  8:180042/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTGCGAGATACCGTCACGCTCGATGCGGTTTGTGCCGCTGGCAGCCACGTCAAGCAACAGATCGACCAAGGTGTCGGGCAGCAGGTTGACCTCATCGACATACAGCACACCGCCGTCGGCCTGGGCCAGCACGCCCGGGGAAAACTGCGCCTTGCCCTGCCCCAGCGCGGCATCCAGGTCGAGGGTGCCGACCAGCCGCTCTTCACTGGCGCCCAATGGCAAGGTGACGAACGGGCCTTCGCCAAGCAGGTCGGCCAGGCCACGGGCCAGGGTGCTCTTGGCCATGCCACGCGG  >  NC_002947/3976459‑3976753

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