Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,933,621	A→G	L18L (CTA→CTG)	unknown →	Phosphate transport regulator (distant homolog of PhoU)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,933,621	0	A	G	100.0%	208.2 / NA	69	L18L (CTA→CTG)	unknown	Phosphate transport regulator (distant homolog of PhoU)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (41/28);  total (41/28)

AACTTTAATTCAATTTCATCCTGCATTATAAAATCCCTTAAAGATCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATCCTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTACAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTGACTCAACAAGATAAATTGGCAAA  >  1/1933378‑1933857                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             aaCTTTAATTCAATTTCATCCTGCATTATAAAATCCCTTAAAGATCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAac                                                                                                                                                                                                                                       <  1:805144/251‑1 (MQ=255)                        cATTATAAAATCCCTTAAAGATCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAAtgt                                                                                                                                                                                                                <  1:210972/251‑1 (MQ=255)                          ttATAAAATCCCTTAAAGATCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAAtgttg                                                                                                                                                                                                              >  1:419762/1‑251 (MQ=255)                            ataAAATCCCTTAAAGATCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAg                                                                                                                                                                                                                           <  1:53322/236‑1 (MQ=255)                                 aTCCCTTAAAGATCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAac                                                                                                                                                                                                                                       <  1:412991/219‑1 (MQ=255)                                        aaaGATCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATcc                                                                                                                                                                                                <  1:524897/251‑1 (MQ=255)                                             tCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCtttt                                                                                                                                                                                           >  1:734813/1‑251 (MQ=255)                                                   ggTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAAtgt                                                                                                                                                                                                                >  1:226886/1‑224 (MQ=255)                                                                                          ccAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGAc                                                                                                                                                                   <  1:673519/230‑1 (MQ=255)                                                                                             ataaataaGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAat                                                                                                                                                                                                         >  1:519207/1‑189 (MQ=255)                                                                                               aaataaGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAg                                                                                                                                                                                                                           <  1:685889/169‑1 (MQ=255)                                                                                                                                    cACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAAtgt                                                                                                                                                                                                                >  1:553105/1‑143 (MQ=255)                                                                                                                                         aCATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAat                                                                                                                                                                                                         >  1:100928/1‑145 (MQ=255)                                                                                                                                         aCATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAAttatta                                                                                                                                                                                                    <  1:238541/150‑1 (MQ=255)                                                                                                                                         aCATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAAttatta                                                                                                                                                                                                    <  1:731254/150‑1 (MQ=255)                                                                                                                                         aCATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACAttt                                                                                                                                                                            >  1:452548/1‑174 (MQ=255)                                                                                                                                         aCATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATAc                                                                                                                                                                       <  1:195927/179‑1 (MQ=255)                                                                                                                                         aCATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGata                                                                                                                                    <  1:277924/214‑1 (MQ=255)                                                                                                                                         aCATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTc                                                                                                                                >  1:456768/1‑218 (MQ=255)                                                                                                                                         aCATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACg                                                                                               >  1:329067/1‑251 (MQ=255)                                                                                                                                            tAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATAc                                                                                                                                                                       <  1:378721/176‑1 (MQ=255)                                                                                                                                             aaTGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAg                                                                                                                                                         <  1:507473/189‑1 (MQ=255)                                                                                                                                               tGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCt                                                                                                             >  1:213083/1‑231 (MQ=255)                                                                                                                                               tGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAAc                                                                                                >  1:777767/1‑244 (MQ=255)                                                                                                                                                          cgTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTaaaa                                                                               <  1:348290/250‑1 (MQ=255)                                                                                                                                                                ttattaACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAattat                                                                                                                                                                                                      <  1:216108/125‑1 (MQ=255)                                                                                                                                                                ttattaACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACAt                                                                                                                                                                              <  1:311832/149‑1 (MQ=255)                                                                                                                                                                ttattaACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATAc                                                                                                                                                                       <  1:93871/156‑1 (MQ=255)                                                                                                                                                                ttattaACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGttt                                                                                                                                         >  1:137436/1‑186 (MQ=255)                                                                                                                                                                ttattaACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACcg                                                                         >  1:507038/1‑250 (MQ=255)                                                                                                                                                                ttattaACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTGACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGGTAAAACGTGAGTTTCGTTTAAAATTACcg                                                                         >  1:176564/1‑250 (MQ=255)                                                                                                                                                                          tccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTAtt                                                              >  1:865429/1‑251 (MQ=255)                                                                                                                                                                          tccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTAtt                                                              >  1:513940/1‑251 (MQ=255)                                                                                                                                                                          tccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTAtt                                                              >  1:275042/1‑251 (MQ=255)                                                                                                                                                                          tccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTAtt                                                              >  1:244142/1‑251 (MQ=255)                                                                                                                                                                           ccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTg                                                                                                                                                                                       >  1:831096/1‑129 (MQ=255)                                                                                                                                                                           ccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAAt                                                                                                                                                                         <  1:205034/143‑1 (MQ=255)                                                                                                                                                                           ccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGAcc                                                                                                                                                                  >  1:95952/1‑149 (MQ=255)                                                                                                                                                                           ccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTAtt                                                              >  1:551898/1‑250 (MQ=255)                                                                                                                                                                           ccttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTAtt                                                              >  1:661512/1‑250 (MQ=255)                                                                                                                                                                            cttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCa                                                            <  1:381216/251‑1 (MQ=255)                                                                                                                                                                            cttAAAGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCa                                                            <  1:571222/251‑1 (MQ=255)                                                                                                                                                                                aaGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGcc                                                        <  1:295880/251‑1 (MQ=255)                                                                                                                                                                                aaGGAGTCACTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGcc                                                        <  1:60575/251‑1 (MQ=255)                                                                                                                                                                                        aCTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCa                                                            >  1:768888/1‑239 (MQ=255)                                                                                                                                                                                         cTTTTTATGGCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTatcgat                                                 >  1:657085/1‑249 (MQ=255)                                                                                                                                                                                                  gCAATGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGcc                                                        >  1:163650/1‑233 (MQ=255)                                                                                                                                                                                                     aTGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTg                                                                                                                              >  1:676676/1‑160 (MQ=255)                                                                                                                                                                                                     aTGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTatcgat                                                 <  1:121152/237‑1 (MQ=255)                                                                                                                                                                                                      tGAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATttgttg                                  <  1:165848/251‑1 (MQ=255)                                                                                                                                                                                                       gAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAg                                                                                                                                                         >  1:94344/1‑131 (MQ=255)                                                                                                                                                                                                       gAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTc                                                             <  1:424151/223‑1 (MQ=255)                                                                                                                                                                                                       gAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTatcgat                                                 >  1:117631/1‑235 (MQ=255)                                                                                                                                                                                                       gAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTgg                                 >  1:364308/1‑251 (MQ=255)                                                                                                                                                                                                       gAATAATATACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTgg                                 >  1:138112/1‑251 (MQ=255)                                                                                                                                                                                                         ataataTACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCa                                                            >  1:628350/1‑222 (MQ=255)                                                                                                                                                                                                          taataTACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGACTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAAt                              >  1:423082/1‑251 (MQ=255)                                                                                                                                                                                                           aataTACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGc                                                                                                                                                                                   >  1:332331/1‑101 (MQ=255)                                                                                                                                                                                                            atatACTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTa                            >  1:403865/1‑251 (MQ=255)                                                                                                                                                                                                                         attTGCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATAAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTGACTCAACAAg               <  1:523604/251‑1 (MQ=255)                                                                                                                                                                                                                             gCCCACTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATAc                                                                                                                                                                       <  1:573297/95‑1 (MQ=255)                                                                                                                                                                                                                                  cTCCCCTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTGACTCAACAAGATAAATTgg      >  1:837847/1‑251 (MQ=255)                                                                                                                                                                                                                                    ccccTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGcc                                                        >  1:405561/1‑199 (MQ=255)                                                                                                                                                                                                                                      ccTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCa                                                            >  1:187903/1‑193 (MQ=255)                                                                                                                                                                                                                                      ccTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTGACTCAACAAGATAAATTGGCaaa  >  1:462528/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTGACTCAACAAGATAAATTGGCaaa  >  1:857851/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTGACTCAACAAGATAAATTGGCaaa  >  1:184631/1‑251 (MQ=255)                                                                                                                                                                                                                                       cTTTAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATCAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTGACTCAACAAGATAAATTGGCaaa  >  1:224502/1‑250 (MQ=255)                                                                                                                                                                                                                                        tttAAAGCCTCTGCAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTg                         <  1:476416/226‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             AACTTTAATTCAATTTCATCCTGCATTATAAAATCCCTTAAAGATCATGAGGTTTTATCTTCTTACCTGATTTATAAGGGATTTATTATCCAATAAATAAGGGGATTTTTCGATTAAAATCGTGAATTTTTCACCAACATAATGTAATATGCGCGTGAATTATTAACCTTCCTTAAAGGAGTCACTTTTTATGGCAATGAATAATATCCTCGGATTATTTGCCCACTCCCCTTTAAAGCCTCTACAAAAACACTCTAATAAAGTTACCGAATGTTGTGAATTATTAATCCCTTTTTTTGAAGCAACATTTAATACCGACTGGGTAAAAGCGGAGCAAATTCGTTTAGATATTTCTGCACACGAACGTCAGGCTGACGCGTTAAAACGTGAGATTCGTTTAAAATTACCGCGCGGTTTATTCATGCCTATCGATCGTACGGATTTGTTGGAATTAGTGACTCAACAAGATAAATTGGCAAA  >  1/1933378‑1933857

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