Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260480	153,689	+G	intergenic (‑67/+93)	h16_B0133 ← / ← h16_B0134	conserved hypothetical protein/outer membrane protein (porin)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260480	153,687	1	.	G	100.0%	139.9 / NA	36	intergenic (‑65/+95)	h16_B0133/h16_B0134	conserved hypothetical protein/outer membrane protein (porin)
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (20/16);  total (20/16)

GTGGTGTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCA‑GGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGTGACGCAGGCCGGCCATGACCGTG  >  AM260480/153564‑153815                                                                                                                             |                                                                                                                                 gtggtgTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCat                                                                                                                          >  2:841960/1‑133 (MQ=255) gtggtgTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCat                                                                                                                          >  2:765297/1‑133 (MQ=255)  tggtgTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCatga                                                                                                                        <  1:521088/134‑1 (MQ=255)    gtgtTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCatgaat                                                                                                                      >  1:996277/1‑134 (MQ=255)      gtTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCatgaatga                                                                                                                    >  2:135626/1‑134 (MQ=255)          gCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAg                                                                                                                <  1:841960/134‑1 (MQ=255)                            gcgTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTctg                                                                                                          >  2:552231/1‑122 (MQ=255)                            gcgTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTctg                                                                                                          <  1:552231/122‑1 (MQ=255)                                     aTCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTtctc                                                                                     >  2:12863/1‑134 (MQ=255)                                        gTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCaa                                                                                  <  2:996277/134‑1 (MQ=255)                                                ttCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTctgc                                                                                                         <  1:751202/103‑1 (MQ=255)                                                ttCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTctgc                                                                                                         >  2:751202/1‑103 (MQ=255)                                                ttCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGAtgtg                                                                          >  2:226675/1‑134 (MQ=255)                                                  ccACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAg                                                                                                                <  2:419903/94‑1 (MQ=255)                                                  ccACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAg                                                                                                                >  1:419903/1‑94 (MQ=255)                                                      gCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGc                                                                    <  1:12863/134‑1 (MQ=255)                                                            gccggcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTc                                                              <  1:765297/134‑1 (MQ=255)                                                              cggcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAg                                                                                                                <  1:81067/82‑1 (MQ=255)                                                              cggcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAg                                                                                                                >  2:81067/1‑82 (MQ=255)                                                              cggcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCAt                                                                                              >  1:109730/1‑100 (MQ=255)                                                              cggcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCAt                                                                                              <  2:109730/100‑1 (MQ=255)                                                              cggcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTtc                                                                                       <  2:433885/107‑1 (MQ=255)                                                              cggcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTtc                                                                                       >  1:433885/1‑107 (MQ=255)                                                                gcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGAt                                                                             <  2:196309/115‑1 (MQ=255)                                                                gcgccggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGAt                                                                             >  1:196309/1‑115 (MQ=255)                                                                     ggcATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCggg                                                      >  1:886734/1‑133 (MQ=255)                                                                               gcgAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCggg                                                      >  2:512621/1‑123 (MQ=255)                                                                               gcgAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCggg                                                      <  1:512621/123‑1 (MQ=255)                                                                                    cAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCaa                                                                                  >  1:380644/1‑90 (MQ=255)                                                                                    cAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCaa                                                                                  <  2:380644/90‑1 (MQ=255)                                                                                           cTGACGACCGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAg                                >  2:664657/1‑133 (MQ=255)                                                                                                   cGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAgg                                                                                                   <  2:625459/58‑1 (MQ=255)                                                                                                   cGCACCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAgg                                                                                                   >  1:625459/1‑58 (MQ=255)                                                                                                      aCCGTTGGCGTGCGGGCTGTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGTGACGc                    <  1:135626/134‑1 (MQ=255)                                                                                                                        gTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGTGACGCAGGCCGGCCATGACCGTg  >  1:433373/1‑134 (MQ=255)                                                                                                                        gTCCAGGGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGTGACGCAGGCCGGCCATGACCGTg  >  1:91168/1‑134 (MQ=255)                                                                                                                             |                                                                                                                                 GTGGTGTTGGCCGTGTTCTCGGGCGGGGCGTTGTGGATCGTGGTCACTTCCACGCTCACGCCGGCGCCGGCATAGGGCGCGAACAGTTCGCTGACGACCGCACCGTTGGCGTGCGGGCTGTCCA‑GGCTTCATGAATGACCAGCGTCTGCTGCAGGGCCATGGCTTTCTCCAAAGGATGTGCCATGCCCGCTCGGCGCGGGCATGGCGTCTACGAAAAGTCAGAACTGGTGACGCAGGCCGGCCATGACCGTG  >  AM260480/153564‑153815

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: ‑