Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,997 A→G 9.3% L403P (CTG→CCG)  PP_4635 ← trans‑2‑enoyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,257,9970AG 9.3% 113.7 / 3.1 43L403P (CTG→CCG) PP_4635trans‑2‑enoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base G (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

ATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTAGTCAGTCAGCTCGAACAGGTTCTCGGTGGTCA  >  NC_002947/5257848‑5258135
                                                                                                                                                     |                                                                                                                                          
aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACa                                                                                                                                            >  2:59525/1‑150 (MQ=255)
aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACa                                                                                                                                            >  1:163728/1‑150 (MQ=255)
   gcgcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCt                                                                                                                                         >  1:114131/1‑150 (MQ=255)
                  ttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACt                                                                                                                          <  2:114131/150‑1 (MQ=255)
                          cAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCg                                                                                                                  >  2:435279/1‑150 (MQ=255)
                                   cAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACg                                                                                                         >  1:220858/1‑150 (MQ=255)
                                    aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGt                                                                                                        >  1:1123806/1‑150 (MQ=255)
                                       gATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCtt                                                                                                     >  2:948175/1‑150 (MQ=255)
                                              tAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTa                                                                                              <  1:948175/150‑1 (MQ=255)
                                                  cTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCa                                                                                          >  2:526265/1‑150 (MQ=255)
                                                   tGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCaa                                                                                         >  1:1144217/1‑150 (MQ=255)
                                                    gTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAAc                                                                                        <  2:163728/150‑1 (MQ=255)
                                                    gTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAAc                                                                                        <  2:1124710/150‑1 (MQ=255)
                                                        tGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCa                                                                                    >  1:368360/1‑150 (MQ=255)
                                                        tGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCa                                                                                    >  2:1221350/1‑150 (MQ=255)
                                                          gCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACGGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgc                                                                                  >  2:95182/1‑150 (MQ=255)
                                                          gCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACGGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgc                                                                                  >  2:95882/1‑150 (MQ=255)
                                                             tctTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACACGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCGACGTCAgcgcg                                                                               >  1:987505/1‑150 (MQ=255)
                                                               ttCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCGGTCGAACTTCACATCTGAGTCCACGTCTTTGTCGTAGTCACAGTCAGCGCGCt                                                                             >  2:102881/1‑150 (MQ=39)
                                                               ttCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTc                                                                                                 >  1:1168296/1‑130 (MQ=255)
                                                               ttCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTc                                                                                                 <  2:1168296/130‑1 (MQ=255)
                                                               ttCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCt                                                                             >  1:252284/1‑150 (MQ=255)
                                                                   cgGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGaa                                                                         <  1:59525/150‑1 (MQ=255)
                                                                         aaCCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGaa                                                                   <  1:526265/150‑1 (MQ=255)
                                                                           cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGTTCACCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACa                                                                 >  1:697081/1‑150 (MQ=255)
                                                                           cccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACa                                                                 >  2:372794/1‑150 (MQ=255)
                                                                            ccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACGGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAg                                                                <  1:95182/150‑1 (MQ=255)
                                                                            ccGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACGGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAg                                                                <  1:95882/150‑1 (MQ=255)
                                                                               cTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGtt                                                             <  2:1123806/150‑1 (MQ=255)
                                                                                      aGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCg                                                                                                <  1:832386/108‑1 (MQ=255)
                                                                                      aGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCg                                                                                                >  2:832386/1‑108 (MQ=255)
                                                                                          ggAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGct                                                  <  1:435279/150‑1 (MQ=255)
                                                                                            tcccccTGTAGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGATTGAGGAACTGcttc                                                <  1:458891/146‑1 (MQ=255)
                                                                                            agaACCTGTGGGAGGGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACCTCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttc                                                <  1:314485/150‑1 (MQ=255)
                                                                                            agaACCTGTGGGAGCGGGTTCACCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttc                                                <  2:697081/150‑1 (MQ=255)
                                                                                              aaCCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGcttctt                                              <  2:777306/150‑1 (MQ=255)
                                                                                                      ggAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGAGTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCgg                                      >  2:930943/1‑150 (MQ=255)
                                                                                                       gAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGc                                     <  2:220858/150‑1 (MQ=255)
                                                                                                       gAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGc                                     >  1:849729/1‑150 (MQ=255)
                                                                                                          cGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTa                                  <  1:372794/150‑1 (MQ=255)
                                                                                                              cTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTagtca                              <  1:756004/150‑1 (MQ=255)
                                                                                                                  cccGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGt                                                                                                                >  2:422911/1‑64 (MQ=255)
                                                                                                                  cccGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGt                                                                                                                <  1:422911/64‑1 (MQ=255)
                                                                                                                           agGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCCAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTAGTCAGTCAGCTCGAACa                 <  2:47528/150‑1 (MQ=255)
                                                                                                                                          tACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTAGTCAGTCAGCTCGAACAGGTTCTCGGTGGTCa  <  1:1221350/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
ATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGGAGAACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCGAACAGGTTGAGGAACTGCTTCTTGTAACCGGCGTAGTCAGTCAGCTCGAACAGGTTCTCGGTGGTCA  >  NC_002947/5257848‑5258135

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