Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,359,314	G→A	Q240* (CAG→TAG)	cadB ←	lysine:cadaverine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,359,314	0	G	A	100.0%	167.1 / NA	46	Q240* (CAG→TAG)	cadB	lysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/25);  total (21/25)

TTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATCACCA  >  NC_000913/4359073‑4359559                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      ttCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACCCCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCat                                                                                                                                                                                                                                               <  2:16138/250‑2 (MQ=255)     gagtACCGTTGCTGGCGACTTCTCCATAAACTTTCGGGAAGTTCCCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACAAACATCACCCCGGAGCCCAGAGAGGTCAGGCTCGCAAAGGCGGGGAAGGCAGAAACCAGCGGCGCAGCCCAGTTCCCGAGGATAGTGGAAGCACTGATTGCAACCGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATg                                                                                                                                                                                                                                          <  1:233023/248‑1 (MQ=255)      aaTACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCCCCTAAGTCGCCGCGATGt                                                                                                                                                                                                                                         >  1:274750/1‑251 (MQ=255)       ataccgtggcttccgccTTCCCCATAAACTTTCGGAAAGTTACCGTCGTTAGCGGCACGTACTCCTGCCTGGCCTACAACCATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCATTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCCCCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTa                                                                                                                                                                                                                                        <  2:48215/235‑1 (MQ=255)                                       cGGGAAGTTACCGTCGTTAGCGGCACGTACCCCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCGGTTCCCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATg                                                                                                                                                                                                         <  1:80326/250‑1 (MQ=255)                                        gggAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACATAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGt                                                                                                                                                                                                       >  1:227954/1‑251 (MQ=255)                                            aGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGt                                                                                                                                                                                                       >  1:278889/1‑247 (MQ=255)                                            aGTTACCGTCGTGAGCGGCACTTACACCTGCCGGGCCTACCAACAGCAGCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTAAAGGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATGGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGt                                                                                                                                                                                                       <  2:278889/247‑1 (MQ=255)                                              atacggtcgtTCGCGGCACGTACACCTGTCTGGCCTCCCCACATCATCCAGGAGCCCAGAGCAGCCAGGCACGCACAGGCGATGAATGCAGAAACCAGCGGAGCAGACCAGTTACCGAGGCTAGTTGAAGCACTGCTTGCAAACGGAGCACCGGAAGCCGCCATTACAGCATACGGATAC‑TACCAGGA‑AAGAAACTAAGTCGCAGCGATGTAAACCATACCAGCTAAAcccatacacagccaggatgcg                                                                                                                                                                                                   >  2:29919/6‑230 (MQ=255)                                               tACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTCCCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAg                                                                                                                                                                                                 <  2:407795/250‑1 (MQ=255)                                                      gtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACgg                                                                                                                                                                                         >  1:448601/1‑251 (MQ=255)                                                         aGCGGCACGTACACCTGCCTGGCCTACCAACAGCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTAc                                                                                                                                                                                      <  1:335551/251‑1 (MQ=255)                                                           cGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGt                                                                                                                                                                                    >  1:133924/1‑251 (MQ=255)                                                           cGGCACGTACACCTGCCTGGCCTACCAACATCACCCAGGAGCCCAGAGAAGTCAGGCCCGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGt                                                                                                                                                                                    <  2:131190/251‑1 (MQ=255)                                                              cACGTCCACCTGCCTGGCCTACCAACATCACCCAGGAGCCAAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTc                                                                                                                                                                                 <  1:398520/251‑1 (MQ=255)                                                                    cacCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGttt                                                                                                                                                                           >  2:296056/1‑251 (MQ=255)                                                                                cTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATAcc                                                                                                                                                                >  1:302379/1‑250 (MQ=255)                                                                                   ccAACATCATCGAGGAGCCAAGAGAAGCCAGGCACGCAACGGCGGTGAATGCAGAAACCAGCGGCGCAGCACAGATACCGAGGATAGTTGAAGCACTAATTGCACACGGGGCACCGGCAGCCGCCATTACAGAAGACGGAAACATAC‑‑GGAGAAGCAGCTAGGTCGAAGCGATGTAAACAAAACCTGCAAAACCAATACCCGACATGGTTGAAAGGGGGACGTAACGATTCGGGTTTTTAACCATaacgta                                                                                                                                                             >  2:452225/1‑245 (MQ=255)                                                                                        atcatcCAGGAGCCCAGAGAAGCCAGGCCCGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTAc                                                                                                                                                       <  2:277837/251‑1 (MQ=255)                                                                                            cccAGGAGCCCAGAGAAGTCAGGCACGCAAAGGGGGTGAATGCAGAAACCAGCGGCGCAGCCCGGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCt                                                                                                                                                   <  2:127339/250‑1 (MQ=255)                                                                                               aGGAGCCCAGAGAAGTCCGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTACACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGCACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCg                                                                                                                                                >  2:397683/1‑251 (MQ=255)                                                                                                                   cACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATg                                                                                                                                                                                                         >  1:48670/1‑174 (MQ=255)                                                                                                                   cACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATg                                                                                                                                                                                                         <  2:48670/174‑1 (MQ=255)                                                                                                                       aaaaGGCTGTGAATGCGGAAACCAGCGGCGCAGCCCAGTTTCCGAGAAAATTGGAAGCAGTGATTGCAACCGGATCACGGGAAGCCGCCATTACAGAAGACGTATACATACC‑GGA‑AAGCATCGAAGTCGCGTCGATGTGAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGTAACGGTACGTTTCTTGTTTTTAACCATACCAGTACTTACAGCGGCTTATTCAACATCCACGAGGGCacac                                                                                                                        <  2:293659/250‑1 (MQ=255)                                                                                                                         aaGGCGGTGAATGCAGAAACCAGCGGCGCAGCCATGTAACCGAGGATAGTTAAAGCACTGAGTGCAAACGGAGCACCGAAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACTCCCACGAAGGCCCACAg                                                                                                                      <  2:227954/251‑1 (MQ=255)                                                                                                                                 gAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                                                                              <  1:296056/251‑1 (MQ=255)                                                                                                                                            ctaGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACtttt                                                                                                    >  2:115485/3‑250 (MQ=255)                                                                                                                                                      aGCCCAGTTCCCGAGGATAGTTGAAGCACTGATTGCAAACGGAGACCCGAAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGGTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACTCCCACGAAGGCCCACAGGCAGAGCAGAATACTTGTAATGATCGCAt                                                                                         <  2:340103/251‑1 (MQ=255)                                                                                                                                                          cAGTTACCGAGGATAGTTGAACCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGAcc                                                                                     <  2:179579/251‑1 (MQ=255)                                                                                                                                                           aGTTACCGAGGCTAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCCTTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAAt                                                                                                                                                                                                                                  >  2:239464/1‑109 (MQ=255)                                                                                                                                                           aGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAAt                                                                                                                                                                                                                                  <  1:239464/109‑1 (MQ=255)                                                                                                                                                           aGTTACCACGGATAGTTGAAGAACTGATTGCAAACGGAGCACCGGAAGCCGAAATTACAGAAGACGGATACATACC‑GGA‑AAGCAACTAAGTCGCAGAGATGTAAACAATACACGCTAAACCAGAACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGATGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGAcc                                                                                     >  2:232553/1‑250 (MQ=255)                                                                                                                                                             ttACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATc                                                                                  >  2:6787/1‑251 (MQ=255)                                                                                                                                                                  gAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTgg                                                                             >  2:12486/1‑251 (MQ=255)                                                                                                                                                                    ggATAGTTGAAGCCCTGATGGCAAACGGAGCACCGGATGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGGTCGCATGACCATCAGTGGTa                                                                           <  1:385092/251‑1 (MQ=255)                                                                                                                                                                                          aaaCGGAGCCCCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTg                                                     <  2:140584/251‑1 (MQ=255)                                                                                                                                                                                          aaaCGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTg                                                     >  1:242884/1‑251 (MQ=255)                                                                                                                                                                                                        aaGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                                                                                 <  1:44486/193‑1 (MQ=255)                                                                                                                                                                                                        aaGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                                                                                 >  2:44486/1‑193 (MQ=255)                                                                                                                                                                                                         aGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCAt                                      <  2:133924/251‑1 (MQ=255)                                                                                                                                                                                                           ccgccATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTATTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCa                                    <  1:397683/251‑1 (MQ=255)                                                                                                                                                                                                            cgccaatcaagaagaCCGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCaa                                   <  2:144790/242‑1 (MQ=255)                                                                                                                                                                                                            cgccATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCa                                    >  2:131157/1‑250 (MQ=255)                                                                                                                                                                                                                             gATACATACC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCcaac                   >  2:106460/1‑250 (MQ=255)                                                                                                                                                                                                                                  atacC‑GGA‑AAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAACCCAATGCCAGCCAACAATAGc             >  2:360383/1‑251 (MQ=255)                                                                                                                                                                                                                                               cACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATcacca  <  1:75436/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      TTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACC‑GGA‑AAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATCACCA  >  NC_000913/4359073‑4359559

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