Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,175,827	C→T	noncoding (74/76 nt)	thrT →	tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,175,827	0	C	T	98.6%	222.3 / NA	70	noncoding (74/76 nt)	thrT	tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (36/33);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGACGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGG  >  NC_000913/4175582‑4176065                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                cGAGACGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCActt                                                                                                                                                                                                                                            <  1:390618/251‑1 (MQ=255)  gagaCGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCGCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCActtc                                                                                                                                                                                                                                           <  1:36964/251‑1 (MQ=255)         acacgacgaactaccgtactaaTAACGATAGCTGCTCTCGATATGCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTGCCAAGCTGATGATGCGGGTTCGATTCCCGATGCCCGTTCCAAGAGGTGCTGATATAGCTCAGTTGGTAGAGCTCACCCTTGGTAAGGGTGAGGTCGGCAGGTCGAATCTGCCTATCAGCTTCACTTCTTTtc                                                                                                                                                                                                                                      <  2:396209/229‑1 (MQ=255)               cGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTccc                                                                                                                                                                                                                               <  1:373919/250‑1 (MQ=255)               cGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGGGAGGTCGGCAGTTCGGATCTGCCTATCAGCATCACTTCTTTTCTCCTcccc                                                                                                                                                                                                                              >  2:275049/1‑250 (MQ=255)                           cGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTc                                                                                                                                                                                                                  >  2:360663/1‑251 (MQ=255)                             aaTAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGCATCTGCCTATCAGCATCACTTCTTTtc                                                                                                                                                                                                                                      >  1:396209/1‑229 (MQ=255)                                 aaGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCAGAGTTCCAGGATGCGGGCATAGAATAATGGCTATTACCTCAGCCTGCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTACGGTCGGCAGTTCGAATCTGCATATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTATGttt                                                                                                                                                                                                             >  2:374684/1‑250 (MQ=255)                                          tctctcGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTc                                                                                                                                                                                                    >  2:403058/1‑250 (MQ=255)                                               cGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcaaca                                                                                                                                                                                                <  1:197019/249‑1 (MQ=255)                                                 atatTAAGTGTATAATGTAAAACAGGTCGCCGATGTCCAGGATGCGGGCATCGTATAATGGCTATTACCCCAGCCTTCCAAGCTGATTATGCGGGTTCGATTCCCGCGGCCCGCTCCAAGATGTGCTGATATTGCTCAGTTTGTAGAGCGCACCCTTGTTAAGTGTGAGGTCGGCAGTTCTAATTTGCCTATCAGCATCTCTTCTTTTCTCCGCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTc                                                                                                                                                                                            <  2:384667/251‑1 (MQ=255)                                                   atgtaGTGCAGAATGAAAATCAGGTAGCCTAGTTACAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAATAAGTCg                                                                                                                                                                                           <  2:143218/246‑1 (MQ=255)                                                    ttCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCggg                                                                                                                                                                                         >  1:469847/1‑251 (MQ=255)                                                      cAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGc                                                                                                                                                                                        <  1:192545/250‑1 (MQ=255)                                                          gCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGtt                                                                                                                                                                                   <  1:360663/251‑1 (MQ=255)                                                           cAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTg                                                                                                                                                                                  >  1:540155/1‑251 (MQ=255)                                                                   aaaTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTg                                                                                                                                                                          <  1:116175/251‑1 (MQ=255)                                                                          gTAGCCGAGGTCCAGGATTCGGGCATCTTATAATGGTTATTACCGCCGCTTTCCAATCTTATGATGCGTGTTCGATTCCCGCTGCCTGCTTCAAGATGTGCTGATATAGTTCAGATGGTAGCGCGCACCCTTTGTAAGGGTGAGGTCGGCATTTCTAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTGTTTTCCTTCTTTGTATTGCATTCAACAAGTCGGGCATGTTGCCTGCTTTATgtgg‑t                                                                                                                                                                  <  2:384580/251‑1 (MQ=255)                                                                            aGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtgg‑tg                                                                                                                                                                 >  1:297072/1‑250 (MQ=255)                                                                             gCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATGCCCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTACAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGat                                                                                                                                                               <  2:11255/251‑1 (MQ=255)                                                                               cGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGatat                                                                                                                                                             >  1:526198/1‑251 (MQ=255)                                                                                 aGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTAACGCTGCCCGCTCCCAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGAGTGAGGTCGGCGGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACCAGTCGGGCATGTTGACTGTTTGCTGTTG‑TGATATcc                                                                                                                                                           >  2:486191/1‑250 (MQ=255)                                                                                     ccAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTtct                                                                                                                                                                                                                                     >  1:181484/1‑174 (MQ=255)                                                                                     ccAGGAGGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCCCCCTTTGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTGtct                                                                                                                                                                                                                                     <  2:181484/174‑1 (MQ=255)                                                                                      cAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGAt                                                                                                                                                      >  1:222074/1‑251 (MQ=255)                                                                                        ggATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTTCTGATATAGCTCAGTTGGGAGAGCGCACCCTTGGTACGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCCTTCAACAAGTCGGGCCTGTTTCCTGGTTGATGTGG‑TGATTTCACCGAttt                                                                                                                                                    >  2:317622/1‑251 (MQ=255)                                                                                         gATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGAttt                                                                                                                                                    >  1:432852/1‑250 (MQ=255)                                                                                            gCGGGCAACGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCAACCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATcc                                                                                                                                                <  1:371152/251‑1 (MQ=255)                                                                                             cGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCg                                                                                                                                               >  1:298417/1‑251 (MQ=255)                                                                                              gggCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCgt                                                                                                                                              >  1:278789/1‑251 (MQ=255)                                                                                                         aaTGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAggg                                                                                                                                   <  1:110030/251‑1 (MQ=255)                                                                                                          aTGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGTTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGa                                                                                                                                  <  2:298417/251‑1 (MQ=255)                                                                                                          aTGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGGGCTGATATAGATCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTATAgggc                                                                                                                                  >  1:289403/1‑250 (MQ=255)                                                                                                                       cAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCTCCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCt                                                                                                                      <  1:317622/250‑1 (MQ=255)                                                                                                                       cAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTc                                                                                                                       <  2:222074/249‑1 (MQ=255)                                                                                                                               cAAGCTGCTGATACGGGTTAGATTCCCGCTGCCCGCTCCAAGATGGGCTGATATAGCTCGGTTGGTCGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTCTTGCATTCAACCAGCCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACACTCGAGGTCTCAAgcaa                                                                                                               >  2:569919/1‑246 (MQ=255)                                                                                                                               cAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGaaa                                                                                                               >  1:70812/1‑249 (MQ=255)                                                                                                                                  gcggaggatgCGGGTTCGATTCCCGCTGACCGATACAAGATGTGCTGATATAGCTCAGTTGGTATAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATTACTTCTTTTCTCCTCCCTGTTTTTTCCTGCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGttt                                                                                                          <  2:262061/245‑1 (MQ=255)                                                                                                                                     gatgatGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTg                                                                                                         >  2:262451/1‑249 (MQ=255)                                                                                                                                                   gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCCGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGCCAATCGATGTCTAAAGAAAAGTTGGAACGTACAAAAcagc                                                                                          >  2:486091/1‑247 (MQ=255)                                                                                                                                                                cgttccatgagttgCTGATTTATCTCTGTTGGTATAGCGCACCCTTGGTAAGTGTGAGGTCGGCAGTTCGCATCTGTCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTGTTCCTTCTGTTTATTGCATTCAACAAGTCGGTCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTCTCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAACGTTTGAACGTACAAAACCGCACGTTAACGTCGGt                                                                            <  2:309091/239‑1 (MQ=255)                                                                                                                                                                  cTCCACGATCTGCTGATATAGATCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGCGGCCGGCAGTTCGAATCTGCCTATCAGCCTCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCAACGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGACAGTACAAAACCGCACGTTAACGTCGGTAc                                                                          >  2:236012/1‑251 (MQ=255)                                                                                                                                                                   tCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTTTGTTTATTCCATTCACCAAGTCGGCCTTGTTGCCTGGTTGA‑GTGGGTGATATCACCGATTTTTCCGTGTCTTAAAGGGACAATGGTTGTCTAAAGAAAGGTTTGAAAGTACAAAACCgacggtaaacaccggaaat                                                                         >  2:310098/1‑233 (MQ=255)                                                                                                                                                                      aaGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATc                                                                      <  1:236012/251‑1 (MQ=255)                                                                                                                                                                       aGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCTAATCTGCCTATCAGCATCACTTCTTTTCGCCTCCCTGTTTTTTTCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCg                                                                     <  2:540155/251‑1 (MQ=255)                                                                                                                                                                             gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACg                                                               >  1:299377/1‑251 (MQ=255)                                                                                                                                                                                gATATATCTCCGTTGGTAGAGCGGACTCGTGGTAAGGGTGAGGTCGGAAGTTCGAATCTGCCTATCAGCATCTCTTCTTTTCTCCTCCCTGTGTTTTCCTTCTGTTTATTGCCTTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGGGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGTACTTTAACGTCGGTACTATCGGCCACGtt                                                             <  2:395490/250‑1 (MQ=255)                                                                                                                                                                                gATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTg                                                            <  1:374684/251‑1 (MQ=255)                                                                                                                                                                                gATATAGCTCAGTTGGTAGAGAGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATAAACATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTTG‑TGATATCACCGAGTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTAACGGCCACGttt                                                            >  2:369643/1‑250 (MQ=255)                                                                                                                                                                                 atatAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGAATCACTTCTTTTCTCCTCCATGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGGACTATCGGCCACGTTg                                                            >  2:32035/1‑250 (MQ=255)                                                                                                                                                                                 atatAGCTCAGTGGGTAGATCGCACCCTTGGTAAGGGTGATGTCTGCAGTTCTAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGa                                                           <  2:427246/251‑1 (MQ=255)                                                                                                                                                                                  tataGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGAc                                                          >  1:509329/1‑251 (MQ=255)                                                                                                                                                                                   ctaGCTCAGTTGGTCGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTATTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATATTGCCTGGTTGATGTGG‑TGCTATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGCAAATTTTGAACGTACAAAACCTCCAGTTAACGTCGGTACTATCGGCCACGTTAAcc                                                         >  2:428135/2‑251 (MQ=255)                                                                                                                                                                                    taGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGTTGTCTAAAGaa                                                                                                                >  2:17062/1‑195 (MQ=255)                                                                                                                                                                                    taGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGaa                                                                                                                <  1:17062/195‑1 (MQ=255)                                                                                                                                                                                           gatGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTaa                                                  >  1:93613/3‑250 (MQ=255)                                                                                                                                                                                                 agagCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGT‑GGATATCACCGATTTATCCGTGTCTTAGAGGGCCAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGGAATATCGGCCACGTTGACCATGGTAAaacaacc                                           >  2:60371/1‑250 (MQ=255)                                                                                                                                                                                                 agagCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAaacaac                                            <  1:60371/250‑1 (MQ=255)                                                                                                                                                                                                 agagCGCACCCTTGGTAAGGGAGAGGGCGGAAGTTCGAATCTGCATATTAGCATCACTTCTTTTCGCCTCCCTGTTTTTTCCTTCTGTTTATTGCAGTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCAACGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGGCTGTACTATCGGCCACGTTGACCATGGTAAaacaac                                            <  2:506428/250‑1 (MQ=255)                                                                                                                                                                                                  gagCGCACCATTGGTAAGGGTGAGGTCGGCAGTTCGAATATGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGCCAATCGATGTCTATAGAAAAGTTTAAACGTACAAAACCGCACGTTAACGTCGGTACTAAAGGACCCGTTGACCTTGGTAAAAGAACgg                                          >  2:503439/1‑250 (MQ=255)                                                                                                                                                                                                          ccTTGTTAAGGTGGAGGTCGGCAGTTCGAATCTGCCTATCAGCTTCACTGTGTTTCTCCTCCCTGTGTTTTCCTTCTGTTTATTTCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGGGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCt                                  <  2:45580/251‑1 (MQ=255)                                                                                                                                                                                                                tAAGGGTGAGGTCGGCATTTCTAATCGTCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGGTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAAt                             <  2:390261/250‑1 (MQ=255)                                                                                                                                                                                                                 aaGGGTGAGTGCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTGCGCCTCCCTTTTGTTTCCTTCTGTTTATGGCATGCAACAAGTCGGGCATGTTGCCTGGTGGATGTGG‑TGATATCACCGATTTATCCGTGGCTTAGAGGGACAATCGATGTCTAAAGAAAATTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTTACCATGGTAAAACATCGTGGACCGCTGCAATCa                           <  2:526198/251‑1 (MQ=255)                                                                                                                                                                                                                  aGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCAc                          >  1:323630/1‑251 (MQ=255)                                                                                                                                                                                                                   gggTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACt                         <  1:32035/251‑1 (MQ=255)                                                                                                                                                                                                                      tGAGGTCGGCAGTTCGAATCTGCATATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATAACCGCTTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCCCGTTGACCATGGGAAAACAACGCTGCCCGATGCCATCACTAcc                      >  2:241343/1‑251 (MQ=39)                                                                                                                                                                                                                       gAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTAcc                      <  1:275049/250‑1 (MQ=255)                                                                                                                                                                                                                             ggCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAACAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTgg               >  1:215523/1‑251 (MQ=255)                                                                                                                                                                                                                                          tGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCACCAATTCGGACATGTTGCCTGTTTGATGAGG‑TGATATCACCGCTTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACAGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAAAGCTGACCGCTGCACTCACTACCGTACTGGCTAAAACCtcagg  >  2:291636/1‑247 (MQ=255)                                                                                                                                                                                                                                           gCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTAcg   <  1:159073/249‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                CGAGACGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGG‑TGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGG  >  NC_000913/4175582‑4176065

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