Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,556,984	G→A	T105I (ACC→ATC)	eutC ←	ethanolamine ammonia‑lyase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,556,984	0	G	A	96.4%	176.1 / ‑6.3	56	T105I (ACC→ATC)	eutC	ethanolamine ammonia‑lyase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/27);  minor base C (2/0);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGATCGCATCAGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGCTTCTGACGGC  >  NC_000913/2556740‑2557229                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      cATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTAATCGCATCAGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAATCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCtt                                                                                                                                                                                                                                                 <  1:318346/251‑1 (MQ=255)    cAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGATCGCATCAGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCg                                                                                                                                                                                                                                              <  2:158468/251‑1 (MQ=255)                 aTCTCTTCGTAGTTCACGGTGATCGCATCAGTTGACAGGCCATCCGTACTAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCCGAGCGTCCTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGCTAACTTTCGCACGAGAGTGATc                                                                                                                                                                                                                                 >  2:35333/1‑251 (MQ=255)                     ctTCGTAGTTCACGGTGATCGCATCAGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCgg                                                                                                                                                                                                                               <  1:3138/249‑1 (MQ=255)                       tCGTAGTTCACGGTGATCGCATCAGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATAGCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTTCGATATCTTTCGAACGAGAGTGATCGGc                                                                                                                                                                                                                              <  2:68191/248‑1 (MQ=255)                              tCACGGTGATCGCATCAGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGAGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGGGATTTTTGTCGCTAATTTAAGCGCGTACTTCCAGCAAACCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGCGTGATCGGCCAGGAAACgc                                                                                                                                                                                                                    >  1:38459/1‑251 (MQ=255)                                aCGGTGATCGCAGCAGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCGTTCACCCACTCTTCCGGCACTTCGTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACgcag                                                                                                                                                                                                                  <  1:283629/251‑1 (MQ=255)                                        cGCATCAGTTGACAGGCCAGCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGACGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCt                                                                                                                                                                                                          <  1:188393/251‑1 (MQ=255)                                              aGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATGGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTGTCACCCACGCTTCCGGCACGGCTTCCAGAACGATATCTTTCGAACGATAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTAc                                                                                                                                                                                                    <  2:16780/251‑1 (MQ=255)                                                  gACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGagag                                                                                                                                                                                                                                      >  1:110015/1‑213 (MQ=255)                                                  gACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGagag                                                                                                                                                                                                                                      <  2:110015/213‑1 (MQ=255)                                                             aGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCgg                                                                                                                                                                                     >  1:184727/1‑251 (MQ=255)                                                                    aCGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGAcc                                                                                                                                                                              >  1:176775/1‑251 (MQ=255)                                                                     cGACCTGAACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGGCAGGTAGAGCTTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGGGAGTGATCGGCCAGGAAACGCAGCAGCGGCTGCGTCCGTGGACGCGGACCGGCACCACCg                                                                                                                                                                             >  2:314365/1‑251 (MQ=255)                                                                       aCCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGt                                                                                                                                                                           <  1:208868/251‑1 (MQ=255)                                                                                         gCAACTCACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCCTACCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTCCCGGCACTTCTTTCAGTACGATATCTTGCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTTCGTACGTAGACGCGGACCGGCACGACCGGTACAAACGCTGGCCACGGt                                                                                                                                                         <  2:320094/251‑1 (MQ=255)                                                                                         gCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGt                                                                                                                                                         <  1:1281/251‑1 (MQ=255)                                                                                           aaCACACTGCGCTTTCAGCGCTTCAACAGCTTTTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGGAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGAGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGCCCGGTACAACCGCGGGCCACggagc                                                                                                                                                       >  1:286982/1‑248 (MQ=255)                                                                                            acacacTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACCCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGGACAAACGCGGGCCACGGtgct                                                                                                                                                      >  2:109162/1‑251 (MQ=255)                                                                                             cacacTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCCCGGtgctg                                                                                                                                                     >  2:157478/1‑251 (MQ=255)                                                                                              acacTGCGCTTTCAGAGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGGGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCCCGACCGGTACAAACGCGGGCCACGGtgctcc                                                                                                                                                    >  2:212075/1‑249 (MQ=255)                                                                                                    cgcTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCa                                                                                                                                              <  1:153101/251‑1 (MQ=255)                                                                                                        ttCAGCGCTTCAACAGCTTCTGCACCCAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGtt                                                                                                                                           <  2:104997/250‑1 (MQ=255)                                                                                                                      aGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTg                                                                                                                            <  2:229722/251‑1 (MQ=255)                                                                                                                          tCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTgcg                                                                                                                          <  1:164898/249‑1 (MQ=255)                                                                                                                                             ccATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACAcc                                                                                                      >  1:70433/1‑250 (MQ=255)                                                                                                                                              cATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCaa                                                                                                    <  1:68956/251‑1 (MQ=255)                                                                                                                                                 atCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCTCTTCGTTCAGTACGATATCTTTCGAACGAGAGTGAGCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATcc                                                                                                 <  2:73874/251‑1 (MQ=255)                                                                                                                                                   ccGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCAc                                                                                               <  1:195940/251‑1 (MQ=255)                                                                                                                                                           gTCAGGAAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGATAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGCCAGGAACAAACGCGGGCCACGGTGCTGCGGAGCAGGTCTGTTTATAAGTCTGCGCGATGCGGCTTTTCAAAACCAATCCACGCtttagg                                                                                       >  1:273048/1‑248 (MQ=255)                                                                                                                                                                        ttttGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCt                                                                          >  1:118851/1‑251 (MQ=255)                                                                                                                                                                        ttgtGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCGCTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCATTTCTGTTAATACGTCTGCGCGATGCTGATTTTCAACACCAATCCACGCTGTTGCTTCAGCGGatcc                                                                           <  2:38459/247‑4 (MQ=255)                                                                                                                                                                          ttGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTa                                                                         >  1:260983/1‑250 (MQ=255)                                                                                                                                                                                   aTTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAgcgc                                                               >  2:306603/1‑251 (MQ=255)                                                                                                                                                                                     ttCAGAGCGTATTTCCAGCCAACCCTGCGCTTACACACACTCTTCCGACACTTCTTTCAGTACGATATCTTTCGACCGAGAGTGATCGGCCAGGAACCGCCGCAGCGCCTGCGTACGTGGCCGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCCCCACCAACCCACGCTTTTGCTTCAGCGGAACCCAAATCCAGCGcac                                                             >  2:189735/1‑251 (MQ=255)                                                                                                                                                                                          agCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTGCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCtt                                                        <  2:70433/251‑1 (MQ=255)                                                                                                                                                                                            cGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTc                                                      >  1:68447/1‑251 (MQ=255)                                                                                                                                                                                             gTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTcc                                                     <  1:157478/251‑1 (MQ=255)                                                                                                                                                                                                cTTCCAGAAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGGGCTGCGGCGCAGTTCTGTTAATACGTCTGAGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCg                                                    >  2:239609/1‑249 (MQ=255)                                                                                                                                                                                                       cAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGCACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCg                                           >  2:144248/1‑251 (MQ=255)                                                                                                                                                                                                        aaGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGCGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGCTTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCCCCgg                                          >  2:28391/1‑251 (MQ=255)                                                                                                                                                                                                          gCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTg                                        >  2:232544/1‑251 (MQ=255)                                                                                                                                                                                                             cTGCGCTTTCACCCACTCTTCCGTCACTTCTTTCAGTACGCTATCTTTCGACCGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCg                                     >  2:220945/1‑251 (MQ=255)                                                                                                                                                                                                                 gcTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGcac                                  >  1:320730/1‑250 (MQ=255)                                                                                                                                                                                                                    ttCACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGtt                              >  2:128497/1‑251 (MQ=255)                                                                                                                                                                                                                      cACCCACTCTTCCGGCACTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTtgg                            <  2:224042/251‑1 (MQ=255)                                                                                                                                                                                                                      cACCCACTCTTCCGGCACGTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTtgg                            <  2:176775/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggCACTTCTTTCAGTACGATATCTCTCGAACGAGAGGGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACtt               <  1:208793/251‑1 (MQ=255)                                                                                                                                                                                                                                    gCACTTATTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACCGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACttt              >  2:210594/1‑251 (MQ=255)                                                                                                                                                                                                                                      aCTTCTTTCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGc            >  1:242634/1‑251 (MQ=255)                                                                                                                                                                                                                                           tttCAGTACGATATCTTTCGACCGAGAGTGATAGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGACCGCCCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTACTAAGTCTGCGCGATGCGGATTTTCAACACCAATCAACGATTTTGCTTACGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGCTTCTg       >  2:312948/1‑251 (MQ=255)                                                                                                                                                                                                                                           tttCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGCTTCTg       >  1:22780/1‑251 (MQ=255)                                                                                                                                                                                                                                            ttCAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGCTTCTg       >  1:197049/1‑250 (MQ=255)                                                                                                                                                                                                                                             tCAGTACGATAGCTTTCGAACGAGAGTGAGCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGCTTCTGa      <  2:229395/250‑1 (MQ=255)                                                                                                                                                                                                                                              cAGTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCGTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGCTTCTGAc     <  2:320730/250‑1 (MQ=255)                                                                                                                                                                                                                                                gTACGATATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGCTTCTGACGgc  <  1:135612/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CATCAGCGGTGGCAGGATCTCTTCGTAGTTCACGGTGATCGCATCAGTTGACAGGCCATCAGAAATAACGACCTGTACGTCCGGATTGGCAACACACTGCGCTTTCAGCGCTTCAACAGCTTCTGCACACAGGCGGCGGCCCATATCCGGGCGAGTCAGGTAGAGATTTTTGTCGCTAATTTCAGAGCGTACTTCCAGCAAGCCCTGCGCTTTCACCCACTCTTCCGGCACTTCTTTCAGTACGGTATCTTTCGAACGAGAGTGATCGGCCAGGAAACGCAGCAGCGCCTGCGTACGTGGACGCGGACCGGCACGACCGGTACAAACGCGGGCCACGGTGCTGCGGCGCAGTTCTGTTAATACGTCTGCGCGATGCGGATTTTCAACACCAATCCACGCTTTTGCTTCAGCGGAACCTAAATCCAGCGCACAGCTTTCCGAGGTCACCGGTGCCGCACAGTTGGTGGTGGCGCACTTTGCTTCTGACGGC  >  NC_000913/2556740‑2557229

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