Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,532 G→A 17.4% V224V (GTG→GTA allE → S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,002,5320GA17.4% 45.1 / 4.6 23V224V (GTG→GTAallES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (8/11);  new base A (2/2);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

ACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGAC  >  NC_002947/4002395‑4002678
                                                                                                                                         |                                                                                                                                                  
accCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAAc                                                                                                                                        >  1:300382/1‑150 (MQ=255)
 cccGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCtt                                                                                                                                           <  2:448601/146‑2 (MQ=255)
 cccGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCtt                                                                                                                                           >  1:448601/1‑145 (MQ=255)
   cGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAg                                                                                                                                     <  2:509300/150‑1 (MQ=255)
                     gACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGTCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGaa                                                                                                                   <  1:359667/150‑1 (MQ=255)
                            gcCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCg                                                                                                            >  1:377179/1‑150 (MQ=255)
                                    aTGCACGTGAACATCGTCACCTTCCATCCTGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCAtg                                                                                                    >  1:68765/1‑150 (MQ=255)
                                        aCGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCtt                                                                                                                                           >  1:425267/1‑106 (MQ=255)
                                        aCGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCtt                                                                                                                                           <  2:425267/107‑2 (MQ=255)
                                                            cAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCaa                                                                            <  2:130050/150‑1 (MQ=255)
                                                                   gCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCt                                                                     <  1:226359/150‑1 (MQ=255)
                                                                            ttCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcggcg                                                            >  2:34037/1‑150 (MQ=255)
                                                                                     cAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCATGATACGCCGGCGGCCCAGGCa                                                   >  2:93219/1‑150 (MQ=255)
                                                                                           cccACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTa                                             <  2:377179/150‑1 (MQ=255)
                                                                                              aCGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcggcg                                                            >  2:17946/1‑132 (MQ=255)
                                                                                              aCGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcggcg                                                            <  1:17946/132‑1 (MQ=255)
                                                                                                           cGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAGGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACac                             <  1:93219/150‑2 (MQ=255)
                                                                                                              ccTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGa                          <  1:34037/150‑1 (MQ=255)
                                                                                                               cTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcgg                                                              >  1:64473/1‑113 (MQ=255)
                                                                                                               cTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcgg                                                              <  2:64473/113‑1 (MQ=255)
                                                                                                                     gTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAAc                   >  2:341625/1‑150 (MQ=255)
                                                                                                                           gAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCAc             <  1:341625/150‑1 (MQ=255)
                                                                                                                                      ggTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGCGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGAc  <  2:68765/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
ACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGAC  >  NC_002947/4002395‑4002678

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