Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	398,073	G→A	G68R (GGG→AGG)	yaiW →	surface‑exposed outer membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	398,073	0	G	A	100.0%	107.3 / NA	34	G68R (GGG→AGG)	yaiW	surface‑exposed outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/18);  total (16/18)

AAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGTCGTGCTGAACGGATGCACATTCCTGCTTTTCTGGTTCATACCGCGCTAAAAATCAAATCGCCAAACGGCAAAAGTTATAGCGAACGGTTAGATTCGGTGCGTACA  >  NC_000913/397828‑398300                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    aaaTTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCGTTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAATGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTg                                                                                                                                                                                                                                >  2:209447/1‑251 (MQ=255)             gaaCCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCGCGCCACCACCTTTGAAAGTCAGAGGCTGGCTCCGACGCTg                                                                                                                                                                                                                    <  2:259149/250‑1 (MQ=255)                 ccATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCCACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGACAGTCAGAGGCTGGCACCGACGCTGGaa                                                                                                                                                                                                                 >  2:93894/1‑250 (MQ=255)                          gAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGGTTTAGCGGGGGGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAAATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGGCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCt                                                                                                                                                                                                       >  2:147557/1‑251 (MQ=255)                              tAAAAAGGAGACGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGAGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGCGCGGGAAAGACCGCGATGCCTGGGCGGAAGATCTCGCCACCAACTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGGTCggg                                                                                                                                                                                                   >  2:17712/1‑250 (MQ=255)                                   aGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTgg                                                                                                                                                                                              <  2:293484/251‑1 (MQ=255)                                        aaGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTg                                                                                                                                                                                         <  1:290185/251‑1 (MQ=255)                                                 cgcgcgTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCCGCTCGGGGCcggcggggggcgggc                                                                                                                                                                                  >  2:278538/1‑240 (MQ=255)                                                  gcgcgTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAgg                                                                                                                                                                               <  1:14431/251‑1 (MQ=255)                                                        cAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGGGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAgggggcc                                                                                                                                                                          >  1:302975/1‑245 (MQ=255)                                                          aTCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGGCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCgggggggcacacgggggcaaa                                                                                                                                                                        >  2:147412/1‑231 (MQ=255)                                                            ccTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTGTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTa                                                                                                                                                                     <  2:302975/251‑1 (MQ=255)                                                              tctATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGGGGGgcagccgggggccaaatatt                                                                                                                                                                   >  2:265153/1‑236 (MQ=255)                                                                 aTCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATc                                                                                                                                                                   <  1:131712/248‑1 (MQ=255)                                                                           tCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGGGGCGCAGCAGGAGGCAAATTATCAGGCCGATCCGGc                                                                                                                                                      >  1:113238/1‑251 (MQ=255)                                                                            cATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGGAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCggt                                                                                                                                                      <  2:319975/250‑2 (MQ=255)                                                                                  tgGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTcc                                                                                                                                                <  1:209447/250‑1 (MQ=255)                                                                                            ggTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGCGCGGCAGAAGATGCCTGCGATCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCGGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGGTCCGTCCTTTCCGGGCTTAAg                                                                                                                                       <  2:193580/249‑1 (MQ=255)                                                                                                            agtagtCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCGGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAag                                                                                                                      <  2:141644/250‑1 (MQ=255)                                                                                                                   aGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGAc                                                                                                              <  1:78440/251‑1 (MQ=255)                                                                                                                        ccGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCa                                                                                                                        <  1:27856/236‑1 (MQ=255)                                                                                                                        ccGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAagatcggaagagggg                                                                                                         >  2:27856/1‑239 (MQ=255)                                                                                                                                             gAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCCAATCCGGCCGTTCCGGGCTTAAGCAAAAAAGCCTGGCAAGCGATTGACCGTCGTGCTGAACGGATGCCCATTcc                                                                                    >  2:195345/1‑251 (MQ=255)                                                                                                                                              aaaaaGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAACGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGGGCAGCAGGAGTCAAATTATCAGGCAGATACGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGTCGTGCTGGACGGATGCACATTCCt                                                                                   >  1:57806/1‑251 (MQ=255)                                                                                                                                                    cGATAGATGTAGCGGGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGTCGTGCTGAACGGATGCACATTCCTGCtttt                                                                             <  2:135439/251‑1 (MQ=255)                                                                                                                                                     gATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTTAGCAAGATCGCCTGGCAAGAGATTGACCCTCGTGCTGAACGGATGCACATTCCTGCttttt                                                                            >  1:202743/1‑250 (MQ=255)                                                                                                                                                            gTAGCGAGCGTCGGGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTGCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGCCCGTCGTGCTGAACGGATGCTCATTCCTGCTTTTCTGGTTc                                                                      <  2:405330/250‑1 (MQ=255)                                                                                                                                                              agcgagcgTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGGCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGGCGTGCTGAACGGATGCACAATCCTGCTTTTCTGGTTcctt                                                                   >  1:397205/1‑248 (MQ=255)                                                                                                                                                                                gaagaTGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGTCGTGCTGAACGGATGCACATTCCTGCTTTTCTGGTTCATACCGCGCTAAaaatcaaat                                                 <  1:209339/251‑1 (MQ=255)                                                                                                                                                                                          gCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTTTCCGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGTCGTGCTGAACGGCTGCACATTCCTGCTTTTCTGGTTCATACCGCGCTAAAAATCAAATCGCCAAACgg                                       >  2:277694/1‑251 (MQ=255)                                                                                                                                                                                                gTGAAAGACCGCGATGCCTGGGCGAAATATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCGGGCGGTGGCGCAGCAGGAGGCAAATTATCAGGCAGATCCGGCCGTTCCGGGGTTAAGCAAGATCTCCTGGCAAGAGAGTGACCGTCGTGCTGAACGGATGCACATTCCTGCTGGTCTGGTGCATACCGCGCTAAAAATCAAATCGCCAAACGGCAAAAg                                 <  1:269829/251‑1 (MQ=255)                                                                                                                                                                                                                ccTGGGCGAACGATCTCGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGCGGCGCAGCAGGAGGCAAATTATCAGGCAGATCCGGCCGGTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGCTTGACCGTCGTGCTGAACGGATTCACATTCCTGCTTTTCTGGTTCATACCGCGCTAAAAATCCAATCGCCAAACGGCAAAAGTTATAGCGAACGGttt                 >  2:48561/1‑250 (MQ=255)                                                                                                                                                                                                                                cGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGTCGTGCTGAACGGATGCACATTCCTGCTTTTCTGGTTCATACCGCGCTAAAAATCAAATCGCCAAACGGCAAAAGTTATAGCGAACGGTTAGATTCGGTGCGTa    <  1:345049/248‑1 (MQ=255)                                                                                                                                                                                                                                cGCCACCACCTTTGAAAGTCAGAGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGTCGTGCTGAACGGATGCACATTCCTGCTTTTCTGGTTCATACCGCGCTAAAAATCAAATCGCCAAACGGCAAAAGTTATAGCGAACGGTTAGATTCGGTGCGTACa  <  1:265153/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    AAATTCAGGAAGTAACCCATACCTTGAGCTAAAAAGGAGAAGTAATGTCGCGCGTCAATCCTCTATCTTCGCTGTCATTGTTGGCGGTTCTGGTTTTAGCGGGGTGTAGTAGTCAGGCTCCGCAACCGTTGAAAAAGGGCGAAAAAGCGATAGATGTAGCGAGCGTCGTGCGGCAGAAGATGCCTGCGAGCGTGAAAGACCGCGATGCCTGGGCGAAAGATCTCGCCACCACCTTTGAAAGTCAGGGGCTGGCACCGACGCTGGAAAACGTCTGCTCGGTGCTGGCGGTGGCGCAGCAGGAGTCAAATTATCAGGCAGATCCGGCCGTTCCGGGCTTAAGCAAGATCGCCTGGCAAGAGATTGACCGTCGTGCTGAACGGATGCACATTCCTGCTTTTCTGGTTCATACCGCGCTAAAAATCAAATCGCCAAACGGCAAAAGTTATAGCGAACGGTTAGATTCGGTGCGTACA  >  NC_000913/397828‑398300

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