Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,215,122	T→G	G281G (GGT→GGG)	metA →	homoserine O‑transsuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,122	0	T	G	100.0%	317.5 / NA	88	G281G (GGT→GGG)	metA	homoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (38/50);  total (38/50)

CACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTAT  >  NC_000913/4214894‑4215337                                                                                                                                                                                                                                     |                                                                                                                                                                                                                        cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTAc                                                                                                                                                                                                                  <  2:436729/236‑1 (MQ=255)  acCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACt                                                                                                                                                                                                                 <  1:506673/236‑1 (MQ=255)     gATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGt                                                                                                                                                                                                               >  1:319004/1‑235 (MQ=255)      aTCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGttt                                                                                                                                                                                                             >  2:353406/1‑236 (MQ=255)                tCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTgg                                                                                                                                                                                                    <  2:86150/235‑1 (MQ=255)                tCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGc                                                                                                                                                                                                   >  1:625258/1‑236 (MQ=255)                  tGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTc                                                                                                                                                                                                 <  2:625258/236‑1 (MQ=255)                  tGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGAATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTc                                                                                                                                                                                                 <  2:433197/236‑1 (MQ=255)                   ggCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCa                                                                                                                                                                                                >  1:486588/1‑236 (MQ=255)                   ggCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATAAACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCa                                                                                                                                                                                                >  2:708775/1‑236 (MQ=255)                            ggAAGGAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTa                                                                                                                                                                                            <  1:564710/231‑1 (MQ=255)                              aagaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGt                                                                                                                                                                                      >  2:69284/1‑235 (MQ=255)                                  agGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTAtt                                                                                                                                                                                          <  1:708775/227‑1 (MQ=255)                                  agGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTAtt                                                                                                                                                                                          <  1:546965/227‑1 (MQ=255)                                  agGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGACTCAACTATTACGTCTAc                                                                                                                                                                                  <  2:319004/235‑1 (MQ=255)                                    gggATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGCGCGAGCTGGCGTCGTCACGGGAATTTACTGCTTACCAACTGGCTCAACTATTACGTCTTccc                                                                                                                                                                                >  1:365518/1‑234 (MQ=255)                                    gggATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACAACACGAGCGAGATGGCGTAGTCACGGGAATTTACTGTTTCCCAACTGGCCCACATATTACGACTACCa                                                                                                                                                                                >  2:366394/1‑235 (MQ=255)                                     ggATGCATAGCGGTTTGCAAGTAAAGATAAGCCCATTGCCTTTGTGACGTGCAATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGa                                                                                                                                                                              <  2:282736/236‑1 (MQ=255)                                        tGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCTCGAGCGAGCTGGCGTAGTCACGGGAATTTa                                                                                                                                                                                                                   >  2:470158/1‑196 (MQ=255)                                        tGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTa                                                                                                                                                                                                                   <  1:470158/196‑1 (MQ=255)                                        tGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATc                                                                                                                                                                            <  2:142741/235‑1 (MQ=255)                                         gctaatCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGGGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCa                                                                                                                                                                           <  2:329935/231‑1 (MQ=255)                                           atatCTGATTGCCAGTAAAGATAAGCGCATTGCCTTTGGGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACg                                                                                                                                                                         <  2:554171/235‑1 (MQ=255)                                               cTGTTTGCCAGTAAAGATAAGCGCATTGCCTTGGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATa                                                                                                                                                                    <  2:30207/236‑1 (MQ=255)                                                tGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATAc                                                                                                                                                                   >  2:163610/1‑236 (MQ=255)                                                    tGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATc                                                                                                                                                               >  1:592535/1‑236 (MQ=255)                                                    tGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATc                                                                                                                                                               >  2:124654/1‑236 (MQ=255)                                                            aaGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCTGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGGCGTCACGGGAATTTACAGTTTACGAACTGGGTCAACAATTACGTCTACCAGAAGAAGCCTTCCGATCTagtggca                                                                                                                                                       >  2:110875/1‑230 (MQ=255)                                                              gATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACgg                                                                                                                                                          <  2:132101/231‑1 (MQ=255)                                                              gATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACgg                                                                                                                                                          <  1:353406/231‑1 (MQ=255)                                                              gATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACAt                                                                                                                                                     <  2:757503/236‑1 (MQ=255)                                                                    cgcATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATcc                                                                                                                                               >  1:2061/1‑236 (MQ=255)                                                                        ttGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATcc                                                                                                                                               >  1:687713/1‑232 (MQ=255)                                                                        ttGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCaa                                                                                                                                             <  1:757139/234‑1 (MQ=255)                                                                              ttGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTc                                                                                                                                                                                                 >  1:496847/1‑176 (MQ=255)                                                                              ttGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTc                                                                                                                                                                                                 <  2:496847/176‑1 (MQ=255)                                                                               tgtgACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                                                                                                                    >  1:314677/1‑236 (MQ=255)                                                                                gtgACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                                                                                                                   <  1:69284/236‑1 (MQ=255)                                                                                 tgACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGAACTACGGCACATGAATCCAACGATGGAGTaa                                                                                                                                  >  2:535102/1‑236 (MQ=255)                                                                                 taaCGGGCCATCCGGAATATGATGCGCAAACGCTGGCGCAGGAATTTTGCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTaa                                                                                                                                  <  2:722177/234‑1 (MQ=255)                                                                                     gggCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCt                                                                                                                                                                                                  >  1:348961/1‑168 (MQ=255)                                                                                     gggCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCt                                                                                                                                                                                                  <  2:348961/168‑1 (MQ=255)                                                                                           tCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtct                                                                                                                               <  2:2061/229‑1 (MQ=255)                                                                                           tCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCCCAACTATTACGGCTACCCGATACCGCCATCAGATCTACGGCACATGAATCCAACGCTGGATTAAACTCATGTGa                                                                                                                        >  2:632912/1‑236 (MQ=255)                                                                                            cccGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                                                                                        <  1:124654/235‑1 (MQ=255)                                                                                                aaTATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                                                                                             >  1:764107/1‑226 (MQ=255)                                                                                                aaTATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                                                                                             <  2:764107/226‑1 (MQ=255)                                                                                                       tgCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                                                                                             >  1:243236/1‑235 (MQ=255)                                                                                                          gcAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACCCCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                                                                                         <  1:366394/236‑1 (MQ=255)                                                                                                          gcAAACGCTGGCGCAGGAATTTTGCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                                                                                         <  1:437714/236‑1 (MQ=255)                                                                                                                  tGGCGCCGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                                                                                                 <  2:486588/236‑1 (MQ=255)                                                                                                                  tGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                                                                                                 <  1:208453/236‑1 (MQ=255)                                                                                                                      gcAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGACCGTTAAGCGATTCAGCAcc                                                                                             <  2:343176/236‑1 (MQ=255)                                                                                                                                    gcgATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATATTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc                                                                               >  1:292051/1‑236 (MQ=255)                                                                                                                                                    aCTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtttttt                                                               >  2:391859/1‑236 (MQ=255)                                                                                                                                                             ggATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCg                                                        <  1:163610/234‑1 (MQ=255)                                                                                                                                                              gATGTTCCGGATAACTATTTCCCGCACAATGATCCGCAAAATACCCCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGaaa                                                     <  1:632912/236‑1 (MQ=255)                                                                                                                                                                    ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTAc                                                <  1:27486/235‑1 (MQ=255)                                                                                                                                                                       tataACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTca                                            >  1:72532/1‑236 (MQ=255)                                                                                                                                                                            cTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAgg                                       <  2:314677/236‑1 (MQ=255)                                                                                                                                                                                 tCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                                                                                              >  2:376636/1‑160 (MQ=255)                                                                                                                                                                                 tCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                                                                                              <  1:376636/160‑1 (MQ=255)                                                                                                                                                                                 tCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg                                                                          <  1:720864/196‑1 (MQ=255)                                                                                                                                                                                 tCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg                                                                          >  2:720864/1‑196 (MQ=255)                                                                                                                                                                                 tCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTGACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAAt                                  <  2:592535/236‑1 (MQ=255)                                                                                                                                                                                     gCAAAATGATCCGGAAAATTCCCCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAAt                              <  1:645588/236‑1 (MQ=255)                                                                                                                                                                                      cacaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGaaa                                                     >  1:238353/1‑212 (MQ=255)                                                                                                                                                                                      cacaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGaaa                                                     <  2:238353/212‑1 (MQ=255)                                                                                                                                                                                              tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGaaa                                                     >  2:309570/1‑204 (MQ=255)                                                                                                                                                                                              tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGaaa                                                     <  1:309570/204‑1 (MQ=255)                                                                                                                                                                                                cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                                                                                         <  2:152675/150‑1 (MQ=255)                                                                                                                                                                                                cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                                                                                         >  1:152675/1‑150 (MQ=255)                                                                                                                                                                                                 gCAAAATCCACCGCAGGCAAGCTGGCGTAGTCACGGGAATTTCCTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCatat                   <  1:110875/235‑1 (MQ=255)                                                                                                                                                                                                 gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCata                    <  2:687713/234‑1 (MQ=255)                                                                                                                                                                                                          acCGGGAGTGAGCTTGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaac         <  2:292051/236‑1 (MQ=255)                                                                                                                                                                                                          acCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                                                                                              >  1:490932/1‑119 (MQ=255)                                                                                                                                                                                                          acCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                                                                                              <  2:490932/119‑1 (MQ=255)                                                                                                                                                                                                          acCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTg                 <  1:93276/228‑1 (MQ=255)                                                                                                                                                                                                          acCGCGAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTg                 >  2:93276/1‑228 (MQ=255)                                                                                                                                                                                                            cgcgAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCt                                                                                                                                                                                                  >  1:300691/1‑49 (MQ=255)                                                                                                                                                                                                            cgcgAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCt                                                                                                                                                                                                  <  2:300691/49‑1 (MQ=255)                                                                                                                                                                                                            cgcgAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                                                                                              <  2:517689/117‑1 (MQ=255)                                                                                                                                                                                                            cgcgAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                                                                                              >  1:517689/1‑117 (MQ=255)                                                                                                                                                                                                            cgcgAGCGAGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTACcaacaa       >  1:484508/1‑236 (MQ=255)                                                                                                                                                                                                             gcgagcgaGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaac         <  1:661882/233‑1 (MQ=255)                                                                                                                                                                                                             gcgagcgaGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAg      >  1:273021/1‑236 (MQ=255)                                                                                                                                                                                                             gcgagcgaGCTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCATCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATcaac         >  2:661882/1‑233 (MQ=255)                                                                                                                                                                                                                 gcgagcTGGCGTAGTCACGGGAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTAt  >  1:56272/1‑236 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                        CACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTTTTTTATTTCCGAAACGTACCTCAGCAGGTGAATAAATTTTATTCATATTGTTATCAACAAGTTAT  >  NC_000913/4214894‑4215337

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