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	94.5%	159.1 / NA	55	E207D (GAA→GAT)	garD	GarD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (21/31);  minor base G (0/3);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTG‑TGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGAAGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCA  >  NC_000913/3275657‑3276141                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 gccgctgtcatcgggggCACCAAAAACCAGCGCGGTATCACCGCCAGCGGCCACTGTG‑TGGCAGGCGTGGTGGACTATGTAGTAAAAATCAGTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAAt                                                                                                                                                                                                                                             <  2:168723/237‑1 (MQ=255)   cGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTG‑TGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGGTTAACGCACCGGCGGCCGTTGTACCTATCCGTCCCATTCACAATATTTCGCTGAATCCTCACTTTGGCCGCGATGTGATgg                                                                                                                                                                                                                                           >  1:347890/1‑251 (MQ=255)    ggtggcatagggggCACCAAAAGCCTGCTCGGTAGCACCACCAACGTCCTATGTG‑TGGCAGGCGTGGTGCACTATGTAGGAAAAATCATGGAAGGCGATCTTCTTCGGAAATACCCGAACGACGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGGCGTGGCG‑TGGCGAGTAACGCACCGGCGGCAGTTGTGCCTAGCGGTACCATACACAATATTTTGCTGAATCCTATCTTTTGCGGCGATTTAATgga                                                                                                                                                                                                                                          <  2:35034/240‑2 (MQ=255)     aTGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTG‑TGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTg                                                                                                                                                                                                                                         <  1:293955/251‑1 (MQ=255)         cAGCGTGGGCACCAAAAACCTGCTCGGTATCCCCACCAGCGTCCACTGTG‑TGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCTGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAAGCCTAACTTTGGCGGCGATGTAATGGTGATTg                                                                                                                                                                                                                                     <  2:369845/251‑1 (MQ=255)          aGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTG‑TGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGGTTAACGCACCGGCGGCAGTTGTACCTATCCGTACCCTTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGGTGTAATgggggttgg                                                                                                                                                                                                                                    >  1:174548/1‑244 (MQ=255)                              gCTCGGTATCACCACCAGCGTCCACTGTG‑TGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACAGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGc                                                                                                                                                                                                                <  1:437247/251‑1 (MQ=255)                                       caccacCAGCGTCCA‑TGTGGTGGCAAGCTTGTTGGACTAGGTAGTAAAAATCATAGATGGCGAGCTGCTTCCGAAATTCCCGAACTTCGATGGCGGGGTGTGGCTGAAGCATTT‑GTACGGTTTTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATGTCGCTGAATCCTAACCTAGGCGGCGAGGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAgaaggat                                                                                                                                                                                                        <  2:17071/250‑7 (MQ=255)                                                gcggcctcttggtgtcAGGCGTGGTGTACTATGTAGTAAAAATCATGTAACGCGTTTTGCTAACGAAATCCCCGAACGTCGATGGCGTGGGGGGGCTGAATCATTT‑GTACGGTTGTTGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCCTTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCttgctg                                                                                                                                                                                               <  2:314236/236‑1 (MQ=255)                                                  ccccctggg‑tgGCTGGCGTGGTGGACTAGGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATGCCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCATTTGTGCCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACt                                                                                                                                                                                            <  2:399129/243‑1 (MQ=255)                                                          g‑tgGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTGCCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACgg                                                                                                                                                                                     <  2:516024/250‑1 (MQ=255)                                                            tgGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCCGCGGCAGTTGTACCTATCCGTAACATTCACAATATTTTGCTGAATCCTAAATTTGGCGGCGATGTAATGGGGATTGGCCTGGGGTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACCGGAACgggt                                                                                                                                                                                   >  1:551856/1‑249 (MQ=255)                                                                         ggACTATGTAGTAAAAGTCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTa                                                                                                                                                                      <  1:361237/251‑1 (MQ=255)                                                                          gACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTTAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGGTGTGAAAAGTTGCAACCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTAt                                                                                                                                                                     >  1:433004/1‑251 (MQ=255)                                                                                gtagtaAAAAGCATTGAACGCGATCTGATACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGt                                                                                                                                                               <  2:55393/251‑1 (MQ=255)                                                                                       aaaTCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGGAGAAAgc                                                                                                                                                        <  1:354022/251‑1 (MQ=255)                                                                                           cATTGAACGCGAGCGGCTACCGACAGACCTGAGCGTCGATGGCGGGGTGGGGCCTAATCATTT‑GTACGGTTGTGTCG‑TGGCGATGAACGCTCCGGCGGCAGTTGTACCTATCCGTACCATGCACAATATTTCGCTGAATCCTAGCTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAATTTGCAGCCTGAGCTCGTTCTGACTGTAACTGATGATGTGCAAGCTATTCCAGTacgaggcgcc                                                                                                                                                     <  2:280685/250‑10 (MQ=255)                                                                                             ttGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATGAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACTGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAg                                                                                                                                                   <  1:313102/250‑1 (MQ=255)                                                                                                     gATCTGCTACCGAAATACCCGACCGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGgat                                                                                                                                                                                   <  2:422814/210‑1 (MQ=255)                                                                                                     gATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGGGATTTGCCTGGGTTGTGAAAAGTTGCAGCCTGGGCGCCTGCTGACTGGAACGgat                                                                                                                                                                                   >  1:422814/1‑210 (MQ=255)                                                                                                      aTCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTTGCGGCGGTGTAATGGGGATTGGCCTGGGGTGTGAAAAGGTGCAGCCTGGGCGCCTGCTGACTGGGAAGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCCGCATTGTCAGt                                                                                                                                         >  1:228040/1‑251 (MQ=255)                                                                                                             ccctgaatgcgcGAACGTCGATGGCGGGGTGGGGCTGAATCATTTCGT‑CGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGGTGTACCTATCCGTTCCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAgg                                                                                                                                  <  2:173899/240‑1 (MQ=255)                                                                                                             aCCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTTCCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTgta                                                                                                                                    <  2:504667/249‑3 (MQ=255)                                                                                                             aCCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCa                                                                                                                                    <  1:149852/249‑1 (MQ=255)                                                                                                                  aaTACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGAGGTTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAAAATTTTGCCGAATCCTAACTTTGGCGGCGATGTTATGGTGATTGGCCTGGGTTTTGGAAAGTTGCAGCCTGAGCGCCCTCTGACTGGAACGGATGATTGGCAAACTATTCCCATAGAAAGCGCCAGCATTGTCAATTTGCAGGATGa                                                                                                                              >  1:376173/1‑250 (MQ=255)                                                                                                                   aTACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGaaa                                                                                                                            <  1:494006/251‑1 (MQ=255)                                                                                                                             gTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGGTGTACCTATCCCTACCATTCAAAAAATTTTGCTGAATCCTAAATTTGGCGGCGATGTAATGGAGGTTGGCCTGGGTTGTGAAAAGTTGCACCCTGAGCGCCTGCTGACTGGAACGGTTGATGTGCAAGCTATTCCAGTAGGAAGCGCCAGCATTGTCAGTTTGCAGGGTGAAAAGccagtcgg                                                                                                                  >  1:570560/1‑244 (MQ=255)                                                                                                                                   ggCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATTCGTACCATTTACAATATTTCGCTGAATCCTAACTTTGGCGGCGGTGTAAAGGGGATTGGCCTGGGGTTTGAAAAGGTGCAGACTGAGCGCCTGCTGACTGGGACGGGTGATGTGCAAGCTATTCCAGTAGAAAGCGCCCGCATTGTCAGTTTTCAAGATGAAAAGGATGTC‑GGTTTTct                                                                                                           >  1:283429/1‑250 (MQ=255)                                                                                                                                   ggCGTGGTGGGGCTGAATCAATT‑GTACGGTTGTGGCG‑TGGCGATGAACGCTCCGGCGGCAGTTGATCCTTTCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGATAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCa                                                                                                           <  2:178895/251‑1 (MQ=255)                                                                                                                                       gggtggGGCAGAATCATTT‑GTTCGGTTGTGGCG‑TGGCGATTAACTCACCGGCGGCAGTTTTTCCTATCCGTACCATTCACAATATTTCGCTGAATTCTAACTTTGGCGGCGATTTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAATGCGCCAGCATTGTCAGGTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTcc                                                                                                       <  1:520843/250‑1 (MQ=255)                                                                                                                                           ggggCTGAATCATTT‑GTACGGTTGTGGCT‑GGGCGATTAACGCACCGGCGGCAGTTGTACCTAGCCGTACCATGCACAATATGTCGCTGAATCCGAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTT‑GGTTTTCAGTCCATgg                                                                                                   <  2:347890/251‑1 (MQ=255)                                                                                                                                                tGAATCATTT‑GTACGGTTGTGGCG‑TGGCGCTTTAAGGACCGGCGGCAGCTGTAGCTATTCGTACCATTAACAATAATTACCTGAATTCTAACTTTGGCGGCGATGTCAAGGGGATTGGCCTGGGGTGTGAAAAGTTGAAGGCTAAACGCCTGCTGACTGGAACAGTTGATGAGCAAGCAGTTACAGTAGAAAAGGGCAGCATTGTCAGGTTGCAAGATGAAAAGCATGTCAGGtttcttgcccagggcg                                                                                                 >  1:506946/1‑236 (MQ=255)                                                                                                                                                 gAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGa                                                                                               <  2:174548/249‑1 (MQ=255)                                                                                                                                                       atgtggacggtttgtgcg‑tggcgATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGACTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATAtt                                                                                        <  2:166353/236‑1 (MQ=255)                                                                                                                                                          tttACGGTTGTGGCGGTGG‑GATTAACGCACCGGCGGCAGTTTTACCTATCCGTACCAGTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGAGGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCTGCATTGTCAGTTTGCAGGATGAAAATCATGTC‑GGTTTTCAGTCCATGGGCGAGTATAtttt                                                                                      <  2:285197/248‑1 (MQ=255)                                                                                                                                                                   gtggcg‑tggcGATTAACGCACCGGCGGCATTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGcc                                                                            <  2:228040/251‑1 (MQ=255)                                                                                                                                                                       cg‑tggcgATTAACGCACCGGCGGCAGTTGTTCCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTGTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTTCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAAc                                                                        <  2:517000/251‑1 (MQ=255)                                                                                                                                                                       cg‑tggcgATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAAc                                                                        >  1:398380/1‑251 (MQ=255)                                                                                                                                                                                  aaCGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGccttcaaa                                                               >  1:237057/1‑244 (MQ=255)                                                                                                                                                                                    cGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACaa                                                             >  1:550335/1‑250 (MQ=255)                                                                                                                                                                                           gcggcAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGaa                                                     <  1:547330/251‑1 (MQ=255)                                                                                                                                                                                                  ttGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACgg                                              >  1:444525/1‑251 (MQ=255)                                                                                                                                                                                                    gTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGc                                             <  1:553393/250‑1 (MQ=255)                                                                                                                                                                                                       ccTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCg                                         <  1:370198/251‑1 (MQ=255)                                                                                                                                                                                                         tATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCg                                         >  1:70998/1‑249 (MQ=255)                                                                                                                                                                                                              gttccATTCCCAATATTGCGCTGAATCCTAAGTGTGGCGGCGATGTAATGGTGAGTGGCCTGGGTTGGGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATTATGTGCAAGCTATTCCAGGAGAAAGCGCCAGCATTGTCAGTTTGCATGATGAAAAGCATGTC‑GGTTGTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACGGAATCAACGGCAGCGAGAAAcc                                  <  2:318538/248‑1 (MQ=255)                                                                                                                                                                                                                  cATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGcc                              >  1:400183/1‑251 (MQ=255)                                                                                                                                                                                                                       acaATATTTCGCTGAATCCTAACTTTGGCGGCGCTGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGGGAAACCTGCCCGGCt                         >  2:526146/1‑251 (MQ=255)                                                                                                                                                                                                                            atTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCATGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTCTTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGACCATGGTCGAGGATATTTTACAGATCGCCGAACGCCATCTACAAAAAATGCATCAACGACAGCGAGAAACCTGCCCGGCTTccg                     >  2:542287/1‑248 (MQ=255)                                                                                                                                                                                                                              ttCGCTGAATCCTAACTGTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCAGGGTCGAGGATATTTTGCAGATCTTATAACGCCATCTACAAAAACTGAATCAATGGCAGCGAGAAACCTGCCCGGCTTCAGaa                   <  1:344139/250‑1 (MQ=255)                                                                                                                                                                                                                                        ccTAACTTTGGCGGCGATGTAATGGAGATTGGCCTGGGTTGTGAAAAGTTGCAGCATGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGCAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTAAGACCATGGTCGCGGATATTTTGAAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCGGCCCGGCtttgagacagtttgatgg        >  2:443024/1‑235 (MQ=255)                                                                                                                                                                                                                                           accTTTGGCGGCGATGTAATGGTGATAGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGGGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTAGTTGGTAt     >  1:284575/3‑251 (MQ=255)                                                                                                                                                                                                                                            aCTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGCGAAACCTGCCCGGCTTCAGCACTGGTCGTTGGTATg    >  1:254346/1‑251 (MQ=255)                                                                                                                                                                                                                                               ttGGCGGCGATGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCa  <  1:569506/250‑1 (MQ=255)                                                                                                                                                                                                                                               ttGGCGGCGAGGTAATGGTGATTGGCCTGGGTTGTTAAAAGTGGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGATTGCAGGATGAAAAGCAAGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCa  <  2:433004/250‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 GCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTG‑TGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTT‑GTACGGTTGTGGCG‑TGGCGATTAACGCACCGGCGGCAGTTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGAAGTAATGGTGATTGGCCTGGGTTGTGAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTC‑GGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTTGGTATGCA  >  NC_000913/3275657‑3276141

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