Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	330,640	A→G	E393G (GAG→GGG)	betT →	choline:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	330,640	0	A	G	100.0%	116.3 / NA	40	E393G (GAG→GGG)	betT	choline:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/25);  total (15/25)

GTCCGGTTGAGTGGATGAATAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAACGGGATATCCGCGCTG  >  NC_000913/330394‑330884                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     gTCCGGTTGAGTGGATGAATAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGGGGAAGCGATGGTCCATCCGGGgcgc                                                                                                                                                                                                                                                  >  1:43662/1‑251 (MQ=255)   ccgggtgggTGGATTAATTACTGGACGGTCTTTTTCTGGGGATGGTGGGGGGGATGGTCGCCGTTTGGTGGCTTGTTCCTGGGGCGTATCTCGCGGGGGCATACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTTACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCgg                                                                                                                                                                                                                                                <  1:90105/243‑1 (MQ=255)     ggTTGAGTGGATGAATAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTGCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCt                                                                                                                                                                                                                                              <  1:182942/251‑1 (MQ=255)                       cTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGgcgc                                                                                                                                                                                                                            >  1:50911/1‑251 (MQ=255)                         ggACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGgcgc                                                                                                                                                                                                                            >  1:210731/1‑249 (MQ=255)                                  tcttctGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGc                                                                                                                                                                                                                  <  2:225001/250‑1 (MQ=255)                                      ctGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGCTGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGttt                                                                                                                                                                                                              >  2:208618/1‑250 (MQ=255)                                                        atggtggttgattGTCGTCTTTTTCCTTGCGGTTATCTCGCGTGGGCTTACCATTCGCCAGTTCGTGCTGGGGCTGTTGATTATTCCGTTTACCTTCACGCGTTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGAGATGGTCCAGCCGGGGCGCGGCTTCTCCAGCCTGCTGGCGCAGTAACCGGCGTTTACCTTTAGCGCCTCCggcg                                                                                                                                                                                           <  2:2448/240‑4 (MQ=255)                                                        aTGGTCGCCGGTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCATGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAAAAGAGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCAGCGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTAAAGCCTGCAGGCGAAGTATCCGGCGTTTGCCTTTAGCGCCACCGTCg                                                                                                                                                                                           >  2:207193/1‑251 (MQ=255)                                                              gCCGTTTTTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTTTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGTATTTGCCGAGGAAGCGATGGTTCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGccacc                                                                                                                                                                                      <  1:158888/250‑1 (MQ=255)                                                                ggTTTGTCGGCTGGTTCCTGGGGCGTATCTCGCGGGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTGCCTTCACGCTGGTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAAGCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCAtt                                                                                                                                                                                   <  2:237727/250‑1 (MQ=255)                                                                    tGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACt                                                                                                                                                                                <  1:233229/250‑1 (MQ=255)                                                                    tGTCGGCTGGTTCCTGGCGCGTATCTCGCGTGGGCGTTCCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTCCCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACt                                                                                                                                                                                <  2:61234/250‑1 (MQ=255)                                                                     gTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCTCGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATACACTGGGGCGCGGCATTTGCCGATGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTg                                                                                                                                                                               <  1:254025/250‑1 (MQ=255)                                                                             gtacatgacgcgTATATCGAGTGGGTGGACCATTCGCCAGTTCGTGCTGAGCACGATGATGATTCAGTTTTCCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCACGCTGTATGAAATCATCCACGGCGGCGCGGCATTTACCGGGGAAGCGATGGTAAATCCGGGGCGCGGCTTCTAAAGCCTGCTGGCGCAGTATCCGGCGTTTACATGTAGCGCCTCCGTAGCCACCATTAAAGGCCTGATGt                                                                                                                                                                      >  2:150465/9‑251 (MQ=255)                                                                                              cgcgTGGGCGTCCCATTCGCCAGCTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGGTATGGCTCTCTGTGTTCGGCAATAGCGCGCTGTAAGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCTGGGCGCGGCTTCGACATCCTGCTGGCGCAGTATCCGGCGTTTACTTTTAGCGCCTCCGTCGCCACCATTGCTGGCCTGCTGGTTGATGTGACCTcg                                                                                                                                                        <  2:164685/248‑1 (MQ=255)                                                                                                 gTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGAc                                                                                                                                                  >  1:192410/1‑251 (MQ=255)                                                                                                      cGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCgg                                                                                                                                             >  1:33299/1‑251 (MQ=255)                                                                                                               cGCCAGTTCGTGCTGGGCACGGTGATTATTCCGTTTGCCTTCACGCTGTTATGGCTCTCGGGGTTCGGCAATAGCGCGCGGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGGCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTgg                                                                                                                                     <  2:246619/250‑1 (MQ=255)                                                                                                                 ccggTTCGTGCGGGGCACGTTGATTATTCCGGTTACCTTCGCGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAAGCAGCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGGCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGc                                                                                                                                  <  2:128639/248‑1 (MQ=255)                                                                                                                      tCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTggg                                                                                                                              <  1:183057/250‑1 (MQ=255)                                                                                                                      tCGTGCTGGGCACGTTGATTATTCCGTGTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCTTCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTTCCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTTCTggg                                                                                                                              <  1:28034/250‑1 (MQ=255)                                                                                                                          gCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAt                                                                                                                          <  1:262375/250‑1 (MQ=255)                                                                                                                             gggCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAAttt                                                                                                                        <  1:235471/249‑1 (MQ=255)                                                                                                                                   gTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGc                                                                                                                >  1:82126/1‑251 (MQ=255)                                                                                                                                      gATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGCCTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGc                                                                                                             >  2:15928/1‑251 (MQ=255)                                                                                                                                             ccGTTTCCCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGa                                                                                                      <  2:210731/251‑1 (MQ=255)                                                                                                                                                        cACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAACCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGCCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCg                                                                                           <  2:193731/251‑1 (MQ=255)                                                                                                                                                              gTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGcc                                                                                      <  2:50911/250‑1 (MQ=255)                                                                                                                                                                                   aaaTAGCGCGCTGGATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGtcttct                                                                <  2:49724/250‑1 (MQ=255)                                                                                                                                                                                                         atccgggggggcgcggcATTGGCCGAGGCAGCGGTGGACCGGCCGGGGTGGGGCTCCTACAGCTCGCTGGCGCAGAACCCGGGGGTTTCCTTTAGCGCCTCCGTCGCCACCATGACGGGGTTGCTGTTTTATGTGGCCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATAGCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTTGTCGGTGGCGATTGGCctgct                                          <  2:113936/242‑1 (MQ=255)                                                                                                                                                                                                          tcCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGCATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTTGTCGGTGGCGATTGGCctgctg                                         >  2:238505/1‑251 (MQ=255)                                                                                                                                                                                                            cACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGAc                                       <  1:49297/251‑1 (MQ=255)                                                                                                                                                                                                             ccggcggGGCGGCATTTGCCGAGGAAGCGAGGGTCCATCCGGGGCGCGGATTCTACAGCCTGCTGGCGCAGTTTCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACg                                      <  2:145084/250‑1 (MQ=255)                                                                                                                                                                                                                            ttGCCGAGGAAGCGATGGTCCATCCCGGGCGCGGCTTCTACTGCCTGCTGGCGCAGTTTCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACGGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTGCTCCGCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCTATTGGCCTGCTGACGCTCGGCATGCtgat                        <  2:187938/250‑1 (MQ=255)                                                                                                                                                                                                                               ccGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACt                    >  1:81399/1‑251 (MQ=255)                                                                                                                                                                                                                                cGAGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACtc                   >  1:35728/1‑250 (MQ=255)                                                                                                                                                                                                                                  aGGAAGCGATGGTCCATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAAc                 <  1:87097/251‑1 (MQ=255)                                                                                                                                                                                                                                                cATCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAACGGGATATCCGCGCt   >  1:205919/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aTCCGGGGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGCCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGACTTTCACCCCGCAGCTTAAAGATATCAACAGCGACGCCCCCCGCTGGCTCCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAACGGGATATCCGCGCTg  >  1:124212/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GTCCGGTTGAGTGGATGAATAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTCGCCACCATTACTGGCCTGCTGTTTTATGTGACCTCGGCGGACTCCGGGGCGCTGGTGCTGGGGAATTTCACCTCGCAGCTTAAAGATATCAACAGCGACGCCCCCGGCTGGCTGCGCGTCTTCTGGTCGGTGGCGATTGGCCTGCTGACGCTCGGCATGCTGATGACTAACGGGATATCCGCGCTG  >  NC_000913/330394‑330884

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