Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,311,337	C→T	A236T (GCA→ACA)	alsB ←	D‑allose ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,311,337	0	C	T	100.0%	193.9 / NA	59	A236T (GCA→ACA)	alsB	D‑allose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/30);  total (29/30)

CGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGCGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTCACCACCT  >  NC_000913/4311092‑4311578                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  cGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTg                                                                                                                                                                                                                                              <  1:107430/251‑1 (MQ=255)  ggCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGc                                                                                                                                                                                                                                             <  2:327784/251‑1 (MQ=255)                      tgTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTACTTTATCCAGCGGGATAACCTTGCCGTATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGCTATCCGCCGGGTTCTGGGCAACCGTCGCGGTCCTTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGATTGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGCGCAACACCCATTGcc                                                                                                                                                                                                                          >  2:497612/2‑250 (MQ=255)                         gACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCACCACCCATTGCCATCg                                                                                                                                                                                                                      >  2:314296/1‑251 (MQ=255)                            cAGGATTGAATCGACCAGTTTAAATTCCGGAGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCgtgt                                                                                                                                                                                                                   >  1:429834/1‑251 (MQ=255)                               gATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCAt                                                                                                                                                                                                                >  1:485721/1‑251 (MQ=255)                                  tGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTc                                                                                                                                                                                                              <  1:537367/250‑1 (MQ=255)                                  tGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTc                                                                                                                                                                                                              <  1:148121/250‑1 (MQ=255)                                  tGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTc                                                                                                                                                                                                              <  2:290269/250‑1 (MQ=255)                                          ccAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAg                                                                                                                                                                                                     <  1:193726/251‑1 (MQ=255)                                              tttAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAAt                                                                                                                                                                                                       >  1:575466/1‑245 (MQ=255)                                                     tCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTaa                                                                                                                                                                                           <  1:195912/250‑1 (MQ=255)                                                           gctttgtcacgcgtgacaatctttccGGATTTCTCAGCGTCAACCATCAGCCTCAGACACTTTGCGCCGTTATCCGCCTGGTTCTGGGAAACCGTCGCGGTCATTTGTCAGGCTTCCACCATTTTGCGGGCGTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCTTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTTACATTCGCGCAATGGATCGCTTTAATATTCga                                                                                                                                                                                    <  2:13851/227‑2 (MQ=255)                                                              tgagccAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATACGCTTCAGACCCGTTGCGCCGATATCCCCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTCCCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGTCTGAGCATCACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGAt                                                                                                                                                                                  <  2:447482/246‑1 (MQ=255)                                                                aTCCAGCGGGATACCCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGACGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCCTCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGTTCGCTTTAATATTCGGATTAc                                                                                                                                                                               >  2:249404/1‑251 (MQ=255)                                                                  ccAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGATCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACg                                                                                                                                                                              <  1:436292/250‑1 (MQ=255)                                                                    aGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTg                                                                                                                                                                           >  1:342470/1‑251 (MQ=255)                                                                          aTAACCTTCACGGATTTCTCAGCGTCAACCATCAGCTACACACCCGTTGAGCCGATATCTGCCGCGTTCTGCCAAAACGTCGCGGACATCTGTCCGACTTCCACCCTTTTGCGGGCTTCCGGAATGCCATCTGTACCGAAGACCAGCACTTTTCCCGTTTTTCCGACGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGATTTATTCTTCGGATTACGttcacacac                                                                                                                                                                     >  2:298961/1‑244 (MQ=255)                                                                          aTAACCGTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAacac                                                                                                                                                                     <  1:101182/251‑1 (MQ=255)                                                                                   ccGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGc                                                                                                                                                            <  1:460827/251‑1 (MQ=255)                                                                                    cGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTTGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCg                                                                                                                                                           >  1:481540/1‑251 (MQ=255)                                                                                     ggATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCg                                                                                                                                                           >  1:42822/1‑250 (MQ=255)                                                                                           ctcAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATcc                                                                                                                                                     >  1:273921/1‑250 (MQ=255)                                                                                              aGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTg                                                                                                                                                 >  1:542623/1‑251 (MQ=255)                                                                                                 gTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCAGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCtt                                                                                                                                              >  1:330712/1‑251 (MQ=255)                                                                                                     aCCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAAt                                                                                                                                          >  1:150111/1‑251 (MQ=255)                                                                                                     aCCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAAt                                                                                                                                          <  1:104508/251‑1 (MQ=255)                                                                                                           aGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTc                                                                                                                                    >  1:148511/1‑251 (MQ=255)                                                                                                             cTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCCTCGTGTCATTCGCGCAATAGTTCGCTTTAATATTCGGATTACGTTGCAACACGTTAGAGGCGACATCCACTGCTTTAATGCGGTccc                                                                                                                                  >  2:478799/1‑251 (MQ=255)                                                                                                                       gTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCTCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGACTGCCATCTGTACCGACGACCAGCCCTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCACCACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTACTATTCTGATTACGTTTCAACACGTTAGTAGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAgg                                                                                                                        >  2:211052/1‑251 (MQ=255)                                                                                                                            gcCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAggctggc                                                                                                                   >  2:529776/1‑251 (MQ=255)                                                                                                                                    ccgccgGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGAc                                                                                                            <  1:104329/250‑1 (MQ=255)                                                                                                                                       ccgGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGGTTGTGACTGTCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGc                                                                                                        <  1:529776/251‑1 (MQ=255)                                                                                                                                         gggttgtgggcacatgTCGCGGTCATTTGTACGGCTTACACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCCCTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATAGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGCTTACGTTGCAACACGTTAGTGGCAACCTCCAGTGCTCTAATGCGGTCCCAGTCGGCAGCCTGGCTGGCGACAAGCtc                                                                                                      >  2:554891/16‑250 (MQ=255)                                                                                                                                                gggCAACCGTCGCGGTCATTTGTCCGGCTGCCACCATTTTGCGGGCTTACGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGATAAGCTTGATCTg                                                                                                <  1:497612/250‑1 (MQ=255)                                                                                                                                                   cAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGgctt                                                                                            >  1:305733/1‑251 (MQ=255)                                                                                                                                                         tCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCtttt                                                                                      >  2:216369/1‑251 (MQ=255)                                                                                                                                                          cgcgGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCACCACCCATTGCCATCGTGTCATTAGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGATAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCttttt                                                                                     >  1:543816/1‑251 (MQ=255)                                                                                                                                                          cgcgGTCATTTGTCCGGCTTCCACCATTGTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATGCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCttttt                                                                                     <  2:383127/251‑1 (MQ=255)                                                                                                                                                                       ccGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCgg                                                                         <  1:278054/250‑1 (MQ=255)                                                                                                                                                                        cGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTg                                                                       <  2:481540/251‑1 (MQ=255)                                                                                                                                                                                              ttCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGCGCAACACCCATTGCCATCGTGTCCTTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGCCAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCAc                                                 >  1:227139/1‑251 (MQ=255)                                                                                                                                                                                                 cGGAATGCCATCTGTACCGACGACCAGCCCTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCgg                                              <  2:330712/251‑1 (MQ=255)                                                                                                                                                                                                          aTCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTAc                                     >  2:235083/1‑251 (MQ=255)                                                                                                                                                                                                            cTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTCATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACCAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCCACGGAGGCGTTACCg                                   >  2:531811/1‑251 (MQ=255)                                                                                                                                                                                                                 ccGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCttt                              <  1:216369/251‑1 (MQ=255)                                                                                                                                                                                                                      gacCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCGTTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCAGTACGACGCGCTTCTCCGGAGGCGTTACCGGCTTTACCCt                         <  2:202956/251‑1 (MQ=255)                                                                                                                                                                                                                        ccAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTCCCCTCa                       <  2:485721/251‑1 (MQ=255)                                                                                                                                                                                                                         cAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGAACGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCaa                      <  1:10478/251‑1 (MQ=255)                                                                                                                                                                                                                           gCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATg                    <  2:445807/251‑1 (MQ=255)                                                                                                                                                                                                                            cACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATg                    >  1:498176/1‑250 (MQ=255)                                                                                                                                                                                                                            cACTTTTCCCGTGTTTCCGGCGGTTGTGCCTGCCTGAGCTACACCCATTGCCAGCGTGTCATTTGCGTAAGAGATCGATTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCGGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGa                   <  1:211052/251‑1 (MQ=255)                                                                                                                                                                                                                             aCTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGAt                  >  2:229596/1‑251 (MQ=255)                                                                                                                                                                                                                                tttCCCGTTTTTCCGGCGTTTGTGACGGCCTGAGCAACACCCATGGCCAGCGTGTCATTCGCGCAATAGATCGTTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCGTGCTACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGTCACCATTACGACGCGCTTCGCCGGAGGCGTTACCGGCTTTCCCCTCAATGAtt                 <  2:417100/249‑1 (MQ=255)                                                                                                                                                                                                                                 ttCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATGCGGATTACGTTGCAACACGTTAGTGGCGACAGCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGc               <  2:102378/250‑1 (MQ=255)                                                                                                                                                                                                                                   cccGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATTATTGCGAc            >  2:445120/1‑251 (MQ=255)                                                                                                                                                                                                                                      gTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTc         <  2:542623/251‑1 (MQ=255)                                                                                                                                                                                                                                           gccGGCGTTTGTGACTGCCTGAGCAACACCCATGGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTGAGTGGCGACGTCCAGTGCTGTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCCTCACCGTAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTcacca     <  2:134383/249‑1 (MQ=255)                                                                                                                                                                                                                                             cGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTCACCACCt  >  1:174000/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  CGGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGCGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTCACCACCT  >  NC_000913/4311092‑4311578

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