Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,275,902	A→T	E207D (GAA→GAT)	garD →	GarD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,275,902	0	A	T	96.5%	176.0 / ‑5.4	58	E207D (GAA→GAT)	garD	GarD
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (23/32);  minor base G (2/0);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGAAG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCAGTGC  >  NC_000913/3275657‑3276145                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAAt                                                                                                                                                                                                                                                  <  1:229068/251‑1 (MQ=255)                ggCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGGTTAACGCCACGGCGGCAGGTGTAACTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATT‑TAATgggggtttggatcgg                                                                                                                                                                                                                                   >  1:587328/1‑238 (MQ=255)                ggCACCAAAAACCTGCGCGGTATCACCACCAGCGTCCACTGTGGGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGGGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTg                                                                                                                                                                                                                                          <  2:587328/244‑1 (MQ=255)                         aaCCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCCATCCGT‑ACCATTCACCATATTT‑CGCTGAATCCTAACTTTTGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑aaa                                                                                                                                                                                                                         >  1:201772/1‑250 (MQ=255)                                           accaGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAGATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCGAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGccc                                                                                                                                                                                                      <  2:479341/251‑2 (MQ=255)                                           accaGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGAGATTTACGCACCGGGGGCAGGTGGTCCTATTCGT‑ACCCTTCAAAATATTTTCG‑TGAATCCTAACTTTGGCGGCGATG‑GAATTGGGATTGGCCT‑GGGGTGTG‑AAAAGTTGCAAGCTGAGCGccc                                                                                                                                                                                                      >  1:531557/1‑250 (MQ=255)                                             caGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGGGGTTAAAGCACCGGCGGCAGTTGTACCTATCCGTTACC‑TTCTAAAAATTT‑TGCTTGATCCTAAATTTGGCGGCGATG‑TACTGGGGGTTGGGCT‑GGGGTGTG‑GAAAGGTGCAGGCTGAGCGCctgc                                                                                                                                                                                                    >  1:350023/1‑251 (MQ=255)                                                   ccACTGTGTGGCAGGCGTGGTGGGCTATGTAGTAAAAATCATTGACCGCGATCTGATACCGAAATACCCGAACGTAGATGGCGTGGTGGGGCTGAATAATTTGTACGGTTGTGGCGTGGCGATTAACGCACCCGCGGCAATTTTTCCCTTCCCC‑ACCATTCTCAATATTT‑TGCTGAATCCCAAATTTTGGGGCGGGGGTAA‑GGGGATTGGCCT‑GGGGTGTG‑GAAAGTTGccaccttaagggctgctgat                                                                                                                                                                                                >  2:386327/1‑230 (MQ=255)                                                   ccACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGAGCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATGCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GTGTTGTG‑AAAAGTGGCAGCCTGAGCGCCTGCTGACTg                                                                                                                                                                                              <  2:353392/251‑1 (MQ=255)                                                        gtgtgGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGaa                                                                                                                                                                                           <  1:93817/249‑1 (MQ=255)                                                                           atATGTAGTAAAAATCATTGAACGCGATCTGCTTCCGAAATACCCGAGCGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTc                                                                                                                                                                      <  2:574528/250‑1 (MQ=255)                                                                             aTGTAGTAAAAATCATTGAACGCGATCTGCTATCGAAATTCCCGACCGTCTATGGCGTGGTTGGGCTGAATCATTTGTACGGTTGTGGTGTGGCGATTAACGCTGCGGCGGCAGTTGTTCCTTTCCGT‑ACCATTCACAATATGT‑CGGTGAATCCTAACTTTGGGGGCGATG‑TAATTGTGATTGGTCT‑GGGTTTTG‑AAAAGCTGCAGGCTGAGCGTCGGCGGACTGGAACTGATGATGTGCAAGCTATTcc                                                                                                                                                                     <  2:525638/250‑1 (MQ=255)                                                                             aTGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCa                                                                                                                                                                    <  2:83748/251‑1 (MQ=255)                                                                               gtagtaAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAAt                                                                                                                                                                                                                                                  <  2:374868/173‑1 (MQ=255)                                                                               gtagtaAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGGACCTATTCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAAt                                                                                                                                                                                                                                                  >  1:374868/1‑173 (MQ=255)                                                                                  gtaAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGACGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGGGCTTTACGCACCGGCGGCAGGTGTACCTATCCGG‑ACCCTTTCCAATCTTT‑CCCTGGATCCTAACTTTTGGGGAGTGG‑TAATGGTGATTGTGCT‑GGGGTGTG‑AAAAGTTTACGCCTGAGCGCCTGCTGGCTGGGACGGAGGATGTGGCAGCTAtttcggtagg                                                                                                                                                               >  2:358943/1‑243 (MQ=255)                                                                                       aagcATTGAACGCGAGCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAgcg                                                                                                                                                          <  1:225789/248‑1 (MQ=255)                                                                                         tCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTTACGCACCGGCGGCAGTTGTACCTATCCGT‑GCCATTCACAATATTT‑CACTGAATCCTAACTTTGGGGGGGTAG‑TAATGGTGGATGGCCTTGGGTT‑TG‑AAAAGTTGCAGCCCGAGCGCCCGCTGAATGGAACGGATGATGTGCAAGCTATTCTAGTAGAAAGCGcc                                                                                                                                                        >  2:119275/1‑251 (MQ=255)                                                                                          cATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTTAAGCACCGGCGGCAGTTGTAACTATCCGT‑ACCCTTCCCCATATTT‑CCCCTAATCCTTACTTTGGCGGCGATG‑TTATGGTGTTTGGGCT‑GGGTTGGGGAAAAGTT‑CAGCCTGGGCGCCTGCTGACTGGAAAGGATGGTGTGCAAGAGATTCCCGGAGAAAGGGccc                                                                                                                                                       >  2:529977/1‑250 (MQ=255)                                                                                              gAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCAt                                                                                                                                                   <  2:339378/251‑1 (MQ=255)                                                                                               aaCGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATATGCAAGCTATTCCAGTAGAAAGCGCCAGCAtt                                                                                                                                                  <  1:511738/251‑1 (MQ=255)                                                                                                    gatatGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGGTTAACGCACCGGCGACAGGTGGACCTATCAGT‑ACCAGTCACAATATTT‑CGCTTAATTCTAAATTAGGCGGGGGTG‑TAATGGGGGTTTGCCT‑GGGGTGTG‑AAAAGGTGCAGCATTAGCGCCTGCTGACTTGAACGGCGGATGTGCAAGAGATTCCCGTAGGAAGCGCCAACATTggccg                                                                                                                                             >  2:22823/5‑247 (MQ=255)                                                                                                      tCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCCGTTGTACCTATCCGT‑AACATTCACAATCTTT‑CGCTGAATCCTAACTTTTGCGGGGGTG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGGCTGGAACGGGTGATGTGCAAGCTATTCCAGTAGAAAGCGCCCGCATTGTCAGtt                                                                                                                                           >  1:44765/1‑251 (MQ=255)                                                                                                       cTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCCCCGGCGGCAGATGTACCTATCCGT‑ACACTTACCAATATTT‑CGCTGAATCCTAAATTTTGCGGAGATG‑GAATGGGGATTGGCCT‑GGGGTGTG‑GAAAGTTGCACCCTGTGCGCCTGCTGATTGGAACGGGTGGTGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGttt                                                                                                                                          >  1:569474/1‑251 (MQ=255)                                                                                                            aCCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAg                                                                                                                                      <  2:43510/250‑1 (MQ=255)                                                                                                            aCCGAAATACCCGAACGGCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAgg                                                                                                                                     <  1:249761/251‑1 (MQ=255)                                                                                                                aaaTACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGa                                                                                                                                 <  1:266458/251‑1 (MQ=255)                                                                                                                 aaTATCCGAACGTCGGTGGCGTGGTGTGGCTGAATCATGTGTACGGTTGGTGCGCGGCGATTAACGCTCCGGCGGCAGTTGTACCTAGCCGT‑AGCATTCACAATATTT‑CGCTGAATCCTAACGTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTTAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCAAGCATTGTCAGTTTTCAGTATGaa                                                                                                                                <  2:587651/251‑1 (MQ=255)                                                                                                                           cGTCGATGGCGTGGTGGTGCTGAATCGTTTGTACGGTTGGTGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCg                                                                                                                      <  2:7525/251‑1 (MQ=255)                                                                                                                                aTGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGCTTAACGCAACGGCGGCAGTTTTTCCTATCCGT‑ACCATTTCCAAAATTT‑TCCTGAATCCTAACTTTGGCGGCGGTG‑TAATAGGGATTGGCCT‑GGGGTGTG‑GAAAGTTGCAGCCTGAGCGCCTGCTGGCTGGAAAGGATGATGTGGAAGCTATTCCCGGAGAAAGAGCCCGCATTGTTAGTTTTCAGGATGAAAAGAATGTCGGtttt                                                                                                                 >  2:421678/1‑251 (MQ=255)                                                                                                                                 tGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCTCCAGCATTGTCAGTTTGCAGGATGa                                                                                                                                 <  2:32462/234‑1 (MQ=255)                                                                                                                                       ggtggGGCTGAATCATTTGTACGGTTGTAGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATAGTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCa                                                                                                          <  1:501180/251‑1 (MQ=255)                                                                                                                                                  aTCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGat                                                                                               <  1:421678/251‑1 (MQ=255)                                                                                                                                                    cATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGatat                                                                                             <  2:310699/251‑1 (MQ=255)                                                                                                                                                    cATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGatat                                                                                             <  1:538742/251‑1 (MQ=255)                                                                                                                                                             ggTTGTGGCGTGGCGATTAACGCACCGGCGGCAGGTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGc                                                                                        >  1:182648/1‑247 (MQ=255)                                                                                                                                                                tgtgGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCg                                                                                  <  1:60193/250‑1 (MQ=255)                                                                                                                                                                      gtggcgATTAACGCAACGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCCGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGGTCGCCGAACg                                                                           >  2:72861/1‑251 (MQ=255)                                                                                                                                                                                 cGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCCTGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACaa                                                                 >  1:413508/1‑250 (MQ=255)                                                                                                                                                                                 agCACCGGCGGAAGATGTACCTATCCGT‑AACATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TTATGGTTATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGCACGGATGATGTTCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGCTGAAAAGCATGTCGGTTTTAAGTCCCTGGTCGAGGTTATTTTGCAGATCGCCGAACGCCACCTACaa                                                                 >  1:413512/2‑250 (MQ=255)                                                                                                                                                                                       ggcggcAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTg                                                           >  1:258091/1‑250 (MQ=255)                                                                                                                                                                                            cAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACGTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAAttaa                                                     <  1:248306/251‑4 (MQ=255)                                                                                                                                                                                              gTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACg                                                   <  1:487819/251‑1 (MQ=255)                                                                                                                                                                                                  tACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGc                                                 >  2:8469/1‑249 (MQ=255)                                                                                                                                                                                                       aTCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑GAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACGGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCAGGTCGGTTTTCAGTCCATGGGCGAGGATATTTTGCAGGTCGCCGAACGCCATCTACAAAAACTGAATCATCGGCAGCga                                            <  2:531557/249‑1 (MQ=255)                                                                                                                                                                                                            gaaCCATTAAAAATCTTT‑CGCTGAATCCTACCTTTGGCGGCGATA‑TAATGGTGATTTGCCT‑GGGATGTG‑AAAAGTTTCAGCCTGAGCGCCTGCTGACTGGAAAGGGTGATGTTCAACCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGCAAAGCAAGCATGTTTTCAGTCCGTGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAACAACTGAATCAACGGCAGCGCGAAAcc                                      >  2:13034/3‑251 (MQ=255)                                                                                                                                                                                                              aCCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTAc                                                                   <  2:168119/220‑1 (MQ=255)                                                                                                                                                                                                              aCCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTAc                                                                   >  1:168119/1‑220 (MQ=255)                                                                                                                                                                                                                 aTTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGGTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATGCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGccc                                 <  2:232425/251‑1 (MQ=255)                                                                                                                                                                                                                           ttt‑CGCTGAATCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCAGACTGGAACGGATGATGTGCAAGCTATTCCAGAAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATATTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAAACTGCCCGGCTTCAGa                        >  1:89370/1‑250 (MQ=255)                                                                                                                                                                                                                                     aTCCTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGtt              >  2:230020/1‑251 (MQ=255)                                                                                                                                                                                                                                       ccTAACTTTGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTgg            >  1:566437/1‑251 (MQ=255)                                                                                                                                                                                                                                       ccTAACTGTGGCGGCGATG‑TAATGGTGATTGGCCT‑TGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCTGTTTTCAGTCCATGGTCGAGGATATTTTGCAGAGCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTgg            <  1:230763/251‑1 (MQ=255)                                                                                                                                                                                                                                          aaCTTTGGCGGGGATG‑TAATGGTGATTGTCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTTCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTAt         <  2:479532/251‑1 (MQ=255)                                                                                                                                                                                                                                             tttGGCGGCGATG‑TAATGGGGATTGGCCT‑GGGTTGTG‑ACAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTCTTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGCTGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCTTTGGTATGc       >  1:92170/1‑250 (MQ=255)                                                                                                                                                                                                                                              ttGGCGGCGATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCAGTGTCAGTTTGCAGGATGAAAAGCATGGCGGTTTTCAGTCCATGGTCGAGGATATTGTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCa      <  1:386327/250‑1 (MQ=255)                                                                                                                                                                                                                                                 gcggcgATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCAGTGc  >  2:127480/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gcggcgATG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCAGTGc  <  1:119275/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGT‑ACCATTCACAATATTT‑CGCTGAATCCTAACTTTGGCGGCGAAG‑TAATGGTGATTGGCCT‑GGGTTGTG‑AAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCAGTGC  >  NC_000913/3275657‑3276145

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