Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,116,497	G→A	H872H (CAC→CAT)	yghJ ←	putative lipoprotein YghJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,116,497	0	G	A	100.0%	238.7 / NA	71	H872H (CAC→CAT)	yghJ	putative lipoprotein YghJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/44);  total (27/44)

CGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTCCCCTGTCTTACTATCGATGGTGTACGGCA  >  NC_000913/3116254‑3116741                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     cGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCGCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATg                                                                                                                                                                                                                                               <  2:383240/251‑1 (MQ=255) cGACATGTTCTGGTACAGACGTGCCAGATCGACGCTGCTCAGACGCTCACCTGTGCGACCATTGGTCCGGAAGTAGAGCTGATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTGTTGCAGTGTCGGCGTGAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCAGCAGTAACCGGAGCTCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATg                                                                                                                                                                                                                                               <  1:296748/251‑1 (MQ=255)  gACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGc                                                                                                                                                                                                                                              <  2:496627/251‑1 (MQ=255)   aCATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGcc                                                                                                                                                                                                                                             <  2:141098/251‑1 (MQ=255)       gTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCAGAATGATATGCCGGa                                                                                                                                                                                                                                          <  1:570034/250‑1 (MQ=255)            ggTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATTATATGCCGGATTa                                                                                                                                                                                                                                       <  2:81719/248‑1 (MQ=255)                cAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACAct                                                                                                                                                                                                                                <  2:342501/251‑1 (MQ=255)                        ccAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAg                                                                                                                                                                                                                         <  1:606171/250‑1 (MQ=255)                         cAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGcc                                                                                                                                                                                                                       <  1:351300/251‑1 (MQ=255)                          aGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGccc                                                                                                                                                                                                                      >  1:477753/1‑251 (MQ=255)                            aTCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCg                                                                                                                                                                                                                    <  2:617442/251‑1 (MQ=255)                             tCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGCTACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCAGGATTAGAACACTCTTTCAGCCCCgg                                                                                                                                                                                                                   >  1:356806/1‑251 (MQ=255)                              cGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCgg                                                                                                                                                                                                                   <  1:482034/250‑1 (MQ=255)                              cGACGCTGCTCAGACGCTCAACTTTGCGACCATTGGTCCGGAAGTAAAGCTAATAATGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATAGCTTTTGAAGCGTCGGCGTTAAGGCTTAGTTGCGTATAATGTGGAACAT‑ACATGCCACCAGTAACCGAACCCCCCGTGCCAGGACGATATTCCACGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCggg                                                                                                                                                                                                                  >  1:415270/1‑251 (MQ=255)                                  gctgctCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAAcg                                                                                                                                                                                                              <  1:351821/251‑1 (MQ=255)                                     gctCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACgcgg                                                                                                                                                                                                           <  1:105362/251‑1 (MQ=255)                                        cAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACgcggc                                                                                                                                                                                                          <  1:340378/249‑1 (MQ=255)                                        cAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACgcggc                                                                                                                                                                                                          <  2:651039/249‑1 (MQ=255)                                                cACCTTTTCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTGAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATTATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTtc                                                                                                                                                                                                <  2:385314/251‑1 (MQ=255)                                                    tttGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCCTCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCt                                                                                                                                                                                            >  1:341679/1‑251 (MQ=255)                                                      tGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGAGTTAAGGCTTAGTTGCGTATACTGTGGAACCT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTc                                                                                                                                                                                          >  2:9241/1‑251 (MQ=255)                                                       gCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTc                                                                                                                                                                                          <  1:168326/250‑1 (MQ=255)                                                        cGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCTCCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGc                                                                                                                                                                                        <  2:223501/251‑1 (MQ=255)                                                                ggTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTc                                                                                                                                                                                 <  2:181171/250‑1 (MQ=255)                                                                         gTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCTCCATCGCTTTTGCAGTGTCGGCGTTAAGGCGTAGTTGCGTATACTGTGGAACAT‑ACATTCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCct                                                                                                                                                                       <  2:395030/251‑1 (MQ=255)                                                                           agagCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCctc                                                                                                                                                                      <  1:184708/250‑1 (MQ=255)                                                                                       gATACAGACGCTGAATGGTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAG‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATc                                                                                                                                                          <  2:116677/250‑1 (MQ=255)                                                                                           cAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCt                                                                                                                                                     <  2:206751/251‑1 (MQ=255)                                                                                                    aaTGTTGGTGCCGAAATTGGCAGCCTGCACCAGCTCTTTTTCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACGGTGGAACATTAGTGGCCAC‑AGTAACCGGAGTACACATGTCAGGAAGATATTCCAGGCAGTCGACCTCATAATGATTTGCCGGATTAGTACTCTCTGTCAGCCCCGGGAACTCGGCGAAAGTTTTCTCCTTCTCAGCCTGCAGAGAATCCTCGGTTTTATTAGCGGCGTCATCAATaa                                                                                                                                            <  2:438550/251‑1 (MQ=255)                                                                                                      tGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCTCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATaa                                                                                                                                            >  2:506033/1‑249 (MQ=255)                                                                                                            tGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATa                                                                                                                                     >  1:394872/1‑250 (MQ=255)                                                                                                             gCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAAc                                                                                                                                   <  1:469337/251‑1 (MQ=255)                                                                                                                  aaTCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAATGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACCCTCTTTCAGCCCGTGGAACGCGGCGAAAATTTTCTCCTTCGCATCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGTCATAACGCGtt                                                                                                                              <  2:259555/251‑1 (MQ=255)                                                                                                                    tcagaatccTGCCCCATCGCTTTTGAAGTGTCGGCGTTAAGGCTTAGATATGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCACCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAATATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGAGCGGCCTCACAAATAAAGGCATAACGCGTTTc                                                                                                                            >  2:325459/8‑251 (MQ=255)                                                                                                                       gCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGt                                                                                                                                                                   >  2:186163/1‑209 (MQ=255)                                                                                                                       gCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGt                                                                                                                                                                   <  1:186163/209‑1 (MQ=255)                                                                                                                       gCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGTCCTCATCAATAAAGGCATAACGCGTTTCCTg                                                                                                                         <  1:169430/251‑1 (MQ=255)                                                                                                                       gCAGCCTGCACCATCGCTTTGTCATTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGATCCCACGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCtt                                                                                                                         <  2:336825/251‑2 (MQ=255)                                                                                                                           ccTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTg                                                                                                                     >  1:397521/1‑251 (MQ=255)                                                                                                                              gCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATCTGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGccc                                                                                                                  >  2:155359/1‑250 (MQ=255)                                                                                                                               cACCATCGCTTTTGCAGTGTCGGCGTTAATGCTTAGTTGCGTATACTGTGTAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGTCCTCATCAATAAAGGCATAACGCGTTTCCTTTTTGCCAt                                                                                                                 <  2:311385/251‑1 (MQ=255)                                                                                                                                aCCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATc                                                                                                                >  1:556157/1‑251 (MQ=255)                                                                                                                                   aTCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTAc                                                                                                             <  2:200348/251‑1 (MQ=255)                                                                                                                                          ttGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTc                                                                                                       <  1:465257/250‑1 (MQ=255)                                                                                                                                                   ccgCGTTAAGGCCTAGTTGCGTTTACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCCGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCTCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCCGCCTCATCAATAAAGGCATATCTCGTTTCCTGTTTTCCATCTACATGTTCCAGCCAGCtt                                                                                             <  2:544686/249‑2 (MQ=255)                                                                                                                                                   cGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTg                                                                                             >  2:631143/1‑251 (MQ=255)                                                                                                                                                     gCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCGCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGc                                                                                           <  1:645993/251‑1 (MQ=255)                                                                                                                                                       gTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCa                                                                                          >  1:609258/1‑250 (MQ=255)                                                                                                                                                       gTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCAGAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCa                                                                                          >  1:573831/1‑250 (MQ=255)                                                                                                                                                            ggCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTcc                                                                                    <  1:571271/251‑1 (MQ=255)                                                                                                                                                                      cGTATACTGTGGCACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCCGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGtt                                                                          >  2:67230/1‑251 (MQ=255)                                                                                                                                                                       gTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGttt                                                                         >  1:344905/1‑251 (MQ=255)                                                                                                                                                                         ataCTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTg                                                                        <  1:596288/250‑1 (MQ=255)                                                                                                                                                                                  aaCAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCATGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGgtt                                                              <  2:337452/251‑1 (MQ=255)                                                                                                                                                                                   acat‑acatTCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCATGCAGTGGACCTCATAATGATATGCCGGATTAGTACACTCTTTCACCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCTTTCAGAGAATCCTCTGTTTTATGATCTGCCTCATTAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCCGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGgttt                                                             <  1:366078/251‑1 (MQ=255)                                                                                                                                                                                   acat‑acatGCCTCCATTAACCGGACCCTCCGTGCCAGGACGATATTCCATGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCGTTCATCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTGCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCGTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGgttt                                                             <  2:365830/251‑1 (MQ=255)                                                                                                                                                                                    cat‑acatGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGgtttg                                                            >  1:389970/1‑251 (MQ=255)                                                                                                                                                                                         catGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAATACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGtttt                                                        <  2:70985/251‑1 (MQ=255)                                                                                                                                                                                            gCCACCAGTAACCGGATCCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCt                                                      <  2:460926/250‑1 (MQ=255)                                                                                                                                                                                               accaGTAACGCGAAACCACGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGCTATACCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTCTTCTACCt                                                  >  2:294959/1‑251 (MQ=255)                                                                                                                                                                                                          gAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTGTTCTACCTGATATTTCCAc                                       <  1:510688/251‑1 (MQ=255)                                                                                                                                                                                                              ccccGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTaa                                   <  1:329600/251‑1 (MQ=255)                                                                                                                                                                                                                  gTGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTc                               >  2:210737/1‑251 (MQ=255)                                                                                                                                                                                                                   tGCCAGGACGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTcc                              >  2:197389/1‑251 (MQ=255)                                                                                                                                                                                                                          aCGATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTAGACTCTTTCAGCCCATGGAACGCGGCGAAACTTTTCTCCTTCGCAGCATTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATCACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCACCTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTATTCCTGCTATTTTCACTTCACTTCCCCTGTCt                       >  2:241808/1‑251 (MQ=255)                                                                                                                                                                                                                            gATATTCCAGGCAGTTGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTCCCCTGTCtt                      >  1:47560/1‑250 (MQ=255)                                                                                                                                                                                                                                   cAGGCAGTTGACCTCATAATGATATGCCGGATTAGCACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTCCCCTGTCTTACTATCg               >  2:657090/1‑250 (MQ=255)                                                                                                                                                                                                                                    aGGCAGTTGACCTCATAATGATCTGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCCGCGAAAATTTTCTCCTTCGCAGCCTTCAAAGCATCCTCTGTTTTATGCTCGGCCTCATCAATAAAGGCCTAACGCGTTTCCTGTTTGCCATCTACATCTTACAGCCAGCTGTCAACTTCCAGATTCGGTTCGTCATCAGTCTTGTTTTCTACCTGATATTTCCACTTCAATTCCCCTGTCTTCCTATCGAt             >  2:199341/1‑251 (MQ=255)                                                                                                                                                                                                                                           tGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTCCCCTGTCTTACTATCGATGGTGTAc      >  1:30779/1‑251 (MQ=255)                                                                                                                                                                                                                                           tGACCTCATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTCCCCTGTCTTACTATCGATGGTGTAc      >  1:545507/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cATAATGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTCCCCTGTCTTACTATCGATGGTGTACGGCa  <  2:389970/249‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CGACATGTTCTGGTACAGACGTTCCAGATCGACGCTGCTCAGACGCTCACCTTTGCGACCATTGGTCCGGAAGTAGAGCTCATGCTGATACAGACGCTGAATGTTGGTGCCTAAATCCGCAGCCTGCACCATCGCTTTTGCAGTGTCGGCGTTAAGGCTTAGTTGCGTATACTGTGGAACAT‑ACATGCCACCAGTAACCGGAACCCCCGTGCCAGGACGATATTCCAGGCAGTTGACCTCATAGTGATATGCCGGATTAGTACACTCTTTCAGCCCCGGGAACGCGGCGAAAATTTTCTCCTTCGCAGCCTTCAGAGAATCCTCTGTTTTATGATCGGCCTCATCAATAAAGGCATAACGCGTTTCCTGTTTGCCATCTACATCTTCCAGCCAGCTGGCAACTTCCAGCTTCGGTTTGTCATCAGGTTTGTTTTCTACCTGATATTTCCACTTAACTTCCCCTGTCTTACTATCGATGGTGTACGGCA  >  NC_000913/3116254‑3116741

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