Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,099,491 C→G 11.3% A76G (GCC→GGC)  PP_2709 → long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,099,4910CG11.3% 101.2 / 3.9 35A76G (GCC→GGC) PP_2709long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (12/19);  new base G (2/2);  total (14/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.89e-01

GCGCAGCGTTACACCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCTTTGCCCGGCAGGCGCCAGCTG  >  NC_002947/3099349‑3099640
                                                                                                                                              |                                                                                                                                                      
gcgcAGCGTTACACCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAAcgc                                                                                                                                                <  2:68686/150‑1 (MQ=255)
  gcAGCGTTACACCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACgcgc                                                                                                                                              >  1:1131217/1‑150 (MQ=255)
            acCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACa                                                                                                                                    <  1:536821/150‑1 (MQ=255)
              cTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGt                                                                                                                                  >  2:240340/1‑150 (MQ=255)
                  cGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTgccg                                                                                                                              >  2:65610/1‑150 (MQ=255)
                              aaCGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAgc                                                                                                                  >  2:45799/1‑150 (MQ=255)
                              aaCGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCGTTCGGCGAACAGCGccg                                                                                                                                            <  1:1065350/124‑2 (MQ=255)
                              aaCGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCGTTCGGCGAACAGCGccg                                                                                                                                            >  2:1065350/1‑123 (MQ=255)
                                  aaGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGa                                                                                                              >  1:1068090/1‑150 (MQ=255)
                                     tCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCGTTCGGCGAACAGCGccg                                                                                                                                            >  1:588198/1‑116 (MQ=255)
                                     tCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCGTTCGGCGAACAGCGccg                                                                                                                                            <  2:588198/117‑2 (MQ=255)
                                                        cgcaagcAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTgca                                                                                        <  2:1200984/145‑1 (MQ=255)
                                                             gcAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCaaa                                                                                   <  1:240340/150‑1 (MQ=255)
                                                               aGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGc                                                                                 <  2:862566/150‑1 (MQ=255)
                                                                gCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACTCAGTTCCAGCACTGCTTGCAGCAAAGCg                                                                                <  1:65610/150‑1 (MQ=255)
                                                                    accaccGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTCTTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAgg                                                                            >  1:125146/1‑149 (MQ=255)
                                                                      caccGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAg                                                                             <  1:471292/147‑1 (MQ=255)
                                                                      caccGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAg                                                                             >  2:471292/1‑147 (MQ=255)
                                                                        ccGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGAcg                                                                        <  2:934743/150‑1 (MQ=255)
                                                                         cGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGAcgg                                                                       >  1:127276/1‑150 (MQ=255)
                                                                           ggAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGAcggcg                                                                     <  2:269865/150‑1 (MQ=255)
                                                                               ccTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACTCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCg                                                                 <  2:1068090/150‑1 (MQ=255)
                                                                                cTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCgt                                                                >  2:544565/1‑150 (MQ=255)
                                                                                cTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCACCACTGCTTGCAGCAAATCGGAGGGACGGCGTTCgt                                                                >  2:315957/1‑150 (MQ=255)
                                                                                      cGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGa                                                          <  2:452995/150‑1 (MQ=255)
                                                                                        tggtgCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAgg                                                        <  2:1131217/150‑1 (MQ=255)
                                                                                               ggATAACGGGCCCGAAGCGCTATTCCGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTg                                                 <  1:444553/150‑1 (MQ=255)
                                                                                                   aaCGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGcc                                             <  1:222664/150‑1 (MQ=255)
                                                                                                          ccGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTggg                                               <  1:690368/141‑1 (MQ=255)
                                                                                                          ccGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTggg                                               >  2:690368/1‑141 (MQ=255)
                                                                                                          ccGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGc                                      <  2:975743/150‑1 (MQ=255)
                                                                                                                       cggggATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCTGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTgg                         <  2:525688/148‑1 (MQ=255)
                                                                                                                           gATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCttt                     <  1:315957/150‑1 (MQ=255)
                                                                                                                                        tGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGcc                                             <  2:291531/113‑1 (MQ=255)
                                                                                                                                        tGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGcc                                             >  1:291531/1‑113 (MQ=255)
                                                                                                                                              ccGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCTTTGCCCGGCAGGCGCCAGCTg  >  1:1023137/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
GCGCAGCGTTACACCTACCGGCAACTGCTCAACGAAGTCGAATCCCGTGCGAGCTGCTTGCGCAGCCAACCACCGGGAACCTTCGCCGTGGTGCTGGATAACGGGCCCGAAGCGCTATTCTGGGATTTGGCGGCGCTGTTCGCCGAACGCGCCAGTGTGACAGTGCCGCCGTTTTTCAGCGCGACCCAGTTCCAGCACTGCTTGCAGCAAAGCGGAGTGACGGCGGTCGTGTGTGAGGCGCAGTGGGCCGAACAGCTCATTGCGCTTGGCTTTGCCCGGCAGGCGCCAGCTG  >  NC_002947/3099349‑3099640

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: