Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260480	507,425	+A	coding (27/564 nt)	h16_B0448 →	fragment of ABC‑type transport system,periplasmic component: APC family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260480	507,423	1	.	A	100.0%	134.5 / NA	35	E9K (GAA→AAA)	h16_B0448	fragment of ABC‑type transport system,periplasmic component: APC family
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (17/18);  total (17/18)

AATACCACGGGCGGCGTTCGCCGCCTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATG‑AACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGCCGTCGGCCGATGTCGGGCGCGACCAATA  >  AM260480/507312‑507550                                                                                                                 |                                                                                                                                aaTACCACGGGCGGCGTTCGCCGCCTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATcg                                                                                                            >  2:153780/1‑134 (MQ=255)    aCCACGGGCGGCGTTCGCCGCCTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATcgccg                                                                                                         >  1:10388/1‑134 (MQ=255)            cggcgTTCGCCGCCTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATcgccgccgc                                                                                                     >  2:291105/1‑130 (MQ=255)            cggcgTTCGCCGCCTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATcgccgccgc                                                                                                     <  1:291105/130‑1 (MQ=255)              gcgTTCGCCGCCTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCgcgctggc                                                                                               >  2:517779/1‑134 (MQ=255)                gTTCGCCGCCTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCgcgctggcgc                                                                                             <  1:319281/134‑1 (MQ=255)                     ccgccTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTc                                                                                        <  1:460318/134‑1 (MQ=255)                               tGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGgcgc                                                                              <  1:153780/134‑1 (MQ=255)                               tGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGgcgc                                                                              <  2:841843/134‑1 (MQ=255)                               tGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGgcgc                                                                              <  1:450535/134‑1 (MQ=255)                                 gggAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTg                                                                            <  2:892078/134‑1 (MQ=255)                                       cGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCCGGCGCGAGCAGCGCCAGCGCTGGCCAGt                                                                      <  2:77345/134‑1 (MQ=255)                                          tCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTg                                                                            >  2:300163/1‑125 (MQ=255)                                          tCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTg                                                                            <  1:300163/125‑1 (MQ=255)                                            aGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGgcgc                                                                              <  1:280994/121‑1 (MQ=255)                                            aGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGgcgc                                                                              >  2:280994/1‑121 (MQ=255)                                                  gcgTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCtgaagtga                                                            <  1:258692/133‑1 (MQ=255)                                                  gcgTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCtgaagtga                                                            <  1:517779/133‑1 (MQ=255)                                                       cGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATcgccgccg                                                                                                      <  2:661662/86‑1 (MQ=255)                                                       cGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATcgccgccg                                                                                                      >  1:661662/1‑86 (MQ=255)                                                       cGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCg                                                      >  1:640277/1‑134 (MQ=255)                                                       cGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCg                                                      >  2:85752/1‑134 (MQ=255)                                                                                gACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGcc                             >  2:1097396/1‑134 (MQ=255)                                                                                gACGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGcc                             >  2:1060086/1‑134 (MQ=255)                                                                                 aCGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTgg                                                                           <  1:335378/87‑1 (MQ=255)                                                                                 aCGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTgg                                                                           >  2:335378/1‑87 (MQ=255)                                                                                 aCGGGCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGCCg                            >  1:659478/1‑134 (MQ=255)                                                                                   gggCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGCCGTc                          <  2:659478/134‑1 (MQ=255)                                                                                     gCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGgcgc                                                                              >  2:576852/1‑80 (MQ=255)                                                                                     gCAATGACAATTCGACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGgcgc                                                                              <  1:576852/80‑1 (MQ=255)                                                                                                  gACGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGCCGTCGGCCGATGTCGGgcg           >  1:1078682/1‑134 (MQ=255)                                                                                                   aCGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGcc                                                                         <  1:89389/71‑1 (MQ=255)                                                                                                   aCGGAGGAGACATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGcc                                                                         >  2:89389/1‑71 (MQ=255)                                                                                                          agaCATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGCCGTCGGCCGATGTCGGGCGCGACCAAt   <  2:1078682/134‑1 (MQ=255)                                                                                                           gaCATGAAACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGCCGTCGGCCGATGTCGGGCGCGACCAATa  >  2:1068200/1‑134 (MQ=255)                                                                                                                 |                                                                                                                                AATACCACGGGCGGCGTTCGCCGCCTGACATGGGGAGCCGTTCAGGCTCGCGTCCGCACGGGCACCACATAGAGTGCCTGACGGGCAATGACAATTCGACGGAGGAGACATG‑AACAGAAATGCGTGGCAATCGCCGCCGCGCTGGCGCTGCTCGCGCCGGCGCTGGCCAGTGCTGAAGTGAAGGTCGGCATTCTTGCCACCCTTTCGGGGCCGTCGGCCGATGTCGGGCGCGACCAATA  >  AM260480/507312‑507550

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