Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 353,067 A→G 100% intergenic (+34/+84) hisF → / ← cbcV Imidazole glycerol phosphate synthase subunit HisF/choline / betaine / carnitine ABC transporter ‑ ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451353,0670AG100.0% 81.6 / NA 24intergenic (+34/+84)hisF/cbcVImidazole glycerol phosphate synthase subunit HisF/choline / betaine / carnitine ABC transporter ‑ ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/6);  total (18/6)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCG  >  AE015451/352860‑353283
                                                                                                                                                                                                               |                                                                                                                                                                                                                        
cTGGGGATTCCGGTGATTGCCGCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                                                                                                              <  1:57371/236‑1 (MQ=255)
  ggggATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGAGCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGccc                                                                                                                                                                                                                 >  2:865575/1‑215 (MQ=255)
                                              aCCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGaa                                                                                                                                                      >  1:920670/1‑230 (MQ=255)
                                                          gCATTCTGGAAGGGCATGCCAGTGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGccc                                                                                                                                                                     <  2:750998/203‑1 (MQ=255)
                                                          gCATTCTGGAAGGGCATGCCAGTGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGccc                                                                                                                                                                     >  1:751001/1‑203 (MQ=255)
                                                                  gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                                                                                                              <  2:778975/170‑1 (MQ=255)
                                                                             cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCGCCCGCTCCCACAGGt                                                                                                                                                                                                     >  2:102527/1‑152 (MQ=255)
                                                                                   ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcac                                                                                                                                                                                                         >  2:379044/1‑142 (MQ=255)
                                                                                                  cATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgc                                                                                                                                                         >  2:934259/1‑175 (MQ=255)
                                                                                                    tCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCcgcg                                                                                                                                                        <  2:815757/174‑1 (MQ=255)
                                                                                                                gCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATccc                                                                              >  1:189262/1‑236 (MQ=255)
                                                                                                                 cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCGTCAACCGTGGTTCTTGCCGa                                                                                                                   >  2:747143/1‑198 (MQ=255)
                                                                                                                    gTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCcca                                                                            >  2:673599/1‑234 (MQ=255)
                                                                                                                                  aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                                                                                                              >  1:730883/1‑186 (MQ=255)
                                                                                                                                       aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTCCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTAtg                                                            >  1:200400/1‑230 (MQ=255)
                                                                                                                                         ggCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTggg                                                                                                                                                                                  >  2:685461/1‑111 (MQ=255)
                                                                                                                                         ggCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGt                                                                                                                               >  1:331302/1‑162 (MQ=255)
                                                                                                                                         ggCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTaa                                                                                                     >  2:372163/1‑188 (MQ=255)
                                                                                                                                         ggCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTAATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa                                                                                     >  2:591508/1‑204 (MQ=255)
                                                                                                                                                                                         gggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATg     <  1:673602/236‑1 (MQ=255)
                                                                                                                                                                                            ctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGt                                                                   >  2:838437/1‑171 (MQ=255)
                                                                                                                                                                                            ctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTg                     >  2:864370/1‑217 (MQ=255)
                                                                                                                                                                                             tgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCg  <  1:934262/235‑1 (MQ=255)
                                                                                                                                                                                               ctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGGGATCTCtg                                                                                                                                                                                      >  1:237556/1‑53 (MQ=17)
                                                                                                                                                                                                               |                                                                                                                                                                                                                        
CTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCG  >  AE015451/352860‑353283

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