Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,272,711	C→T	E387E (GAG→GAA)	uvrA ←	UvrABC excision nuclease subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,272,711	0	C	T	100.0%	199.4 / NA	63	E387E (GAG→GAA)	uvrA	UvrABC excision nuclease subunit A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/37);  total (26/37)

CAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGCATCTGGAAATAATAGAAGTTGCGGCGATCCCAG  >  NC_000913/4272471‑4272945                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           cAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAgc                                                                                                                                                                                                                                  >  2:85377/1‑251 (MQ=255)    aCGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAgcggc                                                                                                                                                                                                                               >  2:17053/1‑251 (MQ=255)             gATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGGGATAGCAGGCAGCGGCGTATTCTCGACATACCCGTGGCGCGCTCCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGAGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCgtat                                                                                                                                                                                                                      >  2:393692/1‑251 (MQ=255)              aTCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTACATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGCCGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTACTTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCTtatt                                                                                                                                                                                                                     >  2:427727/1‑251 (MQ=255)                 tctTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATg                                                                                                                                                                                                                  >  2:413420/1‑251 (MQ=255)                    ttAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATgca                                                                                                                                                                                                                <  1:389481/250‑1 (MQ=255)                     tAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATgcag                                                                                                                                                                                                               <  1:152864/250‑1 (MQ=255)                        agttttttcCTGCAATCTTCTCCCGCTGACCTGCTAGGTTGATATTGTTGAAGAAATCCATCGCATGACCAATGCTCAGTTCGGAGATTGCAGGCAGCTGTGTATTTCCGACATACACGTGGCGCGCTTCCCGACGCTGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTGCGTTATAGCGGCGCTCCATATTATGCAGCACg                                                                                                                                                                                                           <  2:268472/248‑1 (MQ=255)                              tttCTGCAATCTTCGCCCGTTGACCTGCGATTGTGAGATTGTTTAAGAATTCCATCGACTGTCCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTc                                                                                                                                                                                                      <  1:261189/250‑1 (MQ=255)                                   gCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCTTGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAAc                                                                                                                                                                                                  <  1:355771/249‑1 (MQ=255)                                          tCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATgac                                                                                                                                                                                           <  1:212242/249‑1 (MQ=255)                                               cGCTGACCGGCGATGTTGAGATTGTTGAAGAATTCCATCGAATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGaac                                                                                                                                                                                    <  2:99070/251‑2 (MQ=255)                                                gCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATg                                                                                                                                                                                   <  1:44904/251‑1 (MQ=255)                                                 cTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGAAGTTAGCAGGCAGCGGCGAATTCACGACATAAACGTGGCGAGCTTCACGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAACCTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTCTAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGCCGACTgg                                                                                                                                                                                  >  1:465275/1‑251 (MQ=255)                                                     ccTGCGAGTGTGAGATTGTTGAAGAAGTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGCCATACACGTGGCGCGCTTCCCGACGCAGTCGCGTCCCTTCGCAGCTGTCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATGCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCTTTCGAACGGATGACGACGAATGGAGGt                                                                                                                                                                              <  1:199580/251‑1 (MQ=255)                                                       ggCGAGTTTGAGAGTGTTGAAGAATGCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGGATTCGCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACGTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCGTTATAGCGGCGCGCCATAGTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAgggat                                                                                                                                                                            <  2:223480/250‑4 (MQ=255)                                                            gTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGCATGGAGGTATCGccc                                                                                                                                                                       >  1:42484/1‑250 (MQ=255)                                                                         gaagaaTTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATg                                                                                                                                                           <  2:386743/250‑1 (MQ=255)                                                                         gaagaaTTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATg                                                                                                                                                           <  1:345568/250‑1 (MQ=255)                                                                         gaagaaTTCAATCGCATGACCAAGGCGCATGGCGGAGATAGCAGGCAGCGGCGGATGCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTATCGCTCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGt                                                                                                                                                          <  2:143098/251‑1 (MQ=255)                                                                          aagaaTTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGt                                                                                                                                                          >  1:153955/1‑250 (MQ=255)                                                                          aagaaTTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGt                                                                                                                                                          >  1:429688/1‑250 (MQ=255)                                                                           agaaTTCCATCGCAGGACCAATGCTCTTGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTAt                                                                                                                                                        <  2:186337/251‑1 (MQ=255)                                                                             aaTTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTAttt                                                                                                                                                      >  1:217628/1‑251 (MQ=255)                                                                              aTTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTg                                                                                                                                                     <  1:320395/251‑1 (MQ=255)                                                                                     cGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGTCGAATGGAGGTATCGCCACAATCGTTCATGTCTTTGAATTCa                                                                                                                                               >  2:195184/1‑250 (MQ=255)                                                                                      gCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAAt                                                                                                                                             >  1:198477/1‑251 (MQ=255)                                                                                        aTGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGt                                                                                                                                           <  1:387341/251‑1 (MQ=255)                                                                                           aCCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGtttt                                                                                                                                        >  1:267001/1‑251 (MQ=255)                                                                                                     aTGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGa                                                                                                                              <  1:17053/251‑1 (MQ=255)                                                                                                         cGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAAcc                                                                                                                           <  1:24194/250‑1 (MQ=255)                                                                                                         cGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCg                                                                                                                          >  1:75924/1‑251 (MQ=255)                                                                                                         aggAGATAGCAGGCAGCGGCGTATTCTCGAAATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTAATGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGAGGCGCTACATATTATGAAGCACGCCTTCGAACGGATGACGAGGAATGGAGGTAGCGCCACGATAGTTCATGAATTTGAATTCAATGTGTTCTTTGCCAGAACCg                                                                                                                          >  1:75925/2‑251 (MQ=255)                                                                                                             gaTAGCAGGCAGCCGCGTATTCCCGTCATACCCCGGTCGCGCGTCCCGACGCAGACGCGGCCCTTCGCAGCGTGCGCACTGACGATTACTGATAAACTTGGCTAAGTCTTCGCGTACCGCGCTGGATTCCGGTTCTTTATAGCGGCGCTTCATATTATGCAGCACGCCTTCGAACGGATGACGACGTATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTaca                                                                                                                      <  2:405479/251‑1 (MQ=255)                                                                                                                aGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCATACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTacaaca                                                                                                                   <  1:421028/251‑1 (MQ=255)                                                                                                                aGCAGGCAGAGGCGTAATCTCGACATACACGTGGCGCGCTCCCAGACGCAGACGCGTCCCTTCGCCGATGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCGGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGGAGAACGCCTTCGAACGGATGACGACGAATGGAGGTATCGACCCGATCGTTCCTGTATTTGAATTCAAGGTTTTCTTTGACAGCACCgaaaaac                                                                                                                    >  2:248802/1‑244 (MQ=255)                                                                                                                   atgcagCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAAcacca                                                                                                                <  1:320674/249‑1 (MQ=255)                                                                                                                     gaagCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGAAGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACGGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCGTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACt                                                                                                              <  2:246152/249‑1 (MQ=255)                                                                                                                        gcggcgTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCCCGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGcccg                                                                                                                               >  2:256508/1‑229 (MQ=255)                                                                                                                        gcggcgTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAg                                                                                                                               <  1:256508/231‑1 (MQ=255)                                                                                                                         cggcgTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTAt                                                                                                          <  1:380862/251‑1 (MQ=255)                                                                                                                          ggggTATTCTCGACTCACACGTGGCGCGCTTCCCGCCGCAGACGCGACACTTCGCAGCTGGCGCACGGACGATTACTGTTAAACTAGGCTAATTCTTCGCGTACCGCGCTGGCTTCCGTTTCTTTACAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGGCGAATGTAGGTATCGCCACGATCGTTCATTTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTCCAAACCCACtttctg                                                                                                         >  2:17605/4‑248 (MQ=255)                                                                                                                             gTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCAc                                                                                                      >  1:278967/1‑251 (MQ=255)                                                                                                                               attatcGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCTCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCCCGATCGTTCATGTATGTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGt                                                                                                    <  2:42484/247‑1 (MQ=255)                                                                                                                                           cacGTGGCGCGCGTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGc                                                                                        <  1:385550/251‑1 (MQ=255)                                                                                                                                             cctGGCGCTCTTCCCGACTCAGACGCGTCCCTTCGCAGCTGGAGCACGGACGATTACTGATAAACTGGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTg                                                                                      <  2:361254/249‑1 (MQ=255)                                                                                                                                                 gcgcgcTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGcccc                                                                                  <  1:149515/251‑1 (MQ=255)                                                                                                                                                  cgcgcTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGcccc                                                                                  >  1:335022/1‑250 (MQ=255)                                                                                                                                                              cGCAGACGCGGCCCTTCGCAGCGTGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCGTTATAGCGGCGCTCCATATTAAGCAGCACGCCGTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTGCGCGCTCAGGCTGCACCACGGCGCTTCGAc                                                                     <  2:64124/251‑1 (MQ=255)                                                                                                                                                                   aCGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGa                                                                <  1:85377/251‑1 (MQ=255)                                                                                                                                                                   aCGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGa                                                                <  1:493659/251‑1 (MQ=255)                                                                                                                                                                    cgcgTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGAGCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACGTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGaa                                                               <  2:278967/251‑1 (MQ=255)                                                                                                                                                                        tCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGCTGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACttt                                                           >  1:185308/1‑250 (MQ=255)                                                                                                                                                                          ccTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTtat                                                          >  1:449003/1‑250 (MQ=255)                                                                                                                                                                                            ggACGCTTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTCTTTGCATTCAATGTTTTCTTTGCCAGAACCGTACCACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAAATTATAGTGATCTGCCAGCGAttt                                       >  2:181980/1‑251 (MQ=255)                                                                                                                                                                                                gATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGc                                   <  1:398478/251‑1 (MQ=255)                                                                                                                                                                                                   tACTGATAAACTTGGCTAATTCTTCGCGTCCCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGCAt                                 <  2:217628/250‑1 (MQ=255)                                                                                                                                                                                                            aCTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGCATCTGGAaata                        >  1:458215/1‑250 (MQ=255)                                                                                                                                                                                                                 gCTAATTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGCATCTGGAAATAATAGaa                  <  1:102113/251‑1 (MQ=255)                                                                                                                                                                                                                    aaTTCTTCGCGTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGCATCTGGAAATAATAGAAg                 <  2:444898/249‑1 (MQ=255)                                                                                                                                                                                                                              gTACCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGCATCTGGAAATAATAGAAGTTGCGGCGATcc     >  1:270811/1‑251 (MQ=255)                                                                                                                                                                                                                                aCCGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGTACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGCATCTGGAAATAATAGAAGTTGCGGCGATccc    <  2:153955/250‑1 (MQ=255)                                                                                                                                                                                                                                 ccGCGCTGGATTCCGTTTCTTTATAGCGGCGCTCCATATTATGCAGCCCGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACCACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGAGATCTGCCAGCGCTTTCCGCATCTGGACATAATAGAAGTTGCGCCGATCCCac  >  2:503037/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           CAGACGATCGCCGATCTCTTTAAGGATTTTTTCTGCAATCTTCGCCCGCTGACCTGCGAGTTTGAGATTGTTGAAGAATTCCATCGCATGACCAATGCTCATGTCGGAGATAGCAGGCAGCGGCGTATTCTCGACATACACGTGGCGCGCTTCCCGACGCAGACGCGTCCCTTCGCAGCTGGCGCACGGACGATTACTGATAAACTTGGCTAATTCTTCGCGTACCGCGCTGGATTCCGTCTCTTTATAGCGGCGCTCCATATTATGCAGCACGCCTTCGAACGGATGACGACGAATGGAGGTATCGCCACGATCGTTCATGTATTTGAATTCAATGTTTTCTTTGCCAGAACCGTACAACACCACTTTATGCACGTTCGCGCTCAGGCTGCCCCACGGCGCTTCGACGTCGAACTTATAGTGATCTGCCAGCGATTTCAGCATCTGGAAATAATAGAAGTTGCGGCGATCCCAG  >  NC_000913/4272471‑4272945

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