Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	810,116	A→G	Y258C (TAC→TGC)	bioB →	biotin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	810,116	0	A	G	100.0%	215.5 / NA	57	Y258C (TAC→TGC)	bioB	biotin synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/33);  total (24/33)

GCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAAC  >  NC_000913/809971‑810253                                                                                                                                                  |                                                                                                                                          gCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTgc                                                                                                                                      <  1:720864/151‑1 (MQ=255)  cAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTtgcgtgcg                                                                                                                                     <  2:992240/151‑1 (MQ=255)  cAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTtgcgtgcg                                                                                                                                     >  1:968242/1‑151 (MQ=255)   aaaCCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTgcgc                                                                                                                                    >  2:255238/1‑151 (MQ=255)         gCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCt                                                                                                                              <  2:279162/151‑1 (MQ=255)         gCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCt                                                                                                                              <  2:125245/151‑1 (MQ=255)         gCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCt                                                                                                                              <  1:546497/151‑1 (MQ=255)            gACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGcc                                                                                                                           >  1:417508/1‑151 (MQ=255)             aCGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCg                                                                                                                          <  2:121355/151‑1 (MQ=255)              cgccgcCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCgg                                                                                                                         <  2:491057/151‑1 (MQ=255)                ccgccgGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGAc                                                                                                                       <  2:332135/151‑1 (MQ=255)                  gccgGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGAcgc                                                                                                                     >  1:1132455/1‑151 (MQ=255)                   ccgGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGAcgc                                                                                                                     >  2:866255/1‑150 (MQ=255)                    cgGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGAcgcg                                                                                                                    >  2:75392/1‑150 (MQ=255)                     ggAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGAcgcg                                                                                                                    >  2:1175680/1‑149 (MQ=255)                     ggAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAg                                                                                                                  <  2:511921/151‑1 (MQ=255)                               ccAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATgaacg                                                                                                        >  2:812370/1‑151 (MQ=255)                                    cAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATgaacgaac                                                                                                     >  2:1072428/1‑149 (MQ=255)                                    cAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACa                                                                                                    <  2:416059/150‑1 (MQ=255)                                    cAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAg                                                                                                   <  1:210897/151‑1 (MQ=255)                                      aCATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGa                                                                                                  <  1:255182/150‑1 (MQ=255)                                       cATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACt                                                                                                >  2:885490/1‑151 (MQ=255)                                             ggtgaaggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCg                                                                                          <  1:210640/151‑1 (MQ=255)                                                gaaggtgaaAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAtg                                                                                       <  1:80717/151‑1 (MQ=255)                                                       aaaGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTa                                                                                <  1:452289/151‑1 (MQ=255)                                                        aaGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTAt                                                                               <  2:417609/151‑1 (MQ=255)                                                         aGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATg                                                                              <  1:337091/151‑1 (MQ=255)                                                           gCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGc                                                                            >  1:566920/1‑151 (MQ=255)                                                             aCGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAg                                                                          <  2:439507/151‑1 (MQ=255)                                                                     tGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCaaa                                                                   <  1:918052/150‑1 (MQ=255)                                                                       ccGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTc                                                                <  2:820404/151‑1 (MQ=255)                                                                       ccGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTc                                                                <  1:1047845/151‑1 (MQ=255)                                                                             aCGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGAtttt                                                          >  1:427281/1‑151 (MQ=255)                                                                                  gatgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACg                                                     <  1:108091/151‑1 (MQ=255)                                                                                    tgTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGt                                                   >  2:249984/1‑151 (MQ=255)                                                                                     gTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGtt                                                  <  2:164395/151‑1 (MQ=255)                                                                                        gATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGc                                                >  1:808565/1‑150 (MQ=255)                                                                                         aTGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGc                                                >  2:760510/1‑149 (MQ=255)                                                                                         aTGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCaa                                              <  1:60182/151‑1 (MQ=255)                                                                                         aTGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                                               >  2:973714/1‑150 (MQ=255)                                                                                            ccTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAAct                                           >  2:1163277/1‑151 (MQ=255)                                                                                              tttGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAActg                                          >  2:530846/1‑150 (MQ=255)                                                                                                tGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAActgctg                                       >  2:702709/1‑151 (MQ=255)                                                                                                 gATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGa                                      <  2:894987/151‑1 (MQ=255)                                                                                                       aTTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGc                                >  2:833011/1‑151 (MQ=255)                                                                                                           gCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGa                             >  1:684429/1‑150 (MQ=255)                                                                                                            cACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGaa                            >  2:342183/1‑150 (MQ=255)                                                                                                             aCCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATc                          >  1:914706/1‑151 (MQ=255)                                                                                                              ccATTGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATc                          <  1:267149/150‑1 (MQ=255)                                                                                                                  tGCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGga                      <  2:78767/150‑1 (MQ=255)                                                                                                                   gCGGTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGgaa                     <  1:552703/150‑1 (MQ=255)                                                                                                                     ggTCGCGCGGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGgaagaa                  <  1:921753/151‑1 (MQ=255)                                                                                                                        cgcgcgGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGAt               <  1:433327/151‑1 (MQ=255)                                                                                                                          cgcgGATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATa              <  1:1105717/150‑1 (MQ=255)                                                                                                                             ggATCATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATaa             >  1:925254/1‑148 (MQ=255)                                                                                                                                 cATGATGCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTg      <  2:1137609/151‑1 (MQ=255)                                                                                                                                    gatgCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCaa   <  1:812189/151‑1 (MQ=255)                                                                                                                                     atgCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAAc  >  1:1163707/1‑151 (MQ=255)                                                                                                                                     atgCCAACCTCTTGCGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAAc  <  1:839080/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                          GCAAACCTGCCGACGCCGCCGGAAAGCGTGCCAATCAACATGCTGGTGAAGGTGAAAGGCACGCCGCTTGCCGATAACGATGATGTCGATGCCTTTGATTTTATTCGCACCATTGCGGTCGCGCGGATCATGATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGGAAGAAGATAAAGACCTGCAAC  >  NC_000913/809971‑810253

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