Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,928,311	Δ1 bp	coding (529/2061 nt)	ptrB ←	oligopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,928,309	0	A	.	100.0%	229.9 / NA	54	coding (531/2061 nt)	ptrB	oligopeptidase B
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (27/27);  total (27/27)

TATCGTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACAAAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCC  >  NC_000913/1928173‑1928444                                                                                                                                         |                                                                                                                                        tATCGTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACaagc                                                                                                                                      >  2:236023/1‑138 (MQ=255) tATCGTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACaagc                                                                                                                                      >  1:740859/1‑138 (MQ=255)  aTCGTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACaagct                                                                                                                                     >  1:871720/1‑137 (MQ=255)  aTCGTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACaagct                                                                                                                                     <  2:429986/140‑4 (MQ=255)  aTCGTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACaagct                                                                                                                                     <  2:611957/140‑4 (MQ=255)  aTCGTCTTTTTCTTCGTAGATCAGTCTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACaagct                                                                                                                                     >  1:135958/1‑137 (MQ=255)     gTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACaagc                                                                                                                                      >  2:213045/1‑134 (MQ=255)     gTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACaagc                                                                                                                                      <  1:213045/136‑3 (MQ=255)     gTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTggg                                                                                                                                 >  1:915137/1‑140 (MQ=255)     gTCTTTTTCTTCGTAGATCAGTCTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTggg                                                                                                                                 <  2:135958/140‑1 (MQ=255)      tCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGt                                                                                                                                >  2:548269/1‑140 (MQ=255)      tCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGt                                                                                                                                >  1:717437/1‑140 (MQ=255)        tttttCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTc                                                                                                                              >  1:216759/1‑140 (MQ=255)          tttCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCaa                                                                                                                            <  2:717437/140‑1 (MQ=255)          tttCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCaa                                                                                                                            <  1:87544/140‑1 (MQ=255)                    gATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCcag                                                                                                                  <  2:915137/140‑1 (MQ=255)                     aTCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCcagc                                                                                                                 >  1:566192/1‑140 (MQ=255)                          tttATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTc                                                                                                            >  1:54282/1‑140 (MQ=255)                             aTCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCgg                                                                                                         <  2:216759/140‑1 (MQ=255)                               cTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGt                                                                                                       >  2:72953/1‑140 (MQ=255)                               cTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGt                                                                                                       >  1:290742/1‑140 (MQ=255)                                ttGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGta                                                                                                      <  2:167205/140‑1 (MQ=255)                                  gCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGtt                                                                                                             >  1:905593/1‑131 (MQ=255)                                  gCGATGCTCGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGtt                                                                                                             <  2:905593/131‑1 (MQ=255)                                   cGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGGGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGtacca                                                                                                   >  1:534392/1‑140 (MQ=255)                                   cGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGtacca                                                                                                   >  2:835373/1‑140 (MQ=255)                                        cTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGtaccagttac                                                                                              <  1:835373/140‑1 (MQ=255)                                        cTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGtaccagttac                                                                                              <  1:682821/140‑1 (MQ=255)                                        cTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGCTCCGGGtaccagttac                                                                                              <  2:201004/140‑1 (MQ=255)                                          ggCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGtaccagttacca                                                                                            <  2:290742/140‑1 (MQ=255)                                           gCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGtaccagttaccag                                                                                           <  2:323333/140‑1 (MQ=255)                                               aCCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTc                                                                                       <  2:37778/140‑1 (MQ=255)                                                ccGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTATAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGtaccagttac                                                                                              >  2:136891/1‑132 (MQ=255)                                                ccGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTATAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGtaccagttac                                                                                              <  1:136891/132‑1 (MQ=255)                                                ccGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTcc                                                                                      >  1:652616/1‑140 (MQ=255)                                                         tGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGGGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGa                                                                             <  2:534392/140‑1 (MQ=255)                                                                aGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGaa                                                                      <  2:930538/140‑1 (MQ=255)                                                                     tGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCg                                                                 <  1:548269/140‑1 (MQ=255)                                                                      gATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGt                                                                >  2:550426/1‑140 (MQ=255)                                                                            gcagcaGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAAc                                                                         <  2:660056/125‑1 (MQ=255)                                                                            gcagcaGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAAc                                                                         <  2:616505/125‑1 (MQ=255)                                                                            gcagcaGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAAc                                                                         >  1:660056/1‑125 (MQ=255)                                                                            gcagcaGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAAc                                                                         >  1:616505/1‑125 (MQ=255)                                                                                cagcGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGaa                                                                      >  2:598172/1‑124 (MQ=255)                                                                                cagcGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGaa                                                                      <  1:598172/124‑1 (MQ=255)                                                                                cagcGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCg                                                      >  2:195255/1‑140 (MQ=255)                                                                                 agcGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCgg                                                     >  1:1005668/1‑140 (MQ=255)                                                                                       aCCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGgaaaga                                               >  1:246072/1‑140 (MQ=255)                                                                                                   cGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGcc                                   >  1:449149/1‑140 (MQ=255)                                                                                                   cGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGcc                                   <  2:1005668/140‑1 (MQ=255)                                                                                                     aaCATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAg                                 <  1:550426/140‑1 (MQ=255)                                                                                                     aaCATAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAg                                 <  1:72953/140‑1 (MQ=255)                                                                                                        aTAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGc                                    >  2:419419/1‑134 (MQ=255)                                                                                                        aTAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGc                                    <  1:419419/134‑1 (MQ=255)                                                                                                        aTAGTAGAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAgcgc                              >  2:90864/1‑140 (MQ=255)                                                                                                           gtagAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCAt                           >  1:403669/1‑140 (MQ=255)                                                                                                            tagAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGc                                    <  1:781443/130‑1 (MQ=255)                                                                                                            tagAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGc                                    >  2:781443/1‑130 (MQ=255)                                                                                                             agAAAATCCAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATaa                         <  1:90864/140‑1 (MQ=255)                                                                                                                    ccAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAg                                 >  2:351168/1‑125 (MQ=255)                                                                                                                    ccAGGAGTCATTTGCCCAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAg                                 <  1:351168/125‑1 (MQ=255)                                                                                                                                    cAGAC‑AAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGcc  >  2:857465/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        TATCGTCTTTTTCTTCGTAGATCAGTTTATCTTGCGATGCTGGCGTACCGATGGCGTGACGCCAGACCTGATAAGGCAGCAGCGTCACCGGATGCTTGCGAACATAGTAGAAAATCCAGGAGTCATTTGCCCAGACAAAGCTGGGTTCAACGTTATCCAGCAGTTCCGGGTACCAGTTACCAGTTTCCAGATTACGAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCC  >  NC_000913/1928173‑1928444

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