Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,153,836	G→A	Q338* (CAA→TAA)	argE ←	acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,153,836	0	G	A	100.0%	208.7 / NA	55	Q338* (CAA→TAA)	argE	acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/28);  total (27/28)

GATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTTCATTGAGTGTCATGCCAGGCAGCGGACGAATATCCATATGCA  >  NC_000913/4153601‑4154072                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               gATAAGAGGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTc                                                                                                                                                                                                                                 <  2:20509/251‑1 (MQ=255)  aTAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCg                                                                                                                                                                                                                                >  1:210078/1‑251 (MQ=255)     aGATGCGCAAGCATCGCATCTGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCTGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCTTGTTTCCAGATATTCAGCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCTAACTCCATCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTa                                                                                                                                                                                                                             <  1:314411/251‑1 (MQ=38)      gATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTAc                                                                                                                                                                                                                            >  2:133156/1‑251 (MQ=255)            cAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACagtagt                                                                                                                                                                                                                      >  1:97567/1‑251 (MQ=255)               gcatcgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTca                                                                                                                                                                                                                   <  1:133156/251‑1 (MQ=255)                   cgcatcCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTcaccac                                                                                                                                                                                                               >  1:217274/1‑251 (MQ=255)                           gcATTATTTCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTTCCTGGGTTATCAGTTCGCTGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGttt                                                                                                                                                                                                       <  2:451725/250‑1 (MQ=255)                            cagtattGCGGCACTGGAGTTGGGCAACAGCGGCGGATGCGGCGCGAGCGGCGTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTTGGCTTGATAAACCGGGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCGCTTCGGtttt                                                                                                                                                                                                      <  2:333857/248‑1 (MQ=255)                                   gCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTcc                                                                                                                                                                                                <  2:103109/250‑1 (MQ=255)                                    cGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGGTATCAGTTGGCGGGTGGGCTTGATAAACCGTGTTTCCAGATAGCCATCAGGTTGATGAGCCTGATTAATTGAGCGAGGCCCCAGCACCAGCGGCGGGCATAACGGTTAAAGAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGa                                                                                                                                                                                              <  1:73218/251‑1 (MQ=255)                                     ggCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAg                                                                                                                                                                                             <  2:120072/251‑1 (MQ=255)                                        aCTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAACTGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATACACGGCGCTTCGGTACAGTAGTTCACCCCTTCGGTTTTTGCTCCGAGCa                                                                                                                                                                                           >  2:272626/1‑250 (MQ=255)                                          tGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAAtt                                                                                                                                                                                        >  1:96905/1‑251 (MQ=255)                                                     aaCAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAAcccct                                                                                                                                                                             >  2:442753/1‑248 (MQ=255)                                                              ggATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGt                                                                                                                                                                    >  2:81828/1‑251 (MQ=255)                                                                           gcgcCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTg                                                                                                                                                       >  1:212388/1‑251 (MQ=255)                                                                                 tATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCa                                                                                                                                                   <  2:129449/249‑1 (MQ=255)                                                                                  aTCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTc                                                                                                                                                >  1:470959/1‑251 (MQ=255)                                                                                    ccGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCAt                                                                                                                                              <  2:210381/251‑1 (MQ=255)                                                                                      ggtctaagTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTGCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGAGGATTCGGTGGGCATTCATAg                                                                                                                                            <  2:422902/244‑1 (MQ=255)                                                                                       gCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGc                                                                                                                                           <  2:192941/251‑1 (MQ=255)                                                                                            cGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAggg                                                                                                                                      <  2:215160/251‑1 (MQ=255)                                                                                             gTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAggg                                                                                                                                      >  2:233431/1‑250 (MQ=255)                                                                                                  aaTGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATcgg                                                                                                                                 <  1:272626/250‑1 (MQ=255)                                                                                                        aGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGc                                                                                                                          <  1:541091/251‑1 (MQ=255)                                                                                                                   tGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACTAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGAc                                                                                                               >  1:426676/1‑251 (MQ=255)                                                                                                                      aTTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGaatccgctcgttcatcgt                                                                                                            >  2:409805/1‑234 (MQ=255)                                                                                                                                ttATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGAGGCAGCACGTCGCCCCTCTGACGAccc                                                                                                  >  2:106542/1‑251 (MQ=255)                                                                                                                                    cAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGcc                                                                                              <  2:212199/251‑1 (MQ=255)                                                                                                                                          gcgGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGtt                                                                                       <  2:442496/251‑1 (MQ=255)                                                                                                                                                      gATAAACCGTTTTTCCAGATATTCATCAGGTTTATGAGCCTGATTAATTGAGCCAGGCCCCAACTCCATCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCCCCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTGGATGATTCGGTGGGCATTCATAGCCAGGGATCGTCGGATGCAGCGCGTCGACCGTCAGACTACCCGTCC‑TGCGTTCGCTCACCGGAg                                                                           <  2:511582/251‑1 (MQ=255)                                                                                                                                                           aCCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGCGGGCCTTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGCCATGCGGACCGGCC‑AGCGGTCGCGCCCCGGCGCCaa                                                                       >  2:516317/1‑250 (MQ=255)                                                                                                                                                            ccGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGGGCCAATg                                                                     >  1:143445/1‑251 (MQ=255)                                                                                                                                                                gtTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCAt                                                                  <  2:399338/250‑1 (MQ=255)                                                                                                                                                                    ccAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGtt                                                              <  1:81828/250‑1 (MQ=255)                                                                                                                                                                            tcatcaGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAAc                                                     <  1:225926/250‑1 (MQ=255)                                                                                                                                                                              atcaGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGGCAGACGACCCGGCC‑AGCGGTCGCCCCCCGGGGGCAATTCATCGGTGTGC‑AAACCa                                                   >  2:77902/1‑250 (MQ=255)                                                                                                                                                                              atcaGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCCGCCGACGACCCGGAC‑CGCGTTCGCTCCCCCGGGCCCCTGCAACGTTGAGCCAAAcct                                                   >  2:147611/1‑250 (MQ=255)                                                                                                                                                                                caGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTa                                                <  2:58396/251‑1 (MQ=255)                                                                                                                                                                                   gTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGGTCGCTCACCGGAGCCAAAGCATCGTTGAGC‑AAAACATTaaa                                              >  1:375451/1‑249 (MQ=255)                                                                                                                                                                                       atgaGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTTCAt                                         <  1:233431/251‑1 (MQ=255)                                                                                                                                                                                           gCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGGCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTTCATTg                                       >  1:155071/1‑249 (MQ=255)                                                                                                                                                                                             cTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGGTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGATCATTGAgtg                                   >  1:410697/1‑251 (MQ=255)                                                                                                                                                                                               gATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTGTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTTCATTGAgtg                                   <  1:98354/249‑1 (MQ=255)                                                                                                                                                                                                aTTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑ACCGTTCGCTCACCGGAGCCAATGCATCGTTGAAC‑AAACCATTAAGTTCATTGAGGGTc                                 >  1:213880/1‑250 (MQ=255)                                                                                                                                                                                                        aGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACGTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTGCATTGAGTGTCATGCCAGGc                        <  2:453021/251‑1 (MQ=255)                                                                                                                                                                                                         gCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTGCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATGCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTTCATTGAGTGTCATGCCAGGCa                       <  1:241102/251‑1 (MQ=255)                                                                                                                                                                                                           cAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTTCATTGAGTGTCATGCCAGGCAGc                     <  1:116444/251‑1 (MQ=255)                                                                                                                                                                                                           cAGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAAATCATCGTTGAGC‑AAACCAATAAGTTCATTTAGGGTCATGCCAGGCAGc                     >  1:459085/1‑251 (MQ=255)                                                                                                                                                                                                            aGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCCCTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGGCCGTCCGACGACCCGGCC‑AGCGGTCGCTCCCCGGAGCCCATGCCTCGTTGAGA‑AAACCATTAAGGTCATTTAGGGTCATGCCAGGCAg                      >  2:441293/1‑249 (MQ=255)                                                                                                                                                                                                            aGGCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTTCATTGAGTGTCATGCCAGGCAGCg                    <  2:308320/251‑1 (MQ=255)                                                                                                                                                                                                             ggCCCCAACACCAGCGTCGGGCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCAATTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCCGCTCGTCGACCGTCAGACGACCCGGAC‑AACGGTCGCTCACCGGAGCCAATGCCTCGTGGAGC‑AAACCATTAAGTTCATTGAGTGTCATGCCCGGCAGCgg                   >  1:465983/1‑251 (MQ=255)                                                                                                                                                                                                                       ccAGCGTCGGGCATAACGTTTAAATAAAAGGAGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCCTCGGCCGTCAGACGACCCGGCCCAG‑GTTCGCTCACCGGGGCCAAATCATCGTTGAGC‑AAACCCTTAAATTCATTTAGTTGCATACCAGGCCGCGGACGGATATcc         >  2:333956/1‑251 (MQ=255)                                                                                                                                                                                                                               gggCATAACGTTTAAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCCGCCCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCCCCCGAGCCAATGCATCGTTGAGC‑CAACCATTAAGTTCATTTAGTGTCCTGCCACGCCGCGGACGAAAATCCATATGCa  >  2:18799/1‑250 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               GATAAGATGCGCAAGCATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCTCCGAGCAATTTCTCAACCACTTCAACCAGTTGATGATTCGGTGGGCATTCATAGCCAGGGATCGGCGGATGCAGCTCGTCGACCGTCAGACGACCCGGCC‑AGCGTTCGCTCACCGGAGCCAATGCATCGTTGAGC‑AAACCATTAAGTTCATTGAGTGTCATGCCAGGCAGCGGACGAATATCCATATGCA  >  NC_000913/4153601‑4154072

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