Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,217,251 C→A 12.3% P92P (CCC→CCA moaA → cyclic pyranopterin monophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,217,2510CA12.3% 82.8 / 2.9 33P92P (CCC→CCAmoaAcyclic pyranopterin monophosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (15/13);  new base A (2/2);  total (17/16)
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.65e-01

CAGCAAATCCTCAGCCTCGAAGAACTGTTCCAGGTGGCCGAACGCTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCAGCCTGGATAGCCTGGATGCCGACCGTTTCAAGCAGCTGACCCGCACCGGCGAC  >  NC_002947/5217108‑5217398
                                                                                                                                               |                                                                                                                                                   
cAGCAAATCCTCAGCGTCGAAGAACTGTTCCAGGTGGCCGAACGCTTCGTCGCTCTGGGCAGCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCAGCGCCGCCCTGCCCGGCCTg                                                                                                                                               >  2:37501/1‑150 (MQ=255)
      aTCCTCAGCCTCGAAGAACTGTTCCAGGTGGCCGAACGCTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAg                                                                                                                                         <  1:544018/150‑1 (MQ=255)
              ccTCGAAGAACTGTTCCAGGTGGCCGAACGCTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGCGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCt                                                                                                                                 >  2:464385/1‑150 (MQ=255)
                 cGAAGAACTGTTCCAGGTGGCCGAACGCTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGAc                                                                                                                              <  1:344778/150‑1 (MQ=255)
                        cTGTTCCAGGTGGCCGAACGCTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACAAGCAAc                                                                                                                       <  1:82270/150‑1 (MQ=255)
                                ggTGGCCGAACGCTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCa                                                                                                               >  2:397990/1‑150 (MQ=255)
                                      cGAACGCTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCgg                                                                                                         >  1:78574/1‑150 (MQ=255)
                                            cTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCt                                                                                                   >  1:89876/1‑150 (MQ=255)
                                                     gcTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCggg                                                                                                        >  2:207299/1‑136 (MQ=255)
                                                     gcTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCggg                                                                                                        <  1:207299/136‑1 (MQ=255)
                                                      cTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGc                                                                                                                                        >  2:546163/1‑103 (MQ=255)
                                                      cTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGc                                                                                                                                        <  1:546163/103‑1 (MQ=255)
                                                           cACCCGAAAGATCAGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGACTCGCAACTCGGGCGCCTGGCTCAGCCGCTGtt                                                                                    >  1:293065/1‑150 (MQ=255)
                                                              ccGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGa                                                                                 >  1:433926/1‑150 (MQ=255)
                                                                     aTCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCggg                                                                          <  2:169248/150‑1 (MQ=255)
                                                                         gCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGa                                                                      >  2:192818/1‑150 (MQ=255)
                                                                         gCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCAGGCCTgc                                                                                                                                              <  1:40223/78‑1 (MQ=255)
                                                                         gCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCAGGCCTgc                                                                                                                                              <  1:39809/78‑1 (MQ=255)
                                                                         gCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCAGGCCTgc                                                                                                                                              >  2:39809/1‑78 (MQ=255)
                                                                         gCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCAGGCCTgc                                                                                                                                              >  2:40223/1‑78 (MQ=255)
                                                                                gggTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGCGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCt                                                               <  1:464385/150‑1 (MQ=255)
                                                                                       gAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGAc                                                                                >  1:202618/1‑126 (MQ=255)
                                                                                       gAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGAc                                                                                <  2:202618/126‑1 (MQ=255)
                                                                                                     ccAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCagcctggatagcct                                          <  2:78574/150‑1 (MQ=255)
                                                                                                            aTCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGAAGCCGGGGTGGCACGCCTGAACATCAGCCTGGATAGCCTGGATGcc                                   >  1:349223/1‑150 (MQ=255)
                                                                                                                ccgACCTGTGCCGGCGCACCCCCTCCCTCCCTGGACTGCGCGAGCTGTTCCTGACCAGTAACGGCCCGCAACTCGGGCGCCTGTTTCAGCCGCTGTTCTACGCCGGGGTGACACGCCTTAACATCAGCCTGGATAGCCTGGATGCCGAcc                               <  2:240752/149‑1 (MQ=255)
                                                                                                                   aCCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGtt                                                                                    <  2:269046/94‑1 (MQ=255)
                                                                                                                   aCCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGtt                                                                                    >  1:269046/1‑94 (MQ=255)
                                                                                                                     cTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCAGCCTGGATAGCCTGGATGCCGACCGTTTc                          <  1:283105/150‑1 (MQ=255)
                                                                                                                        tgCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCAGCCTGGATAGCCTGGATGCCGACCGTTTCAag                       >  1:303410/1‑150 (MQ=255)
                                                                                                                         gCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCAGCCTGGATAGCCTGGATGCCGACCGTTTCAagc                      <  2:616581/150‑1 (MQ=255)
                                                                                                                          cGGGCGCCTCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCAGCCTGGATAGCCTGGATGCCGACCGTTTCAagca                     <  1:183458/150‑1 (MQ=255)
                                                                                                                                        cccTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCAGCCTGGATAGCCTGGATGCCGACCGTTTCAAGCAGCTGACCCGCACCg       >  2:92816/1‑150 (MQ=255)
                                                                                                                                             cccGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCAGCCTGGATAGCCTGGATGCCGACCGTTTCAAGCAGCTGACCCGCACCGGCGAc  <  1:397990/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAGCAAATCCTCAGCCTCGAAGAACTGTTCCAGGTGGCCGAACGCTTCGTCGCGCTGGGCACCCGAAAGATCCGCCTGACGGGTGGCGAGCCACTGGTACGCCAGGGCATCGTCGACCTGTGCGGGCGCATCGCCGCCCTGCCCGGCCTGCGCGAGCTGTGCCTGACCAGCAACGGCTCGCAACTCGGGCGCCTGGCTCAGCCGCTGTTCGACGCCGGGGTGACACGCCTGAACATCAGCCTGGATAGCCTGGATGCCGACCGTTTCAAGCAGCTGACCCGCACCGGCGAC  >  NC_002947/5217108‑5217398

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