Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,130,169 G→C 16.7% A107G (GCC→GGC)  PP_5378 ← metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,130,1690GC16.7% 52.8 / 6.5 24A107G (GCC→GGC) PP_5378metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (9/11);  new base C (2/2);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

ATCTGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTC  >  NC_002947/6130022‑6130311
                                                                                                                                                   |                                                                                                                                              
aTCTGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCa                                                                                                                                              >  1:531621/1‑150 (MQ=255)
      ggCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTg                                                                                                                                        <  2:531621/150‑1 (MQ=255)
         aTTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGc                                                                                                                                     <  1:366695/150‑1 (MQ=255)
          tttGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCg                                                                                                                                    <  1:194606/150‑1 (MQ=255)
               tCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGCCAAGATTgg                                                                                                                                       >  1:437958/1‑142 (MQ=255)
               tCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGCCAAGATTgg                                                                                                                                       <  2:437958/142‑1 (MQ=255)
                 aaGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAAAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGGGGGGGCGGTGCCgg                                                                                                                             >  1:281236/1‑150 (MQ=255)
                      aaGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATaaa                                                                                                                        <  2:562704/150‑1 (MQ=255)
                       aGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTACCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAg                                                                                                                       <  1:486545/150‑1 (MQ=255)
                           aCAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCt                                                                                                                   <  1:81452/150‑1 (MQ=255)
                                                                         gCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCaa                                                                     >  2:420528/1‑150 (MQ=255)
                                                                                         gATACCGTGTTTAATTGTCCAGAACGCAGCTGCGGGGTCATCGATTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTa                                                     <  2:484568/150‑1 (MQ=255)
                                                                                             ccGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTa                                                 >  2:572600/1‑150 (MQ=255)
                                                                                                                    aactgctgGGTCTTCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGt                          >  2:289948/3‑150 (MQ=255)
                                                                                                                         cGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATaaa                                                                                                                        >  2:111405/1‑51 (MQ=255)
                                                                                                                         cGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATaaa                                                                                                                        <  1:111405/51‑1 (MQ=255)
                                                                                                                                aTCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCAACGGCACAGATAAGGGTTTAATCCTTTAg              >  2:269586/1‑150 (MQ=255)
                                                                                                                                    cTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCTCTGATAAGGGTTTCATCCTTTAGGTTc          <  1:420528/150‑1 (MQ=255)
                                                                                                                                     ttAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCtt                 >  2:390728/1‑142 (MQ=255)
                                                                                                                                     ttAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCtt                 <  1:390728/142‑1 (MQ=255)
                                                                                                                                         aCCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAAc     >  2:350937/1‑150 (MQ=255)
                                                                                                                                          ccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATc                    >  2:140091/1‑134 (MQ=255)
                                                                                                                                          ccAATCTTGCCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATc                    <  1:140091/134‑1 (MQ=255)
                                                                                                                                            aaTCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTc  <  1:572600/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ATCTGCGGCATTTGCTCAAGGAAAGCGACAATACCCCAGATGTACTGATCACCCATGCTCTTGCCCCAAGCGGGCATACCGGTAGCCTTGATACCGTGTTTAATGGTCCAGAACGCAGCTGCGGGGTCATCGCTTAAACCAATCTTGGCAAGATTGGGCGGTGCCGGATAAAGCGCTTGGCTGAGCTCGGTTTTATCCACCCCTGGGGCGAGGTGGCAGCCAATGCACATAGCATTGTAGTTACCGGCACCGGCACGGATAAGGGTTTCATCCTTTAGGTTCGGAACCTC  >  NC_002947/6130022‑6130311

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