Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,629,247	C→T	I11I (ATC→ATT)	ydfZ →	putative selenoprotein YdfZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,629,247	0	C	T	92.1%	98.7 / ‑0.1	38	I11I (ATC→ATT)	ydfZ	putative selenoprotein YdfZ
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (18/17);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTT  >  NC_000913/1629003‑1629485                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               cGCCAGGAAAACAGCGACGGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCGTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACt                                                                                                                                                                                                                                          <  2:393059/251‑1 (MQ=255)          aaCAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATtacc                                                                                                                                                                                                                                             >  2:162069/1‑239 (MQ=255)          aaCAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATtacc                                                                                                                                                                                                                                             <  1:162069/239‑1 (MQ=255)               cGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTAt                                                                                                                                                                                                                            >  1:13759/1‑251 (MQ=255)               cGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTAt                                                                                                                                                                                                                            >  2:146194/1‑251 (MQ=255)                gACTGGTTTAGCGAAGAGTACTTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATg                                                                                                                                                                                                                           >  2:184230/1‑251 (MQ=255)                            gAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCAcc                                                                                                                                                                                                               <  1:398398/251‑1 (MQ=255)                            gAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCAcc                                                                                                                                                                                                               <  1:146194/251‑1 (MQ=255)                                    aTTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATTGCTCATGCATGCATCAAAAAAGATGTGAGCTGGATCAAAAACAAAAAA‑GATTTCACTCGACAGGATTATTTATATTGCGCCCGTTACTTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACGTGCAGCCGTGTTATGGTTTGCGGCACCGGTCacac                                                                                                                                                                                                       <  2:189828/251‑1 (MQ=255)                                             ctctCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCa                                                                                                                                                                                                  >  1:404282/1‑247 (MQ=255)                                                                                           catgcatgCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAACTCcgacgc                                                                                                                                                >  2:236/1‑251 (MQ=255)                                                                                           catgcatgCAACAAAAAATGTGTGAGGTTGA‑CCAAAACAAAAAAGAATTGCACGCGACTGGAGTATTTTTATTGCGCCCTTTACGGGGGCTTCGACTGTAAATCAGAAAGGAGAAAACCCCTATGACGACCTACGAACGTAACCGTAACGCAATCGCCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGAACCAGTCGATTGATACTGAAGGTCGAACCGCGGAGGAAATccgtc                                                                                                                                                  <  2:164689/249‑3 (MQ=255)                                                                                                       aaaaaaGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATcc                                                                                                                                                     >  2:312298/1‑234 (MQ=255)                                                                                                       aaaaaaGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATcc                                                                                                                                                     <  1:312298/234‑1 (MQ=255)                                                                                                                gtgAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGaa                                                                                                                           <  1:237108/251‑1 (MQ=255)                                                                                                                   aGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGt                                                                                                                        <  1:437030/251‑1 (MQ=255)                                                                                                                             aaaacaaaaaA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGc                                                                                                                                                                                                                                       <  1:110762/130‑1 (MQ=255)                                                                                                                             aaaacaaaaaA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGc                                                                                                                                                                                                                                       >  2:110762/1‑130 (MQ=255)                                                                                                                                              cACTCGACTGGAGTATTTATATGGCGCCCGTTACGAGGGCTTCGACTGGAAAACAAAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTgg                                                                                               <  2:119158/250‑1 (MQ=255)                                                                                                                                                            aTTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCg                                                                                >  2:262208/1‑251 (MQ=255)                                                                                                                                                                       gcCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAACCACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAgtt                                                                                                                                 >  1:161104/1‑191 (MQ=255)                                                                                                                                                                       gcCCGTTACGGGGGCTTCGACTGTAAATCAGAAAGGAGAAAACGCCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACGGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAgtt                                                                                                                                 <  2:161104/191‑1 (MQ=255)                                                                                                                                                                         ccGTTACGTGGGCTTCGACTGTAAATCAGAAAAGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGGGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGa                                                                                                                   <  2:233019/203‑1 (MQ=255)                                                                                                                                                                         ccGTTACGTGGGCTTCGACTGTAAATCAGAAAAGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGGGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGa                                                                                                                   >  1:233019/1‑203 (MQ=255)                                                                                                                                                                             tACGTGGGCTTCTACTGAAAAACAGAAAGGAGAAAACACCTATGACGCCCTACGATCGTAACCGTAACGCAATTACCACTGCCAGCCGTGTTATGGTTAGCGGCTCCGGTCACACTGGCAAGATCCTGTCGATGGATACTGAAGGTCGGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgtg                                                                <  2:25872/250‑1 (MQ=255)                                                                                                                                                                              acgggggCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACTACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGGGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgtgt                                                               <  2:453463/246‑1 (MQ=255)                                                                                                                                                                                gTGGGCTTCGACTGTAAATCAGAAAGGAGACAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGCTTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGa                                                             >  2:40373/1‑250 (MQ=255)                                                                                                                                                                                      ttCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATgccgc                                                      >  1:472357/1‑251 (MQ=255)                                                                                                                                                                                          aCTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAAcc                                                                                                                                       <  1:136753/166‑1 (MQ=255)                                                                                                                                                                                                       aaGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGAACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtt                                     >  2:217948/1‑251 (MQ=255)                                                                                                                                                                                                       aaGGAGAAAACACCTATGACGACCTACGACCGTAACCGTAACGCAATTACCCCTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtt                                     <  2:13759/251‑1 (MQ=255)                                                                                                                                                                                                            gaAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAgttg                                                                                                                                <  2:286913/155‑1 (MQ=255)                                                                                                                                                                                                            gaAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAgttg                                                                                                                                >  1:286913/1‑155 (MQ=255)                                                                                                                                                                                                            aaaaaCACCTATGACGACCTACGATCTTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                <  1:246748/250‑1 (MQ=255)                                                                                                                                                                                                              aaaCACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACt                              >  1:348813/1‑251 (MQ=255)                                                                                                                                                                                                                acacCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGACCTAAGCGTGTGAATGCCGCCGATGGCGGCATTTCTTTTTTACTTc                            >  2:285775/1‑251 (MQ=255)                                                                                                                                                                                                                    cTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTTATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGa                                                                                                                            <  1:347426/151‑1 (MQ=255)                                                                                                                                                                                                                    cTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTTATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGa                                                                                                                            >  2:347426/1‑151 (MQ=255)                                                                                                                                                                                                                             acTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGATAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCAGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACAGGCGACAATGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAttgtg               >  2:180515/2‑248 (MQ=255)                                                                                                                                                                                                                                          ccgtaacGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGTTTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGACACGGTCGCttt  >  2:26401/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               CGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAA‑TATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTT  >  NC_000913/1629003‑1629485

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