Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,256,493 A→C 28.1% L178R (CTG→CGG)  PP_3726 ← enoyl‑CoA hydratase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,256,4930AC28.1% 26.4 / 6.1 38L178R (CTG→CGG) PP_3726enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (9/15);  new base C (8/3);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 7.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGG  >  NC_002947/4256367‑4256619
                                                                                                                              |                                                                                                                              
gcgcgTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGGCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAt                                                                                                                       <  1:123158/135‑1 (MQ=255)
             gCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGGCCGCTCGGCGgccccccc                                                                                                          >  1:304850/1‑132 (MQ=255)
              tttGGTCTGCGCATAGGCGATGGTGGGGCGTTTGGCCTGGCGCCCGGCCTGTTCGGGGGCAACGGCAGGCATCTGCTGGTCGGCGCAGTCTTCGCTGACCAACCCGGCGGCCTGGGCGCTTTCGGCGCCCCACAAc                                                                                                         <  2:549462/135‑1 (MQ=255)
              tttGGTCTGCGCATAGGCGAAGGTCGGGCCTTTGGCCAGGCGCCGGGCCTGTTCGGCGGCAACGGCTGGCAACTGCGCGTGGGCGCTGAGTTCGGTGACCAACCCGGCGGCCAGGGCGCGAGCGGGGCCCCACAAc                                                                                                         <  2:628654/135‑1 (MQ=255)
                          aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAc                                                                                             >  2:37917/1‑136 (MQ=255)
                          aTAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCACGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAc                                                                                             >  2:352986/1‑136 (MQ=255)
                           tAGGCGAAGGTCGGGCCTTTGGCCTGGCGCCCGGCCTGTTCGGCGGCAACGGCAGGCAACTGCGCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCTGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAca                                                                                            <  1:556430/136‑1 (MQ=255)
                               cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAcagcag                                                                                        >  2:495903/1‑136 (MQ=255)
                               cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGATCGGCGCCCCCCAACTCGGCGAGGAAcagccc                                                                                        >  1:601158/1‑134 (MQ=255)
                                gAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAgg                                                                                       >  2:510555/1‑136 (MQ=255)
                                  tgggcGGGCCATTGGCCTGGCTCCCGGCCTGTTCGGCGCAAAGGGCAGGCAAGTGCTCGTCGGCGCTTACTTCGCTGACCAACCGGGCGGCCAGGGCGCGATGGGCGCCCCACAACTCGTCGAGGAACAGCAGgcg                                                                                     <  1:484988/132‑1 (MQ=255)
                                    gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCAACGCGTCGGGGGACCGCAGgcgcg                                                                                   >  1:118198/1‑135 (MQ=255)
                                    gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCGATAGGCGCCCCCCCACACGAACAGGAACAACAGgcgca                                                                                   >  1:410290/1‑135 (MQ=255)
                                    gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCCCCACTCGTCGAGGAACAGCAGGCGCt                                                                                   >  1:176839/1‑136 (MQ=255)
                                    gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCt                                                                                   >  1:438275/1‑136 (MQ=255)
                                           cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTg                                                                            >  1:628874/1‑136 (MQ=255)
                                                    ggCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                   <  1:495903/136‑1 (MQ=255)
                                                    ggCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                   <  1:271895/136‑1 (MQ=255)
                                                                        gcAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGcc                                               <  2:262214/136‑1 (MQ=255)
                                                                         cAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAgg                                                                                       >  2:484105/1‑95 (MQ=255)
                                                                         cAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAgg                                                                                       <  1:484105/95‑1 (MQ=255)
                                                                             ggcaggcaACTGCTCGTCGGCGCAGACTACGCTGACCAACCCGGCGGCCTGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTg                                          <  2:588809/136‑1 (MQ=255)
                                                                              gcaggcaACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTgg                                         >  2:425776/1‑136 (MQ=255)
                                                                                      cTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                 <  1:169749/136‑1 (MQ=255)
                                                                                      cTGCTCGGCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                 <  1:320464/136‑1 (MQ=255)
                                                                                      cTGCTCGGCGGCGCAGACTGCGCTCCCCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                 <  2:257519/136‑1 (MQ=255)
                                                                                      cTGCTCGGCGCCGCAGACTGCGCTCCCCACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCACCGGCa                                 <  2:304850/136‑1 (MQ=255)
                                                                                       tGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCAt                                <  2:410290/136‑1 (MQ=255)
                                                                                           gggcgGCGCAGACTGCGCTCACCACCCCGGCGGCCCGGGCGCGAGCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                            <  2:536241/133‑1 (MQ=255)
                                                                                           gggcgGCGCACCCTGCGCTCACCACCCGGGCGGCCCGGGCGCGACCGCCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                            <  1:646600/133‑1 (MQ=255)
                                                                                            gtcgGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                           <  2:438275/136‑1 (MQ=255)
                                                                                            ggcgccgcAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                           <  2:249011/131‑1 (MQ=255)
                                                                                             gcgGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGc                          <  2:662629/135‑1 (MQ=255)
                                                                                                 cgcAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACACCTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCAGtt                      >  2:660833/1‑135 (MQ=255)
                                                                                                         tCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACACGTCGAGGAACACCAGGCGCTTGGCCTGCTCGCTGGCGATCAGCCCTGGCAGGTTCCCGCAGGCGCCGGCCTCGGGGCTGTAGccccagc              >  2:646600/1‑132 (MQ=255)
                                                                                                             tGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCCGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTggg          >  1:136932/1‑135 (MQ=255)
                                                                                                               aCCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGTTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt        <  1:426797/136‑1 (MQ=255)
                                                                                                                     ccGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCgg  >  1:42774/1‑136 (MQ=255)
                                                                                                                              |                                                                                                                              
TCGCGTATCAGCCGCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGG  >  NC_002947/4256367‑4256619

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