Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,591,702 G→A R258H (CGC→CAC)  liuC → methylglutaconyl‑CoA hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,591,7020GA96.4% 85.4 / NA 28R258H (CGC→CAC) liuCmethylglutaconyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/18);  minor base C (1/0);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACCAAGGACTTGCTGCGCGAAGTGGACGACGGCGAACTCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCGCCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATCGCCTGCCGGGTGATGCGCACCGCCAAGGCCATGGGCCTGACCACT  >  NC_002947/4591563‑4591837
                                                                                                                                           |                                                                                                                                       
gcCACCAAGGACTTGCTGCGCGAAGTGGACGACGGCGAACTCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACcctgcct                                                                                                                                 <  5:93361/148‑1 (MQ=255)
      aaGGACTTGCTGCGCGAAGTGGACGACGGCGAACTCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCaaa                                                                                                                           >  3:87909/1‑148 (MQ=255)
           cTTGCTGCGCGAAGTGGACGACGGCGAACTCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGtt                                                                                                                      <  3:35023/148‑1 (MQ=255)
               ctgcGCGAAGTGGACGACGGCGAACTCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGAc                                                                                                                   <  4:83572/147‑1 (MQ=255)
                   gcgAAGTGGACGACGGCGAACTCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACaagaa                                                                                                              <  4:110257/148‑1 (MQ=255)
                   gcgAAGGGGACGACGGCGAACTCAGACCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACaagaa                                                                                                              <  6:30889/148‑1 (MQ=255)
                                   cGAACTCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGa                                                                                              <  2:5628/148‑1 (MQ=255)
                                        tCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAgccgc                                                                                         >  6:20729/1‑148 (MQ=255)
                                            cccGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTGTCCTGGAAAAACGCCCCCCGGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCg                                                                                     >  5:91765/1‑148 (MQ=255)
                                                gCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCttt                                                                                 <  1:169911/148‑1 (MQ=255)
                                                  ccTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGa                                                                               <  4:161002/148‑1 (MQ=255)
                                                   cTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGac                                                                              <  1:155514/148‑1 (MQ=255)
                                                   cTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGa                                                                               >  1:82767/1‑147 (MQ=255)
                                                           ttACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGAGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCctgc                                                                      <  3:116314/148‑1 (MQ=255)
                                                             aCTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCctgctg                                                                    <  4:124189/148‑1 (MQ=255)
                                                                  gAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTc                                                                 <  2:175383/146‑1 (MQ=255)
                                                                  gAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCg                                                                >  6:105396/1‑147 (MQ=255)
                                                                      acacCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAAc                                                           <  3:101239/148‑1 (MQ=255)
                                                                               cccGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGaa                                                  >  5:43035/1‑148 (MQ=255)
                                                                                  gTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATc                                               <  7:115223/148‑1 (MQ=255)
                                                                                          gTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATCGCCTGCCg                                       <  2:72305/148‑1 (MQ=255)
                                                                                             aGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAACCGCCACCCTGCCTGGCCAACCGGTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGCAATCGCCTGCCggg                                     >  7:68311/1‑147 (MQ=255)
                                                                                                                      gcgcTTTCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATCGCCTGCCGGGTGATGCGCACCGCCAAGGCCATgg             <  5:105396/146‑1 (MQ=255)
                                                                                                                       cgcTTTCCTGGAAAAACGCCACCGTGCCTTGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATCGCCTGCCGGGTGATGCGCACCGCCAAGGCCATGGGcc          <  4:87909/148‑1 (MQ=255)
                                                                                                                          tttCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATCGCCTGCCGGGTGATGCGCACCGCCAAGGCCATGGGCCTGa       >  1:151876/1‑148 (MQ=255)
                                                                                                                          tttCCTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTCGTCGCCAACCGCGGCGAACTCGCCTGCCGGGTGATGCGCACCGCCAAGGCCATGGGCCTGa       >  4:52607/1‑148 (MQ=255)
                                                                                                                             ccTGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATCGCCTGCCGGGTGATGCGCACCGCCAAGGCCATGGGCCTGAc      <  2:138146/146‑1 (MQ=255)
                                                                                                                               tGGAAAAACGCCACCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATCGCCTGCCGGGTGATGCGCACCGCCAAGGCCATGGGCCTGACCACt  >  5:97716/1‑148 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GCCACCAAGGACTTGCTGCGCGAAGTGGACGACGGCGAACTCAGCCCGGCCCTGCGCCGTTACTGCGAAAACACCATCGCCCGTATCCGCGTCAGCGCCGAAGGCCAGGAGGGCCTGCGCGCTTTCCTGGAAAAACGCCGCCCTGCCTGGCAAACCGTTGACAAGAAGGAGCCGCGCCCATGAGCCGCCCCGCTTTGACCACCCTGCTGGTCGCCAACCGCGGCGAAATCGCCTGCCGGGTGATGCGCACCGCCAAGGCCATGGGCCTGACCACT  >  NC_002947/4591563‑4591837

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