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	98.5%	213.6 / NA	65	E207D (GAA→GAT)	garD	GarD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (26/38);  minor base G (0/1);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGAAGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCAGTGCGG  >  NC_000913/3275657‑3276147                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gCCGATGGCAGCGTGGGCACCAAAAGCCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAAt                                                                                                                                                                                                                                                   <  2:82589/251‑1 (MQ=255)       ggCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAAGCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGAGGTAATGG‑TGAt                                                                                                                                                                                                                                            <  2:65163/251‑1 (MQ=255)                  cACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGtt                                                                                                                                                                                                                                 <  1:115273/251‑1 (MQ=255)                  cACCAAAAACCGGCTCGTTATCACCTCCAGCGTCCTCTGTGTGGCAGGCGTGGTGGACTAGGGAGTAAAAATCATTGAACGGGAGCTGCTACCGAAATACCCGATCGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCTTGGGtt                                                                                                                                                                                                                                 <  2:131902/251‑1 (MQ=255)                                     aTCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGa                                                                                                                                                                                                              <  1:49948/251‑1 (MQ=255)                                       caccacCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAACTCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCCGTTGTATCTATCC‑GTACCATTTACAATATTTCGCTG‑CATCCTAACCTTGGCGGCGATGTCATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAgg                                                                                                                                                                                                            >  1:44367/1‑250 (MQ=255)                                       caccacCAGCGACCACTGTGTGGAAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGATACCGAAATACCCGCACGTCGATGGCGTGGAGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGGGGCAGTTGTACCTATTC‑GGAACCTTTACAATATTTTGCTG‑TATACTAACTTTGGCGGCGGTGGAATGG‑GGGTTGGGCTGGGGTGGGAAAAGGTTCAGCCGGAgc                                                                                                                                                                                                            >  2:292041/1‑251 (MQ=255)                                              aGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGCACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGGGATTTACGCCCCGGCGGCAGTTGTTACTTTCC‑GCACCATTACCAATATTTCCCTG‑AATCCCCACATTGGCGGCGATGTCATGG‑TGATTGGGCTGGGTTGTGAAAACTTACCGCCTCAGCGgctgct                                                                                                                                                                                                     >  1:67542/1‑251 (MQ=255)                                                        gggtgGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGAGTGCGTGGGGGGGCTGAATCAGTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCGATCC‑GTTCCATGCACTATATTTCGCTT‑AATCCTAACTTTGGCGGCGATGTAATGT‑TGATTGGGCTGGGTTGTGAAAAGTGGCAGCCGGAGCGCTTGCTGGCTGGAAc                                                                                                                                                                                            <  2:232651/248‑1 (MQ=255)                                                            ggCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGgatg                                                                                                                                                                                       <  2:459802/251‑1 (MQ=255)                                                              cAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGgatgat                                                                                                                                                                                     <  1:460622/251‑1 (MQ=255)                                                                ggCGTGGGGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTTAAGCACCGGCGGCAGTTCTACCTATCCCGT‑CAATTAACAATATTTGGCTGTAATC‑TAACTTTTGCGGCGATGTAATGG‑TGATTGGACAGGGTTGGGTACAGCTGGACCCTGAGCGCCTGCTGGCTAGGACGgatgatat                                                                                                                                                                                   >  2:5615/1‑249 (MQ=255)                                                                 gCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGAtgtg                                                                                                                                                                                  <  1:142791/251‑1 (MQ=255)                                                                   gtggtgGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCa                                                                                                                                                                                <  1:251410/251‑1 (MQ=255)                                                                         gACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTAt                                                                                                                                                                          <  1:176441/251‑1 (MQ=255)                                                                           cTATGTAGTAAAAATCATTGCACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGGTTAACGCAACGGCGGCAGTTTTACCTTTCC‑GTAACATTCAAAATATTTCGCTG‑AGTCCTAAATTTTGCGGCGATGTAATGGGTG‑TTGGGCTGGGGTGTGGAAAGGTGGAACATGAGCGCCTGCTGACTGGCACGgagtagtgggcagctattt                                                                                                                                                                        >  2:329301/1‑234 (MQ=255)                                                                             aTGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAAATTTGGCGGCGATGTAATGG‑GGATTGGCCTGGGTTGTGAAAAGGTGCAGCCTGAGGGCCTGCTGACTGGAACGGATGATGTGGAAGCTATTcc                                                                                                                                                                       >  1:353369/1‑250 (MQ=255)                                                                                          cATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTggg                                                                                                                                                                                                                                   <  1:300080/177‑1 (MQ=255)                                                                                                aCGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGGTTAACGCACCGGCGGCAGTTGTACCTATCC‑GGACCATTCCAAATATTTCGCTG‑AATCCTTACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGGTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGGAAGGATGATGGGCAAGGTATTCCAGTAGAAAAAGCCAGCAttt                                                                                                                                                   >  1:366823/1‑250 (MQ=255)                                                                                                     aTCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑CTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑GGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGGACGGGTGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAgg                                                                                                                                              >  1:129430/1‑250 (MQ=255)                                                                                                         gCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGttt                                                                                                                                            <  1:405361/249‑1 (MQ=255)                                                                                                         gCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGAGGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTg                                                                                                                                           <  1:189473/250‑1 (MQ=255)                                                                                                            aCCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCa                                                                                                                                         <  1:288741/249‑1 (MQ=255)                                                                                                            aCCGAAATACCCGAACGGCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCGGACTGGAACGGATGATGTGCaa                                                                                                                                                                               <  1:132368/211‑1 (MQ=255)                                                                                                                         aaCGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGt                                                                                                                          <  2:250624/251‑1 (MQ=255)                                                                                                                               gATGGCGCGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGTTTAAACCACCGGGGGCCGTTGTACCTATCC‑GTACCCCTCACAATCTTTCGCTG‑AATCCTAACCTTGGCGGCGCTGTAATGG‑GGATTGGCCTGGGGTGTGAAAAATTGCACCCTGAGCGCCTGCTGGCTGGGACGGATGGTGTGCAAGCTATTACCGTAGAAAGCGCCAGCATTGTCAGGTTTCAGGATGAAACGCATGTCgggtt                                                                                                                    >  2:397193/1‑248 (MQ=255)                                                                                                                                  ggCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCa                                                                                                                 <  1:65229/251‑1 (MQ=255)                                                                                                                                  ggCGTGGTGGGGCTGAATCATGTGTACGGTCGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTAGCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGATAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCa                                                                                                                 <  1:179494/251‑1 (MQ=255)                                                                                                                                  ggCGTGGTGGGGCGGAATCATTTGTTCGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATGTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATTG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAATCGCCAGCATTGGCAGTTTGCAGTATGAAAATCATGTCGGTTTTc                                                                                                                  <  2:332689/250‑1 (MQ=255)                                                                                                                                          ggggCTGAATCATTTGTACGGTTGTGGCGTTGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATgg                                                                                                         <  2:210210/251‑1 (MQ=255)                                                                                                                                           gggCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAAAGCACCGGCGGCAACTTTAACTATCC‑CGACCATTCCCAATATTTCGCTG‑AAACCTAACTTTTGCGGCGGTGTAAAGG‑TGATTGGCCTGGGGTTTGGAAAGGTGCCGCCCGGGCGCCCGCTGCCTGGAACGGGTGATGTGCAAGGGCTTCCAGGAGAAAACGGCCGCCTTTTCAGTTTGCAGGATGAAAAGCATGGCGGGTTTTAGTccctagg                                                                                                        >  2:5843/1‑246 (MQ=255)                                                                                                                                            ggCTGAATCATGTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTc                                                                                                       <  2:353499/251‑1 (MQ=255)                                                                                                                                             gCTGAATCATTTGTACGGTTGTGGCGTGGCGGTTAACGCACCGGCGGCAGGTGTACCTATTC‑GTACCATTCACAATATTTCGCTG‑AATCCTAAATTTGGCGGGGATGTAATGG‑GGATTGGCCTGGGTTGTGAAAAGGTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCCTTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGGCCATgggcg                                                                                                      >  1:388012/1‑248 (MQ=255)                                                                                                                                             gCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCCGCGGCCGTTGTACCTATCC‑GTACCATTCACAAAATTTCGCCG‑CATCCTAACTTTGGCGGCGCTGGAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGGAGAAAGCGCCAGCATTATCCGTTTGCAGGATGAAAAGGCTGTCGGGTTTTAGTCCCTggg                                                                                                        >  1:204571/1‑248 (MQ=255)                                                                                                                                                 aaTCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGa                                                                                                  <  1:160277/251‑1 (MQ=255)                                                                                                                                                                      gtggcgATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGACAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATCTTTTGCAGATCGCCGAACg                                                                             >  1:166249/1‑251 (MQ=255)                                                                                                                                                                         gcgATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTAGCAGATCGCCGAACGcc                                                                           >  1:248123/1‑250 (MQ=255)                                                                                                                                                                              tAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTa                                                                      <  1:431494/250‑1 (MQ=255)                                                                                                                                                                                aCGCACCGGCGGCAGTTGTACCTATCC‑GTTCCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAAAGG‑TTATTGGCCTGGGGTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGGCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTGTGTAGAGCGCCGAACGCCATCTACaa                                                                   <  2:185978/251‑1 (MQ=255)                                                                                                                                                                                  gCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACaaaa                                                                 <  1:253926/251‑1 (MQ=255)                                                                                                                                                                                   cACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATGTCGCTG‑AATCCTAACTTTGGCGGCTATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGGCGAGGATATTTTGCAGATCGCCGAACGCCATCTACaaaaa                                                                <  2:146377/251‑1 (MQ=255)                                                                                                                                                                                              gTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACg                                                     <  1:61638/251‑1 (MQ=255)                                                                                                                                                                                                 gTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCa                                                  <  2:361036/251‑1 (MQ=255)                                                                                                                                                                                                       aTCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCgag                                             <  1:292041/250‑1 (MQ=255)                                                                                                                                                                                                          c‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCATCGAGAAAc                                         >  1:6412/1‑251 (MQ=255)                                                                                                                                                                                                            gTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAAcc                                        >  1:451322/1‑251 (MQ=255)                                                                                                                                                                                                                 aTTCACAATATTTCGTTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCCGccc                                   >  2:387004/1‑251 (MQ=255)                                                                                                                                                                                                                   tCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCg                                  <  1:5615/250‑1 (MQ=255)                                                                                                                                                                                                                    cacaATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCgg                                 <  1:304732/250‑1 (MQ=255)                                                                                                                                                                                                                     acaATATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑AGATTGGCCTGGGTTGTGCAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCt                               >  2:120786/1‑251 (MQ=255)                                                                                                                                                                                                                       aaTATTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGCTTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTc                             >  2:222819/1‑251 (MQ=255)                                                                                                                                                                                                                         tatTTCGCTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTc                             >  1:194346/1‑249 (MQ=255)                                                                                                                                                                                                                              cGCTG‑AATCCTAACTTTGGCGGCGATGTAAGGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCAGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTg                      <  2:245553/251‑1 (MQ=255)                                                                                                                                                                                                                                cTG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACg                                                                             <  1:459189/194‑1 (MQ=255)                                                                                                                                                                                                                                cTG‑AATCATAACTTTGGCGGCGATGTAATGG‑TGATTGGCATGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGCTGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACg                                                                             >  2:459189/1‑194 (MQ=255)                                                                                                                                                                                                                                 tG‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTc                   >  1:119012/1‑251 (MQ=255)                                                                                                                                                                                                                                  g‑AATCCTAACTTTGGCGGCGATGTAATTG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCg                  <  2:353369/251‑1 (MQ=255)                                                                                                                                                                                                                                  g‑AATCCTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGt                    <  1:397193/249‑1 (MQ=255)                                                                                                                                                                                                                                   gaagccTAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGGTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCg                  <  1:293580/247‑1 (MQ=255)                                                                                                                                                                                                                                    aaTCCTAACTTTGTCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGGGACTGCTGACGGGAACGGATGATGTGCAAGCTAATCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCg                  <  2:230/250‑1 (MQ=255)                                                                                                                                                                                                                                         tAACTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGAGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGAAATTTTGAAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAAATGGACGTTGGt             >  2:464688/1‑250 (MQ=255)                                                                                                                                                                                                                                         tAACTTTAGCGGCGATGTAATGG‑TGATTTGCCTGGGATGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACTAAAACTGAATCAACGGCAGCTAGAAACCTTCCCGGCTTCAGAACTGGTCGTTGGt             >  2:46004/1‑250 (MQ=255)                                                                                                                                                                                                                                          aaCTTTGGCGGCGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTa            <  1:120786/250‑1 (MQ=255)                                                                                                                                                                                                                                                ggcggcGATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCAGTg     >  2:37093/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ggcgATGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCAGTGCgg  >  1:453935/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAGTTGTACCTATCC‑GTACCATTCACAATATTTCGCTG‑AATCCTAACTTTGGCGGCGAAGTAATGG‑TGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCAGTGCGG  >  NC_000913/3275657‑3276147

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