Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,823,802	A→G	intergenic (‑33/+47)	recA ← / ← pncC	DNA recombination/repair protein RecA/NMN aminohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,823,802	0	A	G	100.0%	252.4 / NA	71	intergenic (‑33/+47)	recA/pncC	DNA recombination/repair protein RecA/NMN aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/32);  total (39/32)

GAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATCCGGCCCGGCGATACCACTAATAGACA  >  NC_000913/2823558‑2824048                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       gAAGATTCCGGTCCGTAGAATTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAat                                                                                                                                                                                                                                                  <  1:435818/251‑1 (MQ=255)    gATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATg                                                                                                                                                                                                                                               <  1:270703/251‑1 (MQ=255)       tcgggTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCCGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGGCGATAGCCATTTTTACTCCTGTCATGCCGGGCAATACCGGATGGTCAATCTGTTc                                                                                                                                                                                                                                            >  2:338896/4‑251 (MQ=255)        ccgttCCGTATATTTCGACGATACGGCCCATCGTCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCCCATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCt                                                                                                                                                                                                                                           <  1:291680/247‑1 (MQ=255)               gTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAg                                                                                                                                                                                                                                    <  2:449882/251‑1 (MQ=255)                         acgaTACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGTGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACt                                                                                                                                                                                                                          <  1:89285/251‑1 (MQ=255)                           gaTACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGt                                                                                                                                                                                                                        <  2:451343/251‑1 (MQ=255)                           gaTACGGCACATAGGAAGACAACCTGCCACAAGAGCGATATCCAGTGACCGCAAACCGGTAGAGCTGCTTTCCCCATCCATGGAACGGTCTTCACACAGGCGCCTGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGt                                                                                                                                                                                                                        >  2:309134/1‑251 (MQ=255)                              aCGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGGTCTGTGGAAGCAATTATACTGTATg                                                                                                                                                                                                                     <  1:450488/251‑1 (MQ=255)                                ggCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCt                                                                                                                                                                                                                   >  2:5537/1‑251 (MQ=255)                                                    gcaccAAGCGCGCTATCCAGTGAAAGCGAACGTGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCACCGCTTTCTGTTTGTTTTCGTCGATTGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGtt                                                                                                                                                                                                 <  2:144328/246‑1 (MQ=255)                                                    gCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGttt                                                                                                                                                                                                <  1:282886/250‑1 (MQ=255)                                                               gATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATTCCGGATGGTCAATATGTTCTGTTTAAGCAATGATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgt                                                                                                                                                                                    <  1:300956/251‑1 (MQ=255)                                                               gATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCAGTTTTACTCCTGTCATGCCGGGTAATACCGGAGGGTCAATAGGTTCTGTTGAAGCAATGATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAttgt                                                                                                                                                                                    <  1:49155/251‑1 (MQ=255)                                                                      aGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCcac                                                                                                                                                                              >  1:272267/1‑250 (MQ=255)                                                                           aaGCGAACCGGTCGAGATGGTGTACACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATGGTTTCTCAATCTGGCCCAGGTCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGACATTGTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAAGATGTTCTGTTGAGGAAATTATACTGTATGCTCATACAGTAGCAAGTGTTTTGTAGAAAGTTTTGCCACAAGGTc                                                                                                                                                                        <  2:414994/251‑1 (MQ=255)                                                                            aGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCTCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCCCAAGGTCt                                                                                                                                                                       <  2:232130/251‑1 (MQ=255)                                                                             gCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTg                                                                                                                                                                      >  1:384610/1‑251 (MQ=255)                                                                               gAACCGGTAGAGATGGTTTCCACTTCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATGGTTGCCACAATGTCTGCa                                                                                                                                                                    <  2:393263/251‑1 (MQ=255)                                                                               gAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGc                                                                                                                                                                     >  1:309102/1‑250 (MQ=255)                                                                                aaCCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCa                                                                                                                                                                    >  2:550869/1‑250 (MQ=255)                                                                                   cGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGCTGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGc                                                                                                                                                                >  2:203717/1‑251 (MQ=255)                                                                                     gTAGAGATGGTTTCCACATCCATGGCACGGTCTTCACCCCGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGCGCTGCCGCCAACGCTTTATGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCtt                                                                                                                                                              >  2:553869‑M1/1‑249 (MQ=255)                                                                                     gTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCtt                                                                                                                                                              <  2:91681‑M1/251‑3 (MQ=255)                                                                                                  ccACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggccttt                                                                                                                                                 <  2:113968‑M1/251‑16 (MQ=255)                                                                                                     caTCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAAGTGTTGCCACAAGGTCTGCAATGCttgagacttgcctttctg                                                                                                                                              <  2:128975‑M1/251‑19 (MQ=255)                                                                                                          aTGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCCAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgac                                                                                                                                            >  1:281081‑M1/1‑228 (MQ=255)                                                                                                          aTGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctgtctgacaga                                                                                                                                         <  1:203717‑M1/251‑24 (MQ=255)                                                                                                            ggAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagat                                                                                                                                        <  2:527318‑M1/250‑25 (MQ=255)                                                                                                                     ttCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCAAACGCTTTCTGTTCGTTTTGGTCGATCGCCCTTTTTACTCCTGTCACGCCGGGTAATACCGGCTGGTCAATATGTTCTGTTGAAGCACTTCTACTGTATGCTCATACAGTATCAAGTGTTTTGTAGCAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttac                                                                                                                              >  1:551707‑M1/1‑217 (MQ=255)                                                                                                                             ggCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCACAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCctgagactggcctttctgacagatgcttacttactcgcgga                                                                                                                       >  1:574910‑M1/1‑209 (MQ=255)                                                                                                                              gcgcATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaac                                                                                                                     <  1:275938‑M1/251‑44 (MQ=255)                                                                                                                                 cATGATGGAGCATTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGGCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCGGCAATGCgtgagactggcctttctgacagatgcttacttactcgcggaacagc                                                                                                                  <  1:359283‑M1/251‑47 (MQ=255)                                                                                                                                  atgatgGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcg                                                                                                                 >  2:461010‑M1/1‑204 (MQ=255)                                                                                                                                    gatgGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgct                                                                                                               >  2:362329‑M1/1‑202 (MQ=255)                                                                                                                                    gatgGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgct                                                                                                               >  1:554142‑M1/1‑202 (MQ=255)                                                                                                                                     atgGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgctt                                                                                                              >  2:365174‑M1/1‑201 (MQ=255)                                                                                                                                                aCCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagc                                                                                                   >  1:112817‑M1/1‑190 (MQ=255)                                                                                                                                                     aTTGTTTCTCAATCTGGCCCAGTGCTGCCGCCACCGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttcctcgcggaccagcgcttcgttattcagcccctg                                                                                              >  2:552298‑M1/1‑185 (MQ=255)                                                                                                                                                        gTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgt                                                                                           >  2:232872‑M1/1‑182 (MQ=255)                                                                                                                                                           tctcAATCTGGACCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttg                                                                                        >  1:459368‑M1/1‑179 (MQ=255)                                                                                                                                                              cAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgc                                                                                       <  2:24936‑M1/249‑74 (MQ=255)                                                                                                                                                              cAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagccccggcgtttgc                                                                                       <  2:189326‑M1/249‑74 (MQ=255)                                                                                                                                                                  cTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccct                                                                                               >  1:81215‑M1/1‑172 (MQ=255)                                                                                                                                                                  cTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccct                                                                                               <  2:81215‑M1/237‑66 (MQ=255)                                                                                                                                                                   tGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggattt                                                                                >  1:167489‑M1/1‑171 (MQ=255)                                                                                                                                                                    tgcCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGACAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttc                                                                               >  2:512279‑M1/2‑170 (MQ=255)                                                                                                                                                                     gCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTAGATAGCCATTTTTACTCCTGTCATGCCCGGTAATACCGGATGGTCAAAATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGATTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcg                                                                              >  1:238816‑M1/1‑169 (MQ=255)                                                                                                                                                                              tgcCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacgg                                                                     >  2:222531‑M1/1‑160 (MQ=255)                                                                                                                                                                                  gccAACGCTTTCTGTTTGTTTTCGTCGATAGTCATTTTTACTCCTGTAATGCCGGGTAATACCGGATGGTCACTCTGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCAACAAGGTCTGCAATGCttgagactggcctttcttacagatgcttacttactcgcggaacaccgcttcgctattcagcccctgcgtttgcaggatttcgcccaaacggcccc                                                                 >  2:18811‑M1/1‑156 (MQ=255)                                                                                                                                                                                     aaCGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggc                                                              >  2:164722‑M1/1‑153 (MQ=255)                                                                                                                                                                                           ttCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttc                                                          <  2:74018‑M1/249‑103 (MQ=255)                                                                                                                                                                                            tCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaa                                                        <  1:550869‑M1/250‑105 (MQ=255)                                                                                                                                                                                               gtttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctg                                                    >  1:519833‑M1/1‑143 (MQ=255)                                                                                                                                                                                                 ttgtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctga                                                   >  1:373664‑M1/1‑141 (MQ=255)                                                                                                                                                                                                   gtttTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaccggcgcaggccttcaacctgaatc                                                >  2:326731‑M1/1‑139 (MQ=255)                                                                                                                                                                                                      ttCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTCCTACCGGATGGTCAATATGTTCTGTTGAAGCACTTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctgg                                             >  1:539117‑M1/1‑136 (MQ=255)                                                                                                                                                                                                          tcgATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGCCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcga                                          >  1:441979‑M1/1‑132 (MQ=255)                                                                                                                                                                                                             aTAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAAGATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTGGCCACAAGGTCTGCAATGCgtgagactggcctttctgacagatgcttacttactcgcggaacagcgctgcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaaccggaatctggcgaacac                                      <  2:12603‑M1/251‑123 (MQ=255)                                                                                                                                                                                                                  cATTTTTACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttca                                 <  1:5537‑M1/251‑128 (MQ=255)                                                                                                                                                                                                                     ttttACTCCTGTCATGCCGGGTAATACCGGATGGACAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagctgcttacttactcgcggcacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgaaggccttcaacctgaatctggcgaacacgttcacgg                              >  2:347785‑M1/1‑121 (MQ=255)                                                                                                                                                                                                                        tACTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggt                            <  1:362329‑M1/250‑133 (MQ=255)                                                                                                                                                                                                                         aCTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtg                           <  2:309102‑M1/250‑134 (MQ=255)                                                                                                                                                                                                                         aCTCCTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttcggacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggctcaggccttcaacctgaatctggcgaacacgttcacgggtg                           <  2:58527‑M1/250‑134 (MQ=255)                                                                                                                                                                                                                            ccTGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggctaacacgttcacgggtgag                         <  1:431715‑M1/249‑136 (MQ=255)                                                                                                                                                                                                                             cTGTCATGCCGGGTAATACCGGATGGTCACTATCTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGAAATGAttgagactggcctttctgacagatgcttacttcctcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaccggcgcaggccttcaacctgcatctggcgaacccgttcacgggtgaggcc                      >  2:599617‑M1/1‑113 (MQ=255)                                                                                                                                                                                                                              tGTCATGCCGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgttttcaggctttcgcgcaaacggctcaggcattcaacctgaatctggcgaacacgttcacgggtgaggca                      >  2:153213‑M1/1‑112 (MQ=255)                                                                                                                                                                                                                                     ccGGGTAATACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggccaatttcac              >  1:484977‑M1/1‑105 (MQ=255)                                                                                                                                                                                                                                             tACCGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGAATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggccaatttcacgacctac       >  1:359287‑M1/1‑97 (MQ=255)                                                                                                                                                                                                                                                cGGATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggccaatttcacgacctacatct   >  1:102769‑M1/1‑94 (MQ=255)                                                                                                                                                                                                                                                 ggATGGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCttgagactggcctttctgacagatgcttacttactcgcggaacagcgcttcgatattcagcccctgcgtttgcaggatttcgcgcaaacggcgcaggccttcaacctgaatctggcgaacacgttcacgggtgaggccaatttcacgacctacatctt  >  1:117202‑M1/1‑93 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       GAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTGTCATGCCGGGTAATACCGGATAGTCAATATGTTCTGTTGAAGCAATTATACTGTATGCTCATACAGTATCAAGTGTTTTGTAGAAATTGTTGCCACAAGGTCTGCAATGCATACGCAGTAGCCTGACGACGCACCGCATCACGGTCGCCGCTGAAGCATTCCCGCCGGGTAATGCCTTCACCGCGGGCAGTGGCAAAAGCAAACCAGACGGTGCCGACAGGCTTCTCTTCACTGCCGCCATCCGGCCCGGCGATACCACTAATAGACA  >  NC_000913/2823558‑2824048

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