Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,241,585 A→C 100% 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,5850AC94.4% 135.9 / ‑1.0 46V96G (GTG→GGG) hmgCmaleylacetoacetate isomerase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (21/22);  total (23/23)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

AGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGG  >  NC_002947/5241459‑5241730
                                                                                                                               |                                                                                                                                                  
aGCCCTTGGCTGATCCAGTGGCCGAGCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGc                                                                                                                              <  1:224192/149‑1 (MQ=255)
     ttGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGACGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCACCCCCCCAACCCTGGCCCGTGCCTcgg                                                                                                                         >  7:149296/1‑149 (MQ=255)
       ggCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTcgg                                                                                                                         <  2:18904/147‑1 (MQ=255)
            aTCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTgg                                                                                                                   >  3:293250/1‑149 (MQ=255)
            aTCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGc                                                                                                                  >  7:21561/1‑149 (MQ=255)
                       gATCCACTGGTTGACCTGCCCCTCCTCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAAc                                                                                                       >  6:320906/1‑148 (MQ=255)
                       gATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACg                                                                                                      >  3:30976/1‑149 (MQ=255)
                           cACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCGCCCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCa                                                                                                    <  4:293250/148‑1 (MQ=255)
                            aCTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCACGGCGGTGGCCGG‑CAGCAACGCAGGc                                                                                                 >  8:86089/1‑149 (MQ=255)
                                gTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGg‑ca                                                                                                            >  8:188583/1‑134 (MQ=255)
                                gTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGg‑ca                                                                                                            <  7:188583/134‑1 (MQ=255)
                                gTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCtgtg                                                                                             <  4:30976/149‑1 (MQ=255)
                                         ccccTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAAc                                                                                      <  3:309819/147‑1 (MQ=255)
                                          cccTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACctc                                                                                   <  8:21561/149‑1 (MQ=255)
                                            cTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAACCCTCtt                                                                                 >  1:125059/1‑149 (MQ=255)
                                               aTCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTcc                                                                               <  5:296658/148‑1 (MQ=255)
                                                tCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAg                                                                             <  6:129175/149‑1 (MQ=255)
                                                   tGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTa                                                                          >  8:235793/1‑149 (MQ=255)
                                                    gACCGGCCTGGCGCAGCCGTTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGGTACGAATTATGGCCG‑CCCCCCCCCGATCGCTGGCCGGTGGCTCGGCGGTGGCCGGCCAGCAACGCAGGCTGTGGATAAACCTCTTCCAg                                                                             <  1:38125/146‑1 (MQ=255)
                                                    gACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑ACCCCCCCCCAAACCTGGCCCGTGCCTCGGCGGTGGACGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAg                                                                             >  2:38125/1‑145 (MQ=255)
                                                     aCCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCCGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTa                                                                          <  8:52338/147‑1 (MQ=255)
                                                         gCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCgat                                                                    <  5:320906/149‑1 (MQ=255)
                                                              gcgcctccGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATc                                                                <  7:268116/142‑1 (MQ=255)
                                                                  aGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCgg                                                           <  7:235793/149‑1 (MQ=255)
                                                                            gCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCAc                                                  <  7:86089/148‑1 (MQ=255)
                                                                             cacaCTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGc                                                          <  8:70835/139‑1 (MQ=255)
                                                                             cacaCTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGc                                                          >  7:70835/1‑139 (MQ=255)
                                                                                acTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTAc                                                                         <  8:244418/121‑1 (MQ=255)
                                                                                acTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTAc                                                                         >  7:244418/1‑121 (MQ=255)
                                                                                  tGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGTTc                                                                <  2:60355/128‑1 (MQ=255)
                                                                                  tGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGGTc                                                                >  1:60355/1‑128 (MQ=255)
                                                                                   gACGTTGTGCATCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCa                                                   <  3:209563/140‑1 (MQ=255)
                                                                                   gACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACAAGATGATCGCCGGCGACAGCa                                                   >  4:209563/1‑140 (MQ=255)
                                                                                      gTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTcac                                       >  5:307177/1‑149 (MQ=255)
                                                                                           gCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCg                                  >  6:14441/1‑149 (MQ=255)
                                                                                            cAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGt                                 >  4:278606/1‑149 (MQ=255)
                                                                                                    gAATGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTgcg                          >  3:180178/1‑148 (MQ=255)
                                                                                                      aTGTCGCAACCAATGATGGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCaa                       <  4:180178/149‑1 (MQ=255)
                                                                                                               ccAATGATGGCCG‑CCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTAc              >  6:173946/1‑149 (MQ=255)
                                                                                                               ccAAAGATGGCCG‑CCCCCCCCCGAACCCAGGCCCCTGCAACGCCCGTGGCCGG‑CAGCAACGCAGGCAGATGAAAAACCTAATCCAGCAACACGATGAACGCCCGCGACAGCCCCAGCACCTCACCACCCTACGTGCGCAACGCCGCTAc              >  6:185369/1‑149 (MQ=255)
                                                                                                                     atgGCCG‑CCCCCCCCCGAACCCTGGCCCTTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGAAAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGcc        >  7:233035/1‑149 (MQ=255)
                                                                                                                     atgGCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGcc        >  7:237396/1‑149 (MQ=255)
                                                                                                                     atgGCCG‑CCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGcc        <  5:173946/149‑1 (MQ=255)
                                                                                                                       ggCCG‑CCCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCt      >  4:40324/1‑149 (MQ=255)
                                                                                                                           g‑ccCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCGGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCgg  <  1:273756/149‑1 (MQ=255)
                                                                                                                           g‑ccCCCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCgg  <  4:305307/149‑1 (MQ=255)
                                                                                                                               |                                                                                                                                                  
AGCCCTTGGCTGATCCAGTGGCCGATCCACTGGTTGACCTGCCCCTCATCCTGACCGGCCTGGCGCAGCCGGTTGAGCACACTGACGTTGTGCAGCGGATGAATGTCGCAACCAATGATGGCCG‑CCACCCCCCGAACCCTGGCCCGTGCCTCGGCGGTGGCCGG‑CAGCAACGCAGGCTGTGGATAAACCTCTTCCAGGTACTCGATGATCGCCGGCGACTGCACCAGCACCTCACCACCGTCGGTGCGCAAGGCCGGTACGCGGCCCTGCGG  >  NC_002947/5241459‑5241730

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