Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260480	2,843,580	+G	intergenic (‑1/+68)	pip ← / ← h16_B2516	prolyl aminopeptidase/transcriptional regulator, LysR‑family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260480	2,843,579	1	.	G	100.0%	170.4 / NA	44	M1M (ATG→CTG) †	pip	prolyl aminopeptidase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (28/16);  total (28/16)

CGTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑GATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACG  >  AM260480/2843447‑2843701                                                                                                                                      |                                                                                                                           cgTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           <  1:823170/134‑1 (MQ=255) cgTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           <  1:1038634/134‑1 (MQ=255) cgTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           <  1:381784/134‑1 (MQ=255) cgTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           <  2:601907/134‑1 (MQ=255)     cAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           <  2:218297/130‑1 (MQ=255)     cAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           >  1:218297/1‑130 (MQ=255)                 cgGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACg                                                                                                              >  2:640679/1‑132 (MQ=255)                   gCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACg                                                                                                              >  2:700118/1‑130 (MQ=255)                       ggCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCg                                                                                                       <  1:353982/133‑1 (MQ=255)                         cgcCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                    <  1:700118/134‑1 (MQ=255)                                cGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAg                                                                                             <  1:640679/134‑1 (MQ=255)                                                 agtagATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCcaca                                                                                              <  1:913067/116‑1 (MQ=255)                                                 agtagATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCcaca                                                                                              >  2:913067/1‑116 (MQ=255)                                                             gCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAg                                                                >  2:605501/1‑134 (MQ=255)                                                               ccGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACg                                                                                                              >  1:455836/1‑86 (MQ=255)                                                               ccGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACg                                                                                                              <  2:455836/86‑1 (MQ=255)                                                               ccGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAgcc                                                              >  2:131821/1‑134 (MQ=255)                                                               ccGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAgcc                                                              >  2:709322/1‑134 (MQ=255)                                                               ccGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAgcc                                                              >  1:578991/1‑134 (MQ=255)                                                                 gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGCGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGga                                                                                                            >  1:240002/1‑86 (MQ=255)                                                                 gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGga                                                                                                            <  2:240002/86‑1 (MQ=255)                                                                 gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                    <  2:369452/94‑1 (MQ=255)                                                                 gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                    >  1:369452/1‑94 (MQ=255)                                                                 gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAgcccg                                                            >  2:787729/1‑134 (MQ=255)                                                                 gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAgcccg                                                            >  1:92557/1‑134 (MQ=255)                                                                 gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           >  2:336527/1‑70 (MQ=255)                                                                 gTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           <  1:336527/70‑1 (MQ=255)                                                                       gACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           <  2:622019/64‑1 (MQ=255)                                                                       gACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑g                                                                                                                           >  1:622019/1‑64 (MQ=255)                                                                          aTCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑gga                                                                                                                         >  2:556121/1‑61 (MQ=255)                                                                          aTCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑gga                                                                                                                         <  1:556121/63‑3 (MQ=255)                                                                             cAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCgg                                                >  2:654857/1‑134 (MQ=255)                                                                                 gTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGAGCAAGGCCGGCATGGATATGGGTCTACGgaga                                                                                                          <  2:17416/72‑1 (MQ=255)                                                                                 gTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGAGCAAGGCCGGCATGGATATGGGTCTACGgaga                                                                                                          >  1:17416/1‑72 (MQ=255)                                                                                    ccGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                  >  2:7958/1‑109 (MQ=255)                                                                                    ccGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                  <  1:7958/109‑1 (MQ=255)                                                                                      ggTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCg                                       >  2:1149238/1‑134 (MQ=255)                                                                                        tttCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                    >  1:521537/1‑71 (MQ=255)                                                                                        tttCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                    <  2:521537/71‑1 (MQ=255)                                                                                                      aTTTCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCt                       >  1:537271/1‑134 (MQ=255)                                                                                                         tCGGGATAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAg                    >  2:102853/1‑134 (MQ=255)                                                                                                              aTAGAGGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAAcc               >  2:1129596/1‑134 (MQ=255)                                                                                                                agagGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCca                                                                                                >  2:185295/1‑51 (MQ=39)                                                                                                                agagGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCca                                                                                                <  1:185295/51‑1 (MQ=39)                                                                                                                 gagGGTGCGCAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACcgcg            <  1:787729/134‑1 (MQ=255)                                                                                                                        cgcAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTc                                  <  2:352036/105‑1 (MQ=255)                                                                                                                        cgcAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTc                                  >  1:352036/1‑105 (MQ=255)                                                                                                                          cAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGAc   >  2:683025/1‑134 (MQ=255)                                                                                                                          cAAGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGAc   >  1:1052037/1‑134 (MQ=255)                                                                                                                           aaGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                    <  2:1041399/36‑1 (MQ=255)                                                                                                                           aaGGCCGGCATGGATATGGGTCTACGGAGATCGGGa                                                                                                    >  1:1041399/1‑36 (MQ=255)                                                                                                                           aaGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                  <  2:1040838/70‑1 (MQ=255)                                                                                                                           aaGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGt                                                                  >  1:1040838/1‑70 (MQ=255)                                                                                                                           aaGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACg  >  1:230123/1‑134 (MQ=255)                                                                                                                            aGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCa                             >  2:2688/1‑106 (MQ=255)                                                                                                                            aGGCCGGCATGGATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCa                             <  1:2688/106‑1 (MQ=255)                                                                                                                                      |                                                                                                                           CGTGCAGGAACACCGCCGGCTTGGCGCCGGGCGTGCCCACACGCTCGTAGTAGATGACATGCCCGTCGCCGACATCCAGCGTGCCGGTTTCATAGGGCTCGATTTCGGGATAGAGGGTGCGCAAGGCCGGCAT‑GATATGGGTCTACGGAGATCGGGAGCCACAGTGTACGCCGGAGGGAACGTCAAACAGTAGCCCGCCCCTTAACCGGCGAACACCGCCTTCACCCAGTCCCTCAGCAACCGCGTGGCGGGACG  >  AM260480/2843447‑2843701

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