Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,795,913 A→T 18.7% V192D (GTC→GAC)  PP_3355 ← beta‑ketothiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,795,9130AT18.7% 45.4 / 6.7 27V192D (GTC→GAC) PP_3355beta‑ketothiolase
Reads supporting (aligned to +/- strand):  ref base A (14/8);  new base T (3/2);  total (17/10)
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.04e-01

GCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATAGGCACGGGCCAGGTTTTCCGCAGTACC  >  NC_002947/3795782‑3796037
                                                                                                                                   |                                                                                                                            
gcAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGAt                                                                                                                          >  2:398515/1‑136 (MQ=255)
          cAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcgtcg                                                                                                                <  2:446289/136‑1 (MQ=255)
                                ggtggggggATGTCCTGGCTGACCTTACTTTCGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGTCGCTGACGTTTTCCTCGTCGAACCAGCCCTGCACCTGGGCtt                                                                                          <  2:282118/135‑1 (MQ=255)
                                    ggtggCGGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAg                                                                                      <  1:251793/136‑1 (MQ=255)
                                     gtggCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAgg                                                                                     >  2:308087/1‑136 (MQ=255)
                                           gTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCAc                                                                                                  >  1:18220/1‑117 (MQ=255)
                                           gTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCAc                                                                                                  <  2:18220/117‑1 (MQ=255)
                                            tCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCt                                                                              >  2:337449/1‑136 (MQ=255)
                                                     aCCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTgc                                                                     <  1:308087/136‑1 (MQ=255)
                                                         cactcacGCCACGCGGCAATGCAAAGCCGCGTGGCTGGCAGCCTTCTACCTCCACTGCCTGCGCGGTGACGCTGACGACTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGGCCGCTGATGGCAGCGCac                                                                 >  2:305735/1‑135 (MQ=255)
                                                                    cgcgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGt                                                      >  2:13762/1‑136 (MQ=255)
                                                                       ggCCATGCAAAGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCACCCCAGCACCTGGGCTTGCCCGGCCCGCTGATGGCTGCGCAGGGCCCCCCCGTccc                                                   >  2:352540/1‑135 (MQ=255)
                                                                       ggCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTccc                                                   >  2:361731/1‑135 (MQ=255)
                                                                       ggCAATGCAATGCCGCGTGGCAGGCAGCCTTCTACCTGCAGTGCCTGCGCCGGGACGCTGTCGATTTCCTCGTCGAACCAGCCCAGCACCGGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCAcgcctgcc                                                   >  2:409208/1‑131 (MQ=255)
                                                                                 ggcCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg                                         <  1:154248/135‑1 (MQ=255)
                                                                                    cgcgTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACgcgc                                      >  2:260184/1‑136 (MQ=255)
                                                                                         ggctggcAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGc                                 >  1:349146/1‑136 (MQ=255)
                                                                                         ggctggcAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGc                                 >  2:43738/1‑136 (MQ=255)
                                                                                           ctggcAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGTCGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCt                                                                       >  1:386872/1‑96 (MQ=255)
                                                                                           ctggcAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGTCGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCt                                                                       <  2:386872/96‑1 (MQ=255)
                                                                                           ctggcAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGTCGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCt                                                                       <  2:389253/96‑1 (MQ=255)
                                                                                           ctggcAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGTCGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCt                                                                       >  1:389253/1‑96 (MQ=255)
                                                                                                       tACCTCCAGTGCCTGCGCGGTGACGCTGACGAATTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATAGGCACGGGcc                   >  2:324490/1‑136 (MQ=255)
                                                                                                              aGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATAGGCACGGGCCAGGtttt            <  1:43738/136‑1 (MQ=255)
                                                                                                                 gCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcgc                                        >  1:441343/1‑105 (MQ=255)
                                                                                                                 gCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcgc                                        <  2:441343/105‑1 (MQ=255)
                                                                                                                        cgGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCCGGCCATAGGCACGGGCCAGGTTTTCCGCAGTAcc  >  1:232460/1‑136 (MQ=255)
                                                                                                                                   |                                                                                                                            
GCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGCGCCAGGCCATAGGCACGGGCCAGGTTTTCCGCAGTACC  >  NC_002947/3795782‑3796037

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: