Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,962,620 Δ1 bp 9.1% coding (774/1152 nt) PP_3492 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,962,6200A. 9.1% 129.1 / 7.2 44coding (774/1152 nt)PP_3492acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (18/22);  new base . (2/2);  total (20/24)
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.60e-01

CTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGCCGCCCGCCTGCGCA  >  NC_002947/3962476‑3962768
                                                                                                                                                |                                                                                                                                                      
cTGCGCAGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGa                                                                                                                                                   >  2:939547/1‑150 (MQ=255)
 ggcgcGGGCTCCTGCGACAACTGCCGCATTCCTGCCGCCAATATCCCCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGaa                                                                                                                                                  <  2:768252/149‑1 (MQ=255)
  gcgcGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGCCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑g                                                                                                                                                <  1:1044557/150‑1 (MQ=255)
    gcgGTCACCCTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑gcg                                                                                                                                              <  2:598564/150‑1 (MQ=255)
         cACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGC‑TTCGAATGCGgcgcg                                                                                                                                         <  2:985846/150‑2 (MQ=255)
                 aCAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGtgatgc                                                                                                                                 >  2:855503/1‑145 (MQ=255)
                         cGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGt                                                                                                                                      <  1:64862/137‑1 (MQ=255)
                         cGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGt                                                                                                                                      >  2:64862/1‑137 (MQ=255)
                                   ccgccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCTCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGt                                                                                                              >  2:235918/1‑150 (MQ=255)
                                    cgccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGtt                                                                                                             <  2:215020/150‑1 (MQ=255)
                                     gccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGttt                                                                                                            <  1:279193/150‑1 (MQ=255)
                                     gccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGttt                                                                                                            <  1:537370/150‑1 (MQ=255)
                                       cAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTgg                                                                                                          >  2:891585/1‑150 (MQ=255)
                                       cAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGC‑TTCGAATGCGgcgcgg                                                                                                                                        <  2:946114/121‑3 (MQ=255)
                                       cAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGC‑TTCGAATGCGgcgcgg                                                                                                                                        >  1:946114/1‑119 (MQ=255)
                                            tCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc                                                                                                     <  2:279050/150‑1 (MQ=255)
                                                    gAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGGCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAAc                                                                                             <  1:446304/150‑1 (MQ=255)
                                                          ggCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAacac                                                                                       <  2:542919/150‑1 (MQ=255)
                                                              aGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCaga                                                                                   <  1:1036795/150‑1 (MQ=255)
                                                                     ggCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc                                                                                                     <  2:503196/125‑1 (MQ=255)
                                                                     ggCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc                                                                                                     >  1:503196/1‑125 (MQ=255)
                                                                        cATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCaa                                                                         <  1:939547/150‑1 (MQ=255)
                                                                                 gTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGc                                                                <  1:746948/150‑1 (MQ=255)
                                                                                    cAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGa                                                             >  2:736571/1‑150 (MQ=255)
                                                                                             aGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACAtgcaggt                                                    >  2:867589/1‑150 (MQ=255)
                                                                                                    cGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGcaggtgcaggtg                                             <  2:133414/150‑1 (MQ=255)
                                                                                                      tATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGaa                                           >  1:1127293/1‑150 (MQ=255)
                                                                                                        tCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc                                                                                                     >  1:869097/1‑90 (MQ=255)
                                                                                                        tCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc                                                                                                     <  2:869097/90‑1 (MQ=255)
                                                                                                        tCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATgcaggtgcaggt                                              <  1:1085632/145‑1 (MQ=255)
                                                                                                        tCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATgcaggtgcaggt                                              >  2:1085632/1‑145 (MQ=255)
                                                                                                             aTTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAAcgccgc                                     >  1:626708/1‑149 (MQ=255)
                                                                                                             aTTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAAcgccgc                                     <  2:626708/149‑1 (MQ=255)
                                                                                                             aTTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCa                                    >  2:1161316/1‑150 (MQ=255)
                                                                                                                gccgccCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGt                                 <  1:1161316/150‑1 (MQ=255)
                                                                                                                    cccAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGc                             <  1:736571/150‑1 (MQ=255)
                                                                                                                          ccTTGGGCATCGCCCGCGCCGCATTCGCA‑GCGGGCCTGGTATAATCGCGTTACTGCAG‑ACAGTTTGGCAAGCCGATCAACCACGACCATAGCATGGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCACACGTGGGCTGATcc                       >  1:931335/1‑150 (MQ=255)
                                                                                                                           cTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc                                                                                                     >  2:875462/1‑71 (MQ=255)
                                                                                                                           cTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGc                                                                                                     <  1:875462/71‑1 (MQ=255)
                                                                                                                            ttGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTg                     >  1:641361/1‑150 (MQ=255)
                                                                                                                                       cagcgcCGC‑TTCGAATGCGGCGCTGGTCTACTCGCGTGACCGCATCACCG‑TTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGACATACAGGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGccgcccg          >  1:874282/3‑150 (MQ=255)
                                                                                                                                               cATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCt                            <  2:402813/124‑1 (MQ=255)
                                                                                                                                               cATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCt                            >  1:402813/1‑124 (MQ=255)
                                                                                                                                               cATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGCCGCCCGCCTGCGCa  >  2:441436/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CTGCGCGGTCACCTTCGACAACTGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCAT‑ACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATTGCCAACCTGCTGGCCGACATGCAGGTGCAGGTGAACGCCGCACGTTTGCTGATCCTGCATGCCGCCCGCCTGCGCA  >  NC_002947/3962476‑3962768

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: ‑