Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,813	Δ2 bp	coding (142‑143/1161 nt)	metB →	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,811	0	C	.	100.0%	175.3 / NA	38	coding (140/1161 nt)	metB	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (17/21);  total (17/21)
*	NC_000913	4,128,812	0	G	.	100.0%	175.3 / NA	38	coding (141/1161 nt)	metB	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (17/21);  total (17/21)

ACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACA  >  NC_000913/4128571‑4129054                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                   aCGCAAAGAAGTTTAGATGTCCAGATGTATTGCAGTCCATTAACACAATGTTTCCTCTGGTGCCTGACATTTAACCGACACAGACCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCa                                                                                                                                                                                                                                          >  1:276911/1‑250 (MQ=255)   gCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCaccc                                                                                                                                                                                                                                       >  2:102473/1‑248 (MQ=255)       agaagTTTAGATGTCCAGATGTATTGACGGCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAAc                                                                                                                                                                                                                                   <  2:560/251‑1 (MQ=255)                                    tCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGc                                                                                                                                                                                                        >  2:153735/1‑249 (MQ=255)                                     ccATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGa                                                                                                                                                                                                     <  1:164147/251‑1 (MQ=255)                                                           ggTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggt                                                                                                                                                                                  <  2:205489/248‑1 (MQ=255)                                                                   aCATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtacttac                                                                                                                                                                       <  2:51402/251‑1 (MQ=255)                                                                        tCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAAt                                                                                                                                                                   <  1:41298/250‑1 (MQ=255)                                                                          aCCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATAcc                                                                                                                                                                >  2:229861/1‑251 (MQ=255)                                                                               cataGCCCAGGGAACTTCATCACATGCCGCGTAAACAGGCCACCATCGCAGTGCGTAGCTGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCAt                                                                                                                                                           <  2:155597/248‑1 (MQ=255)                                                                               cAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCAt                                                                                                                                                           <  1:123420/251‑1 (MQ=255)                                                                                aaaGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATg                                                                                                                                                          >  2:116199/1‑251 (MQ=255)                                                                                aaaGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTCGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATg                                                                                                                                                          <  1:113330/251‑1 (MQ=255)                                                                                    cccAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCcg                                                                                                                                                      >  1:211833/1‑251 (MQ=255)                                                                                             cTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAACTGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCAcc                                                                                                                                             >  2:136769/1‑251 (MQ=255)                                                                                                    ccaTGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAAc                                                                                                                                      <  1:208747/250‑1 (MQ=255)                                                                                                    acaTGACGCGTAAACAGGCCACCATCGCAGCGCGTAGAGGGATAACTGACGACGAACAGTATGGATGCGTTGTCCCACCGATCCATATTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAAc                                                                                                                                      >  1:266451/1‑251 (MQ=255)                                                                                                        gACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtg                                                                                                                                                                                 <  1:50216/204‑1 (MQ=255)                                                                                                        gACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtg                                                                                                                                                                                 >  2:50216/1‑204 (MQ=255)                                                                                                                         ccaTCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTgg                                                                                                                 >  2:198946/1‑251 (MQ=255)                                                                                                                                 gTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGGGGTGCTGGTGCAGTACTTACTAATACCGGCAGGTCAGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATctgc                                                                                                         >  2:145786/1‑251 (MQ=255)                                                                                                                                           ggTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgcc                                                                                               <  1:134526/251‑1 (MQ=255)                                                                                                                                             ttAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgccgc                                                                                             >  1:277321/1‑251 (MQ=255)                                                                                                                                              tAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTCTAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgccgc                                                                                             >  1:277821/1‑250 (MQ=255)                                                                                                                                                         gaACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcg                                                                                 >  1:120650/1‑251 (MQ=255)                                                                                                                                                                     tGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGGAGccttc                                                                     >  1:213290/1‑247 (MQ=255)                                                                                                                                                                        gTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGcc                                                                  <  1:116199/251‑1 (MQ=255)                                                                                                                                                                        gTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCCATCGcc                                                                  >  1:57200/1‑251 (MQ=255)                                                                                                                                                                           gTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGcccgt                                                               >  1:31888/1‑248 (MQ=255)                                                                                                                                                                                 ccGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACa                                                         <  1:136769/251‑1 (MQ=255)                                                                                                                                                                                   gATCCATCTGTCCAGCACCTATAACTGTACCGGATTTAATGAACCGCGCGCGCATGAGTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCGCGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGt                                                       <  1:44328/251‑1 (MQ=255)                                                                                                                                                                                                aGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAAcgcg                                          <  2:231104/251‑1 (MQ=255)                                                                                                                                                                                                         aaCTTGACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCGCGCCTGCTACGTCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTAt                                  <  2:63952/250‑1 (MQ=255)                                                                                                                                                                                                                    aTTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGttgttt                      <  1:229861/251‑1 (MQ=255)                                                                                                                                                                                                                         aTGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGa                 <  1:102473/251‑1 (MQ=255)                                                                                                                                                                                                                                 cgcgcgcATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCg         <  1:198339/251‑1 (MQ=255)                                                                                                                                                                                                                                   cgcgcATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGt                                                                                                     <  1:121149/157‑1 (MQ=255)                                                                                                                                                                                                                                        aTGATTACT‑‑CGTCGCGGCAACCCAACGGGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCAGGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACa  <  2:149404/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                   ACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAAGGCGATGAACA  >  NC_000913/4128571‑4129054

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