Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,241,585 A→C V96G (GTG→GGG)  hmgC ← maleylacetoacetate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,241,5850AC96.4% 92.6 / ‑3.3 28coding (287/633 nt)hmgCmaleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (15/12);  minor base . (0/1);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGC  >  NC_002947/5241447‑5241728
                                                                                                                                          |                                                                                                                                                
tCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGaa                                                                                                                                         >  4:156837/1‑148 (MQ=255)
        cgcgAGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGcc                                                                                                                                 >  4:82993/1‑148 (MQ=255)
             gCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGGGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGc                                                                                                                            <  7:137626/148‑1 (MQ=255)
              cccTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGcc                                                                                                                           >  4:170106/1‑148 (MQ=255)
               ccTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCt                                                                                                                          >  7:122235/1‑148 (MQ=255)
                 ttGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCCGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTcg                                                                                                                        >  3:171438/1‑148 (MQ=255)
                    gCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGGTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTcgg                                                                                                                       <  1:62314/146‑1 (MQ=255)
                     cTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTcggcgg                                                                                                                    >  1:20024/1‑148 (MQ=255)
                                      ccACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCa                                                                                                   >  8:11772/1‑148 (MQ=255)
                                        aCTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGACCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAgg                                                                                                 >  8:131086/1‑148 (MQ=255)
                                                      cccTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAAc                                                                                     <  4:178508/146‑1 (MQ=255)
                                                      cccTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACct                                                                                   >  4:152948/1‑148 (MQ=255)
                                                        cTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGCCCCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCCCCGCAGGCTGTGGATAAACctct                                                                                 >  1:174354/1‑148 (MQ=255)
                                                         tCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCtt                                                                                >  1:16282/1‑148 (MQ=255)
                                                            tCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTcc                                                                              <  8:153269/147‑1 (MQ=255)
                                                            tCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCACCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCATTCTGTGGATAAACCTCTTc                                                                               <  6:16617/146‑1 (MQ=255)
                                                                     gCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCACCCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCga                                                                    >  7:133113/1‑148 (MQ=255)
                                                                                       aGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCGTCCAGGTACTCGATGATCGCCGGCGCCTGcc                                                  >  6:104299/1‑147 (MQ=255)
                                                                                        gCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCAc                                                 <  5:104299/148‑1 (MQ=255)
                                                                                        gCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCAc                                                 <  6:47532/147‑1 (MQ=255)
                                                                                         cacaCTGCCGTTGGGCAGCGGATTAATGGCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCAcc                                                <  7:15935/148‑1 (MQ=255)
                                                                                          acacTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCa                                               >  3:76504/1‑148 (MQ=255)
                                                                                           cacTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAg                                              <  8:122235/148‑1 (MQ=255)
                                                                                            acTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGc                                             <  3:156837/148‑1 (MQ=255)
                                                                                                 cGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCAcc                                          <  2:16282/146‑1 (MQ=255)
                                                                                                         aGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGt                                <  4:76504/148‑1 (MQ=255)
                                                                                                                                      cgccCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAGGGCCGCGTACCGGCCTTGc  <  3:82993/148‑1 (MQ=255)
                                                                                                                                      cgccCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCGGCCAGGTCCTCGATGAGCGCCGGCGCCTGCCCCAGCACCTCCCCACCGTCGGTGCGCAAGGCCGGTCCGCGGCCCTGc  >  1:58113/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                
TCCACTGCCGCGAGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGGCAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGC  >  NC_002947/5241447‑5241728

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