Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 878,918 C→A 15.8% G148G (GGC→GGA PP_0763 → medium‑chain‑fatty acid CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947878,9180CA15.8% 60.0 / 3.0 25G148G (GGC→GGAPP_0763medium‑chain‑fatty acid CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (13/8);  new base A (2/2);  total (15/10)
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.85e-01

GCCTGTCGCCCGAGCAGATCCTCTACACCATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCAGCCCGCGCTACGACTTCCCGGATTTCGACGAGAATTCGGTGGCCACTACCTTCTACACCACTGGCACCACCGGTAA  >  NC_002947/878776‑879058
                                                                                                                                              |                                                                                                                                            
gcCTGTCGCCCGAGCAGATCCTCTACACCATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTAGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACa                                                                                                                                       >  2:512089/1‑150 (MQ=255)
   tGTCGCCCGAGCAGATCCTCTACACCATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGa                                                                                                                                    <  1:313861/150‑1 (MQ=255)
        cccGAGCAGATCCTCTACACCATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGcc                                                                                                                               <  1:613606/150‑1 (MQ=255)
           gAGCAGATCCTCTACACCATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGaa                                                                                                                            >  1:97661/1‑150 (MQ=255)
                           ccATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACaa                                                                                                                                      <  1:124821/124‑1 (MQ=255)
                           ccATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACaa                                                                                                                                      >  2:124821/1‑124 (MQ=255)
                            cATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTAGCGGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGAACCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCgg                                                                                                           >  1:419935/1‑150 (MQ=255)
                                        cGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGCCTGGTc                                                                                                             <  1:466180/136‑1 (MQ=255)
                                        cGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGCCTGGTc                                                                                                             >  2:466180/1‑136 (MQ=255)
                                              ccGCGTGGTGCTGGTCGTTAGTGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAgctgag                                                                                         <  2:419935/150‑3 (MQ=255)
                                               cgcgTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGa                                                                                                        <  2:65306/134‑1 (MQ=255)
                                               cgcgTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGa                                                                                                        >  1:65306/1‑134 (MQ=255)
                                                           gTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGAACACTGTCGACAAGACCCTGCTACTGACCGGTGGCCCGTACAGGACTGCCTGACTGCCCGATCTGTGGTGGGTGTATGAGCAGCTGCTGGCTGCTGCCAGc                                                                            >  2:130477/1‑150 (MQ=255)
                                                                          gtcgGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCAGCCCGCGCTACGACTTc                                                             >  2:481445/1‑150 (MQ=255)
                                                                                cTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCAGCCCGCGCTACGACTTCCCGGAt                                                       >  1:434206/1‑150 (MQ=255)
                                                                                             tCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGATGTCTGCTGCCAGCCCGCGCTACGACTTCCCGGATTTCGACGAGAAtt                                          >  1:67248/1‑150 (MQ=255)
                                                                                                             cACTGTCGACAAGACCCTGCTACTGACCGATGGACCGGACAAGACTGc                                                                                                                                >  2:12102/1‑48 (MQ=255)
                                                                                                             cACTGTCGACAAGACCCTGCTACTGACCGATGGACCGGACAAGACTGc                                                                                                                                >  2:321153/1‑48 (MQ=255)
                                                                                                             cACTGTCGACAAGACCCTGCTACTGACCGATGGACCGGACAAGACTGc                                                                                                                                <  1:321153/48‑1 (MQ=255)
                                                                                                             cACTGTCGACAAGACCCTGCTACTGACCGATGGACCGGACAAGACTGc                                                                                                                                <  1:12102/48‑1 (MQ=255)
                                                                                                                tGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCAGCCCGCGCTACGACTTCCCGGATTTCGACGAGAATTCGGTGGCCACTACCTTCTa                       <  1:512089/150‑1 (MQ=255)
                                                                                                                    gACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCAGCCCGCGCTACGACTTCCCGGATTTCGACGAGAATTCGGTGGCCACTACCTTCTAcacc                   >  2:470132/1‑150 (MQ=255)
                                                                                                                                     gACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCa                                                                              >  1:234884/1‑74 (MQ=255)
                                                                                                                                     gACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCa                                                                              <  2:234884/74‑1 (MQ=255)
                                                                                                                                     gACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCAGCCCGCGCTACGACTTCCCGGATTTCGACGAGAATTCGGTGGCCACTACCTTCTACACCACTGGCACCACCGGTaa  >  2:66217/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GCCTGTCGCCCGAGCAGATCCTCTACACCATGAACCATGCCGAAGACCGCGTGGTGCTGGTCAACAGCGACTTCGTCGGCCTGTACCAGGCCATCGCCGGGCAGCTGACCACTGTCGACAAGACCCTGCTACTGACCGATGGCCCGGACAAGACTGCCGAACTGCCCGGTCTGGTCGGCGAGTATGAGCAGCTGCTGGCTGCTGCCAGCCCGCGCTACGACTTCCCGGATTTCGACGAGAATTCGGTGGCCACTACCTTCTACACCACTGGCACCACCGGTAA  >  NC_002947/878776‑879058

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: