Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,284,838 C→T 11.1% A329V (GCG→GTG)  fadH → 2,4‑dienoyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,284,8380CT11.1% 103.9 / 3.9 36A329V (GCG→GTG) fadH2,4‑dienoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base T (2/2);  total (20/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 = 7.12e-01

CAAGTTGCGCGGCGTGGTGAGCATTCCGCTGATCACCACCAACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTGGACCATACCTTCGGCGGCAAGCTGACCAGT  >  NC_002947/2284728‑2284965
                                                                                                              |                                                                                                                               
cAAGTTGCGCGGCGTGGTGAGCATTCCGCTGATCACCACCAACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGAc                                                                                                        <  2:338665/136‑1 (MQ=255)
      gcgcGGCGTGGTGAGCATTCCGCTGATCACCACCAACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTc                                                                                                  <  2:242364/136‑1 (MQ=255)
          ggcgTGGTGAGCATTCCGCTGATCACCACCAACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTcaaca                                                                                              <  1:316971/136‑1 (MQ=255)
               ggtgAGCATTCCGCTGATCACCACCAACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGcc                                                                                         <  2:10700/136‑1 (MQ=255)
                           gCTGATCACCACCAACCGCATCAACGCCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGt                                                                             >  1:301193/1‑136 (MQ=255)
                           gCTGATCACCACCAACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGt                                                                             >  2:284628/1‑136 (MQ=255)
                                tCACCACCAACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGa                                                                        >  2:14552/1‑136 (MQ=255)
                                       caACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATc                                                                 >  2:355083/1‑136 (MQ=255)
                                              aTCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCt                                                          <  2:317869/136‑1 (MQ=255)
                                               tCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTg                                                         <  1:405560/136‑1 (MQ=255)
                                               acAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTg                                                         <  1:337074/135‑1 (MQ=255)
                                                 aaCACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCa                                                       >  2:21493/1‑136 (MQ=255)
                                                    acCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCg                                                    <  1:14552/136‑1 (MQ=255)
                                                      cccGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGc                                                  >  2:200516/1‑136 (MQ=255)
                                                        cGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTg                                                >  1:36541/1‑136 (MQ=255)
                                                         ggAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGc                                               <  1:284628/136‑1 (MQ=255)
                                                           aaGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCaa                                             >  1:396736/1‑136 (MQ=255)
                                                           aaGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCaa                                             >  2:360793/1‑136 (MQ=255)
                                                           aaGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCaa                                             >  2:15750/1‑136 (MQ=255)
                                                           aaGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCaa                                             >  1:398709/1‑136 (MQ=255)
                                                             gTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAAcc                                           <  1:355083/136‑1 (MQ=255)
                                                             gTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAAcc                                           <  1:219733/136‑1 (MQ=255)
                                                                  cGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGgcc                                      <  1:155190/136‑1 (MQ=255)
                                                                   gAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGCCTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGgccc                                     >  1:140217/1‑135 (MQ=255)
                                                                   gAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGCAATCACCCCCTGCATCGGCTGCAACCAGgccg                                     >  1:72403/1‑135 (MQ=255)
                                                                   gAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGgcct                                     >  1:367804/1‑136 (MQ=255)
                                                                   gAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGgcct                                     >  1:398583/1‑136 (MQ=255)
                                                                        aGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGcctgcctg                                >  1:386099/1‑136 (MQ=255)
                                                                                  gAGGGCGATGCGGACATGGTCTCGATGGTGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTgg                               <  2:397629/127‑1 (MQ=255)
                                                                                  gAGGGCGATGCGGACATGGTCTCGATGGTGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTgg                               <  2:395665/127‑1 (MQ=255)
                                                                                  gAGGGCGATGCGGACATGGTCTCGATGGTGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTgg                               >  1:395665/1‑127 (MQ=255)
                                                                                  gAGGGCGATGCGGACATGGTCTCGATGGTGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTgg                               >  1:397629/1‑127 (MQ=255)
                                                                                  gAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTGGACCATACCt                      >  2:366808/1‑136 (MQ=255)
                                                                                       cGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTGGACCATACCTTcggc                 >  1:87378/1‑136 (MQ=255)
                                                                                          tGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTGGACCATACCTTCggcggc              <  1:366808/136‑1 (MQ=255)
                                                                                                      ctcGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTGGACCATACCTTCGGCGGCAAGCTGACCAGt  <  2:228302/136‑1 (MQ=255)
                                                                                                              |                                                                                                                               
CAAGTTGCGCGGCGTGGTGAGCATTCCGCTGATCACCACCAACCGCATCAACACCCCGGAAGTGGCCGAGGCAGTGCTGGCCGAGGGCGATGCGGACATGGTCTCGATGGCGCGACCGTTTCTCGCCGACCCGGACTTCGTCAACAAGGCCGCTGCCGGTCGTGCGGATGAAATCAACACCTGCATCGGCTGCAACCAGGCCTGCCTGGACCATACCTTCGGCGGCAAGCTGACCAGT  >  NC_002947/2284728‑2284965

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