Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,247,959 2 bp→CA coding (723‑724/813 nt) h16_A3001 → transcriptional regulator, LysR‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,247,9590AC100.0% 78.6 / NA 21G241G (GGA→GGCh16_A3001transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (9/12);  total (9/12)
*AM2604793,247,9600TA100.0% 76.8 / NA 21F242I (TTC→ATC) h16_A3001transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (9/12);  total (9/12)

TTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTG  >  AM260479/3247838‑3248080
                                                                                                                         ||                                                                                                                         
ttGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGc                                                                                                                >  1:725549/1‑134 (MQ=255)
                     aCAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGccgagcgc                                                                                           <  2:725549/134‑7 (MQ=255)
                                 cgcTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCtg                                                                            <  1:748132/134‑2 (MQ=255)
                                                         tCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTg                                                   >  2:536085/1‑134 (MQ=255)
                                                              cTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTgg                                              >  2:950432/1‑134 (MQ=255)
                                                                ggAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCg                                                                                                                  >  1:563033/1‑68 (MQ=255)
                                                                ggAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCg                                                                                                                  <  2:563033/68‑1 (MQ=39)
                                                                       aTCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAAcgt                                     <  2:298641/134‑3 (MQ=255)
                                                                          gACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGc                                                                                                           >  2:1070816/1‑65 (MQ=255)
                                                                          gACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGc                                                                                                           <  1:1070816/65‑1 (MQ=255)
                                                                                         aaGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTca                  <  2:586596/133‑1 (MQ=255)
                                                                                                    ggAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCa                                                       <  1:395520/87‑1 (MQ=255)
                                                                                                    ggAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCa                                                       >  2:395520/1‑87 (MQ=255)
                                                                                                     gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCAtc                                                                  <  2:988843/75‑1 (MQ=255)
                                                                                                     gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCAtc                                                                  >  1:988843/1‑75 (MQ=255)
                                                                                                     gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTca                  >  2:203785/1‑121 (MQ=255)
                                                                                                     gagaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTca                  <  1:203785/121‑1 (MQ=255)
                                                                                                       gaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTAcc    <  2:409425/134‑1 (MQ=255)
                                                                                                         tCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTg  <  1:536085/134‑1 (MQ=255)
                                                                                                              cTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTca                  >  2:6573/1‑112 (MQ=255)
                                                                                                              cTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTca                  <  1:6573/112‑1 (MQ=255)
                                                                                                                         ||                                                                                                                         
TTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTG  >  AM260479/3247838‑3248080

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