Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	40,450	A→G	8.9%	F494F (TTT→TTC)	caiT ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,450	0	A	G	8.9%	137.7 / 4.5	45	F494F (TTT→TTC)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (23/18);  new base G (2/2);  total (25/20)
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.97e-01

TCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTG  >  NC_000913/40164‑40727                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                      tccaCACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGTTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTtaataac                                                                                                                                                                                                                                                                                     >  1:507057/1‑288 (MQ=255)          cggcgAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCg                                                                                                                                                                                                                                                                             <  2:436914/288‑1 (MQ=255)           ggcgAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGt                                                                                                                                                                                                                                                                            <  1:399072/288‑1 (MQ=255)                  gTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTAGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATaa                                                                                                                                                                                                                                                                    >  1:81867/1‑289 (MQ=255)                           cccGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACga                                                                                                                                                                                                                                                           >  2:491990/1‑289 (MQ=255)                             cacgTCCAGACCACCGTGCCCTTCAGGGACAAGCAGACGGTCGATACCCATAGCGGCCGGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACgaag                                                                                                                                                                                                                                                         <  1:248141/287‑1 (MQ=255)                                                  ttCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCg                                                                                                                                                                                                                                    >  2:127653/1‑289 (MQ=255)                                                            agcagACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCgg                                                                                                                                                                                                                          >  1:39214/1‑289 (MQ=255)                                                                             cATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAg                                                                                                                                                                                                          >  1:51554/1‑288 (MQ=255)                                                                               tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCCTCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAACATCGATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTCATAAAGGAGAGCGTCACCATAATGTTGACGAAGACCAGCGGGTATCGTCCGGTGATAATGGCGGTTTGAACCGGTTTCAGgc                                                                                                                                                                                                        >  1:207082/1‑288 (MQ=255)                                                                               tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAgg                                                                                                                                                                                                         >  1:138540/1‑287 (MQ=255)                                                                               tatCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGgc                                                                                                                                                                                                        >  1:324302/1‑288 (MQ=255)                                                                                      cAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATGAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAgc                                                                                                                                                                                                >  2:98160/1‑289 (MQ=255)                                                                                      cAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATGAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAgc                                                                                                                                                                                                >  2:98157/1‑289 (MQ=255)                                                                                          tcTTTGACCAAACGTTCCCGGTAGAGGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCAGTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAgcgc                                                                                                                                                                                              <  1:345571/286‑1 (MQ=255)                                                                                                   aaaCGTTCCGGTTAGACGCTGTCACGGTCGCACTAGGCAAAATAGGCCTTCCAGTTTTCGCTGGCCATCAGTTCGCGGATACTGGCGACAAAATGTTCCTGAGCATCATTTAAATTAAAATCCATCTTTCAAAATCTTGATATTTTTGGGGTGAATGAATCTTTCTAGTTCTGTTTCTCGTCTTAAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACATCGGGCATCCTCCGGCTTTACTGGCGGTTTGCATCGGTTTCTGTCCGACGTGCGCTCGCAGCACaa                                                                                                                                                                                   >  1:310352/1‑289 (MQ=255)                                                                                                          ccGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGAt                                                                                                                                                                            >  2:48420/1‑289 (MQ=255)                                                                                                                  aCGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCa                                                                                                                                                                     >  2:87242/1‑288 (MQ=255)                                                                                                                                      ggCCAAATAGGCCTCCCAGTTTGTGCGGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAGTTAAAATCCATCCTTCAACGTCTTGATATTTTGGAGGTGAATCAATCTGTCCAGTGCTGGTTCTCGTCTGTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGGATCGGTTTCAGTCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGGTa                                                                                                                                                <  1:292112/289‑1 (MQ=255)                                                                                                                                      ggCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCTCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATa                                                                                                                                                <  2:138540/289‑1 (MQ=255)                                                                                                                                      ggCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCAAACAAGAATTGACCAACCCAt                                                                                                                                                 >  2:459775/1‑288 (MQ=255)                                                                                                                                                gCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACcagca                                                                                                                                      >  2:493166/1‑289 (MQ=255)                                                                                                                                                     ccAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAggt                                                                                                                                 <  2:81867/289‑1 (MQ=255)                                                                                                                                                          tttCACTGGCCATCAGTTCCCGGATAACGGCGACAAACAGTTACTGCTAATCATTTAAATTAAAATCCATCTTTCAGCCTCTTGATATTTTGGGGGTTAATTAATCGTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATTTAGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCTCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAggtgg                                                                                                                               <  1:522079/286‑1 (MQ=255)                                                                                                                                                              gCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATACAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGttcttc                                                                                                                         >  1:221050/1‑288 (MQ=255)                                                                                                                                                                      tCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATcgcg                                                                                                                <  2:39214/289‑1 (MQ=255)                                                                                                                                                                                             aaCAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGtt                                                                                                                             >  1:220279/1‑253 (MQ=255)                                                                                                                                                                                             aaCAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGtt                                                                                                                             <  2:220279/253‑1 (MQ=255)                                                                                                                                                                                              aCAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGc                                                                                        <  1:87242/289‑1 (MQ=255)                                                                                                                                                                                              aCAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGc                                                                                        <  1:491990/289‑1 (MQ=255)                                                                                                                                                                                                aGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATGAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCa                                                                                      <  1:98157/289‑1 (MQ=255)                                                                                                                                                                                                aGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCa                                                                                      <  1:127653/289‑1 (MQ=255)                                                                                                                                                                                                aGTTCCTGCGCATTATTTAAATTAAAAGCCATCTTTCAACCTTCTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATGAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCa                                                                                      <  1:98160/289‑1 (MQ=255)                                                                                                                                                                                                   tCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCa                                                                                      <  2:51554/286‑1 (MQ=255)                                                                                                                                                                                                                      taaaATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAAc                                                                <  2:324302/289‑1 (MQ=255)                                                                                                                                                                                                                                    cAACCTCTTGATATTTTTGGGGTTAATTAATCTGTCCAGTTCTGTTTCGCGTCGTTAAGAAAGGAGGGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTTCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAAAAATACCGATAATACCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCAACATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTgg                                                  <  2:310352/289‑1 (MQ=255)                                                                                                                                                                                                                                          ctTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATaa                                            >  2:137592/1‑289 (MQ=255)                                                                                                                                                                                                                                                     gggggTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATg                                 >  1:267895/1‑289 (MQ=255)                                                                                                                                                                                                                                                                    tttCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACAtg                  >  2:316520/1‑289 (MQ=255)                                                                                                                                                                                                                                                                             cTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTg         <  1:48420/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                tttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGc                                                                          <  2:419314/221‑1 (MQ=255)                                                                                                                                                                                                                                                                                tttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGc                                                                          >  1:419314/1‑221 (MQ=255)                                                                                                                                                                                                                                                                                tttCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTg      >  2:293370/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                     cgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCTATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGCGGTGGTTCTTCACCATCGCGTACCTCGCGGCAACTCGACATCCCCCGGGTATACGAGCAGGCGTTACCCAGCGTCACGGTGGCCATAAAGCAGAGGATGAAGAATCCCCACCTGGTGGCGGCGCTGAGTg  >  1:414039/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                     cgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTg  >  2:452038/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                     cgTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTg  >  2:264436/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                      TCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGGGTTAATTAATCTTTCCAGTTCTGTTTCGCGTCTTTAATAAAGGAGAGCGTCACCATAATGTTGACGAAGAACAGCGGGCATCCTCCGGCGATAATGGCGGTTTGAATCGGTTTCAGGCCGCCGAGCGCCAGCAGAACAATACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTG  >  NC_000913/40164‑40727

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