Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,878,931	Δ1 bp	coding (842/918 nt)	cas1 ←	multifunctional nuclease Cas1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,878,928	0	G	.	100.0%	329.8 / NA	78	coding (845/918 nt)	cas1	multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (38/40);  total (38/40)

TGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGA  >  NC_000913/2878793‑2879062                                                                                                                                        |                                                                                                                                       tGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTgggcc                                                                                                                                    >  1:1154713/1‑138 (MQ=255)   gCTAATCTGCCTCGTAAGCGCGTAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTgggccgg                                                                                                                                  <  2:453703/140‑5 (MQ=255)      aaTCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCcggcgg                                                                                                                              <  1:1031285/140‑1 (MQ=255)      aaTCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCcggcgg                                                                                                                              <  2:555693/140‑1 (MQ=255)            ccTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGTCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTAt                                                                                                                        <  2:152522/140‑1 (MQ=255)            ccTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGTCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTAt                                                                                                                        >  1:152522/1‑140 (MQ=255)            ccTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTAt                                                                                                                        >  1:690953/1‑140 (MQ=255)             cTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTAtt                                                                                                                       >  1:499392/1‑140 (MQ=255)             cTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTAtt                                                                                                                       >  1:1147262/1‑140 (MQ=255)             cTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTAtt                                                                                                                       >  2:633944/1‑140 (MQ=255)                 tAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTtctc                                                                                                                   >  1:1201660/1‑140 (MQ=255)                 tAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTtctc                                                                                                                   >  2:717679/1‑140 (MQ=255)                          aGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCg                                                                                                              <  1:532289/136‑1 (MQ=255)                          aGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCg                                                                                                              >  2:532289/1‑136 (MQ=255)                             tACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCa                                                                                                       <  2:70320/140‑1 (MQ=255)                             tACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCa                                                                                                       <  1:789043/140‑1 (MQ=255)                             tACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCa                                                                                                       <  1:898510/140‑1 (MQ=255)                             tACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCa                                                                                                       <  1:973112/140‑1 (MQ=255)                             tACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCa                                                                                                       <  2:1201660/140‑1 (MQ=255)                             tACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCa                                                                                                       <  2:1147262/140‑1 (MQ=255)                               cATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACg                                                                                                     >  2:340532/1‑140 (MQ=255)                               cATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACg                                                                                                     <  2:238195/140‑1 (MQ=255)                               cATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACg                                                                                                     <  2:964159/140‑1 (MQ=255)                                aTTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGt                                                                                                    <  1:185632/140‑1 (MQ=255)                                 ttttCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTc                                                                                                   >  2:1145250/1‑140 (MQ=255)                                 ttttCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTc                                                                                                   <  2:779778/140‑1 (MQ=255)                                    tCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCctc                                                                                                >  1:844324/1‑140 (MQ=255)                                      aGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTc                                                                                                   <  2:99417/135‑1 (MQ=255)                                      aGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTc                                                                                                   >  1:99417/1‑135 (MQ=255)                                      aGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCta                                                                                              <  2:499392/140‑1 (MQ=255)                                         gACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATaa                                                                                           <  1:1145250/140‑1 (MQ=255)                                          aCCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCAGCGGGCAGCACGACCTCTATaaa                                                                                          >  2:293646/1‑139 (MQ=255)                                          aCCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCCGCGGCAAGAACGTCCTCTCTaac                                                                                          >  2:423012/1‑139 (MQ=255)                                          aCCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAg                                                                                          <  2:792897/140‑1 (MQ=255)                                            cACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTc                                                                                                   >  1:226809/1‑129 (MQ=255)                                            cACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTc                                                                                                   <  2:226809/129‑1 (MQ=255)                                                aCCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCCGCGGCAAGCACGTCCTCTATAAACgggaa                                                                                    >  2:480607/1‑137 (MQ=255)                                                     cATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTGTTTTTCCAGCGGCAAACACGGCCTCCATAAGCGGGAACAAtt                                                                               >  2:1210443/1‑140 (MQ=255)                                                      aTACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCCGCGGCAAGCACGGCCCCTATAAACGGGATCAAttt                                                                              >  2:589884/1‑140 (MQ=255)                                                         cTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGTGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGc                                                                           <  2:913452/140‑1 (MQ=255)                                                          tCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCt                                                                          <  1:668937/140‑1 (MQ=255)                                                          tCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCt                                                                          <  1:633944/140‑1 (MQ=255)                                                           cATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTa                                                                         <  2:690953/140‑1 (MQ=255)                                                               tCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTa                                                                         <  2:1067792/136‑1 (MQ=255)                                                               tCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTa                                                                         >  1:1067792/1‑136 (MQ=255)                                                                 aGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCCGCGGCAAGCACGCCCCCTATAAACGGGATCAATTTTGCTAATGttt                                                                   >  2:914321/1‑140 (MQ=255)                                                                   cTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCt                                                                          >  2:981570/1‑131 (MQ=255)                                                                   cTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCt                                                                          <  1:981570/131‑1 (MQ=255)                                                                   cTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTta                                                                 <  2:1154713/140‑1 (MQ=255)                                                                   cTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTta                                                                 >  1:234606/1‑140 (MQ=255)                                                                      cTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCCGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTtacta                                                              >  2:1108938/1‑140 (MQ=255)                                                                      cTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTtacta                                                              <  1:717679/140‑1 (MQ=255)                                                                      cTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGGCCTCTATAAGCGGAATCAATTTTGCTAATGTTTttcta                                                              >  1:152595/1‑140 (MQ=255)                                                                         cGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTg                                                           <  1:589884/140‑1 (MQ=255)                                                                           aTGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTta                                                                 >  2:335636/1‑132 (MQ=255)                                                                           aTGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTta                                                                 <  1:335636/132‑1 (MQ=255)                                                                                cTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGGATTTCTCCCGCGGGGAGCACGTCCCCTATAAACGGGAACAATTTTGCTCATGTTTTTCTACTTCGaaaaa                                                    >  2:555895/1‑140 (MQ=255)                                                                                 tGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGGATCAATTTGGCTAATGTTTTACTACTGCGaaaaaa                                                   >  2:1038747/1‑139 (MQ=255)                                                                                   cATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCCGCGGCCAGCACGTCCTCTATAAGCGGGATCAATTTTGCTAATGTTTTTCTACTGCGgaaaaaaa                                                 >  2:1133669/1‑137 (MQ=255)                                                                                   cATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAatat                                                 <  1:555895/140‑1 (MQ=255)                                                                                   cATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAatat                                                 <  1:480607/140‑1 (MQ=255)                                                                                        cccAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCAGCGGCCAGCACGTCCTCCTTAAACGGAATCAATTTTGGTAATGTTTTTCTACTGCGGAAAAAAAccccc                                            >  2:418864/1‑138 (MQ=255)                                                                                        cccAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTg                                            <  1:1210443/140‑1 (MQ=255)                                                                                          cAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCa                                          <  2:107328/140‑1 (MQ=255)                                                                                                aaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCAc                                         <  1:1132291/135‑1 (MQ=255)                                                                                                aaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCAc                                         >  2:1132291/1‑135 (MQ=255)                                                                                                  cAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAAc                                  <  1:591319/140‑1 (MQ=255)                                                                                                   aGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTTGCTAATGTTTTACTACTGCGAAAAAAATCCCCGCACGCCCAACg                                 >  2:548809/1‑140 (MQ=255)                                                                                                          ggAATGGCAACAGGCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTcc                          <  1:1189113/140‑1 (MQ=255)                                                                                                                      ggCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCCAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCa              >  1:40673/1‑140 (MQ=255)                                                                                                                      ggCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCgg                  >  1:201857/1‑136 (MQ=255)                                                                                                                      ggCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCgg                  <  2:201857/136‑1 (MQ=255)                                                                                                                      ggCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAACGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCa              >  1:2401/1‑140 (MQ=255)                                                                                                                       gCTGTGCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCAc             <  2:844324/140‑1 (MQ=255)                                                                                                                         tgtgCATCTTCAGGT‑GGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTTCTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCa           >  2:47083/1‑140 (MQ=255)                                                                                                                          gtgCATCTTCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAg          >  1:1043512/1‑140 (MQ=255)                                                                                                                                 ttCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTacg   >  1:194083/1‑140 (MQ=255)                                                                                                                                  tCAGGT‑GGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGGATCAATTTGGCTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTacga  >  2:1028008/1‑140 (MQ=255)                                                                                                                                  tCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTacga  <  1:1133669/140‑1 (MQ=255)                                                                                                                                  tCAGGT‑GGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTacga  <  1:914321/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       TGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGA  >  NC_000913/2878793‑2879062

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