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	97.6%	126.7 / NA	41	I11I (ATC→ATT)	ydfZ	putative selenoprotein YdfZ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/22);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCG  >  NC_000913/1629003‑1629489                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   cGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACt                                                                                                                                                                                                                                              <  2:270194/251‑1 (MQ=255)  gCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCAGTGACTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCAACAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACGCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAAGCAGAAAGGAGAAAACACCTTTGTCGACCTACGATCGTAACCGTAACGCAAGTACCACTg                                                                                                                                                                                                                                             <  2:801973/251‑1 (MQ=255)          aaCAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTCCGATCGTAACCGTAACGCAATTACCACTGGCAGCCgt                                                                                                                                                                                                                                     >  2:201834/1‑251 (MQ=255)             aGCGACTGGTTTAGCGAAGATCAATTCATTTCCTCTCATCCCATCCTGGGTTAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATAAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATGTCACTCGACAGGAGTATTTATAGTGCGCCCGTTACGTGGGCTTCGACTGTAAAGCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGtt                                                                                                                                                                                                                                  <  1:614743/251‑1 (MQ=255)                           cGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCAc                                                                                                                                                                                                                    <  1:89753/251‑1 (MQ=255)                                      tCATTTCCTCTCATCCCATCCGAGGTGAGAGTCTTTTCCCCCACCTAATGGCTCATGCGTGCATCAAAAAAGATGTCACCTTGCTCAAAAACAAACACTATTTCCCTCGCCAGGAGTATTTAAGTTGCGCCCGTTACGTGGGCTTAGTCTGTAAATCAGAAAGGAGAAAACACCTATGACAACCTACGCACGTAACCGTAACGCAATTACCCATGAGAGCCGTGTTATGGttgactgaaccgtccaaaccg                                                                                                                                                                                                         >  2:500138/1‑232 (MQ=255)                                       cATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTgg                                                                                                                                                                                                        <  2:544929/251‑1 (MQ=255)                                             ctctCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACACAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCa                                                                                                                                                                                                      >  2:704058/1‑247 (MQ=255)                                                ccatcccatccGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCt                                                                                                                                                                                               <  1:420101/250‑1 (MQ=255)                                                  atcccatccGGGGTGAGAGTCTTTTCCTCAGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTCACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGt                                                                                                                                                                                             >  1:172638/1‑251 (MQ=255)                                                                    gTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAAGTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTc                                                                                                                                                                           <  1:71744/251‑1 (MQ=255)                                                                                                atgcatCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGaa                                                                                                                                                >  1:94143/1‑250 (MQ=255)                                                                                                       aaaaaaGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGGCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTa                                                                                                                                        <  1:562490/251‑1 (MQ=255)                                                                                                         aaaaGATGGGAGCTTGATCAAAAACAAAAAATATTCCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAgt                                                                                                                                      <  2:601946/251‑1 (MQ=255)                                                                                                                        gATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGa                                                                                                                       <  1:245337/251‑1 (MQ=255)                                                                                                                         aTCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTg                                                                                                                                                                                                         <  2:606497/168‑1 (MQ=255)                                                                                                                         aTCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTg                                                                                                                                                                                                         >  1:606497/1‑168 (MQ=255)                                                                                                                             aaaacaaaaaATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAag                                                                                                                     >  1:757215/1‑248 (MQ=255)                                                                                                                               aacaaaaaATATTTCACTCGACAGGAGGATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGaa                                                                                                                 >  1:806333/1‑250 (MQ=255)                                                                                                                                   aaaaaTATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGTCGGTAAATCAGAAAGGAGAAAACATCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGAGCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTgg                                                                                                            <  1:407580/251‑1 (MQ=255)                                                                                                                                   aaaaaTATTTCACTCGACAGAAGTATTTAAAGTGCTCCCGTTCCGTGGGCTTCGACGGTAAATCAGAAGGGATAAAACCCCTAGAACTACCTACGATCGTAACCGAAACGCAATTACCACTGGCGTCCGTGTTATGGTTAGCGGCACCGGTCCCACTGGCAAGATCCGGTCGATGGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAA‑CCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGc                                                                                                           <  2:283446/251‑1 (MQ=255)                                                                                                                                     aaaTATTTCACGCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATAAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGca                                                                                                          <  1:519723/251‑1 (MQ=255)                                                                                                                                              aCTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGAc                                                                                                 >  1:761207/1‑251 (MQ=255)                                                                                                                                                      ggAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACAACTGGCAGCCGTGTTATGGTTAGc                                                                                                                                                                                                                         >  1:736511/1‑123 (MQ=255)                                                                                                                                                      ggAGTATGTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACTCCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGc                                                                                                                                                                                                                         <  2:736511/123‑1 (MQ=255)                                                                                                                                                                 atTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGAAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATg                                                                               >  1:739768/1‑250 (MQ=255)                                                                                                                                                                   tGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGaa                                                                             <  1:201834/250‑1 (MQ=255)                                                                                                                                                                          gTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCgt                                                                     >  2:458769/1‑251 (MQ=255)                                                                                                                                                                                     ttCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACTGTAACGCAATTCCCACTTGCAGCCGTGTTATGGTTCGCGGCACCGGTCACACGGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGGTTTTGAAGGTTGTGAAGAGAAGCTGGCATCACTGGACCTGATTCGTCTCGGCATGAATTAAGCGTGTGAATgccg                                                           <  1:352104/250‑1 (MQ=255)                                                                                                                                                                                         aCTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAt                                                      >  1:302505/1‑251 (MQ=255)                                                                                                                                                                                          cggaaaaTCAAAAAGGAGAATACACCTACAAAGACTTTCGATTCTAACCGAAACGCAATTACCACAGGCAGCCGGGTTATGGTTAGCTGCCCCGGTCACACTGGCAAGATCCCGTCGATTAATCCTGAAGGTCTGCCCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACCGGACCTTATTCGTCGAGGCATGAACTAAGCGTGTGAATTCCGCCGATg                                                     <  2:713016/247‑1 (MQ=255)                                                                                                                                                                                          cTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCcgcg                                                                                                                                                   >  1:644344/1‑157 (MQ=255)                                                                                                                                                                                          cTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCcgcg                                                                                                                                                   <  2:644344/157‑1 (MQ=255)                                                                                                                                                                                                               acacCTATGACGACCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTc                                <  1:458769/251‑1 (MQ=255)                                                                                                                                                                                                                 cccgatgccgacCTACGTACGTAAGCGTAACGCACGTACCACTGGATGCTGTGTTATGGTTAGCGGCACCGGCCACACTGGCAAGTTCCTGTCGATTGATTCTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACTTGAGTCGTCTCGGCATTAACTAAGCTTGTGAAGGCCGCCGATGGCGGCATTGCTTTTTTACTTCAc                              <  2:454603/243‑1 (MQ=255)                                                                                                                                                                                                                           acCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGCAGGTCTGACCg                                                                                                                                                                     >  1:374540/1‑106 (MQ=255)                                                                                                                                                                                                                           acCTACGATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCg                                                                                                                                                                     <  2:374540/106‑1 (MQ=255)                                                                                                                                                                                                                                 gATCGTAACCGTAACGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACg              >  1:592819/1‑251 (MQ=255)                                                                                                                                                                                                                                            taacGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgc   >  1:109485/1‑251 (MQ=255)                                                                                                                                                                                                                                            taacGCAATTACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCAAGGAATATTTTGCCACGGGCGCTTTcgc   >  2:559603/1‑251 (MQ=255)                                                                                                                                                                                                                                             aactcAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATCGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg  <  1:408054/247‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   CGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTATGACGACCTACGATCGTAACCGTAACGCAATCACCACTGGCAGCCGTGTTATGGTTAGCGGCACCGGTCACACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCG  >  NC_000913/1629003‑1629489

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