Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,626,944	C→T	L665L (CTG→CTA)	pdeF ←	cyclic di‑GMP phosphodiesterase PdeF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,626,944	0	C	T	100.0%	219.8 / NA	64	L665L (CTG→CTA)	pdeF	cyclic di‑GMP phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/35);  total (29/35)

TCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCA  >  NC_000913/2626816‑2627076                                                                                                                                 |                                                                                                                                     tCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCAtt                                                                                                                           <  2:920492/140‑1 (MQ=255)   gTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTca                                                                                                                         <  1:612969/140‑1 (MQ=255)   gTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTAAATTTTAAGTAGATCGGCATTca                                                                                                                         <  1:366936/140‑1 (MQ=255)          aCGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACAttttt                                                                                                                  <  2:732422/140‑1 (MQ=255)             tACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAg                                                                                                               <  2:2472/140‑1 (MQ=255)                   gCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCgcgc                                                                                                         <  2:993122/140‑1 (MQ=255)                   gCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCgcgc                                                                                                         >  1:912366/1‑140 (MQ=255)                   gCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCgcgc                                                                                                         <  1:597032/140‑1 (MQ=255)                    ccACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCa                                                                                                        >  2:32618/1‑140 (MQ=255)                       cTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAg                                                                                                     <  1:201355/140‑1 (MQ=255)                       cTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAg                                                                                                     <  2:717081/140‑1 (MQ=255)                       cTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAg                                                                                                     >  2:440766/1‑140 (MQ=255)                       cTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAg                                                                                                     >  2:315201/1‑140 (MQ=255)                         aCCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCt                                                                                                   >  1:376718/1‑140 (MQ=255)                              cATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCcgcg                                                                                                          >  2:643478/1‑128 (MQ=255)                              cATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCcgcg                                                                                                          <  1:643478/128‑1 (MQ=255)                              cATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGt                                                                                              >  1:813246/1‑140 (MQ=255)                              cATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGt                                                                                              >  2:221854/1‑140 (MQ=255)                                      tcATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTg                                                                                      <  2:681125/140‑1 (MQ=255)                                        aTTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGcc                                                                                    >  2:67145/1‑140 (MQ=255)                                        aTTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCAGGCGTAGCCGGTGcc                                                                                    <  2:632860/140‑1 (MQ=255)                                             cgcCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAat                                                                               >  1:221852/1‑140 (MQ=255)                                                             gATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTg                                                               >  1:470574/1‑140 (MQ=255)                                                              aTGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTgg                                                              <  2:813246/140‑1 (MQ=255)                                                              aTGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTgg                                                              >  2:135877/1‑140 (MQ=255)                                                                    cTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCt                                                        <  2:470574/140‑1 (MQ=255)                                                                    cTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCt                                                        <  1:960962/140‑1 (MQ=255)                                                                      aTCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTaa                                                      <  1:32618/140‑1 (MQ=255)                                                                           aTAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCt                                                 >  2:564791/1‑140 (MQ=255)                                                                             aaTCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGa                                               <  1:135877/140‑1 (MQ=255)                                                                               tCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAg                                             <  1:221854/140‑1 (MQ=255)                                                                               tCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAg                                             >  1:847666/1‑140 (MQ=255)                                                                               tCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAg                                             >  1:354338/1‑140 (MQ=255)                                                                                ccAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGa                                            <  2:376718/140‑1 (MQ=255)                                                                                 cAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCa                                                            <  1:601339/123‑1 (MQ=255)                                                                                 cAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCa                                                            >  2:601339/1‑123 (MQ=255)                                                                                 cAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAt                                           >  1:714352/1‑140 (MQ=255)                                                                                 cAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAt                                           >  1:301900/1‑140 (MQ=255)                                                                                  aGACTATTTGACACAATATTGCGGATAAACCTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAtg                                          <  2:108294/140‑1 (MQ=255)                                                                                   gACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAtgc                                         >  1:810601/1‑140 (MQ=255)                                                                                     cTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGAtgctg                                       >  2:856903/1‑140 (MQ=255)                                                                                        tttGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCaa                                    >  1:78061/1‑140 (MQ=255)                                                                                         ttGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAg                                   >  2:257338/1‑140 (MQ=255)                                                                                         ttGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAg                                   <  2:147538/140‑1 (MQ=255)                                                                                          tGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAgg                                  <  1:315201/140‑1 (MQ=255)                                                                                          tGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAgg                                  <  1:440766/140‑1 (MQ=255)                                                                                             cacaATATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGa                               >  1:623950/1‑140 (MQ=255)                                                                                                aaTATTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGAttt                            >  1:2227/1‑140 (MQ=255)                                                                                                  tatTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgc                          >  1:689109/1‑140 (MQ=255)                                                                                                  tatTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgc                          >  1:281402/1‑140 (MQ=255)                                                                                                   atTGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcg                         <  2:931801/140‑1 (MQ=255)                                                                                                    ttGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                        <  2:301900/140‑1 (MQ=255)                                                                                                    ttGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                        <  2:148818/140‑1 (MQ=255)                                                                                                    ttGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                        <  2:714352/140‑1 (MQ=255)                                                                                                    ttGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                        <  2:912366/140‑1 (MQ=255)                                                                                                    ttGCGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTgcgc                        <  1:317033/140‑1 (MQ=255)                                                                                                       cGGATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTg                     >  1:987667/1‑140 (MQ=255)                                                                                                         gATAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCt                   <  2:810601/140‑1 (MQ=255)                                                                                                          aTAAAACTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCtt                  >  2:519128/1‑140 (MQ=255)                                                                                                            aaaaCTGCCGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTaa                >  2:154409/1‑140 (MQ=255)                                                                                                                   ccGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGt         <  2:623950/140‑1 (MQ=255)                                                                                                                   ccGTCAATTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGt         <  2:987667/140‑1 (MQ=255)                                                                                                                        aaTTTTAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCagag    <  1:539633/140‑1 (MQ=255)                                                                                                                          ttttAAGTAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCa  <  1:519128/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                     TCGTTTTCAACGTACTCTGCCACTACCAGCATTTTCTTCATTCGCGCCAGGTGGCAAATCGATGCCACTATCTGATAATCCAGACTATTTGACACAATATTGCGGATAAAACTGCCGTCAATTTTAAGCAGATCGGCATTCACATTTTTAAGCCGCGCATAGCTGGCGTAGCCGGTGCCGAAATCATCAATCGCAATCTGGCAGCCTAATTCCTGAAGATGCTGCAAGGTGATTTGCGCCTGCTTAACATTGGTCAGAGCA  >  NC_000913/2626816‑2627076

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