Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,907,139	G→T	A147E (GCG→GAG)	phoU ←	negative regulator of the pho regulon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,907,139	0	G	T	100.0%	222.2 / NA	64	A147E (GCG→GAG)	phoU	negative regulator of the pho regulon
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (44/20);  total (44/20)

CATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCATCGATCGCCA  >  NC_000913/3906894‑3907382                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    cATCGCCACCGACGTGCAGGAAATCCTGCCCCTTCACGTACTAGAAGATCAACTCACAAATATTCTGCCAGCGGTCGCCAATCCGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCGTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGTCTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAAc                                                                                                                                                                                                                                                >  2:93142/1‑251 (MQ=255) cATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAAc                                                                                                                                                                                                                                                >  1:371559/1‑251 (MQ=255)        aCCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTcc                                                                                                                                                                                                                                          >  1:121317/1‑250 (MQ=255)         ccGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAg                                                                                                                                                                                                                                        <  2:439712/251‑1 (MQ=255)            aCGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCAc                                                                                                                                                                                                                                     <  1:154085/251‑1 (MQ=255)              gTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCAcgt                                                                                                                                                                                                                                   >  1:200119/1‑251 (MQ=255)                     aaaTCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTgc                                                                                                                                                                                                                              >  1:388157/1‑249 (MQ=255)                                 ttCACGTAGTAGAAGATAAACTCAAAAATATTCTGGCCGCGCTCGCCAATACGTTCGATACAACGCGCACAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGCACAATACCTCCGTATTCCTGATCGACTTTTTTATCTTCCCGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCCGCATCTGGATGGt                                                                                                                                                                                                                >  2:208502/1‑251 (MQ=255)                                     cGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATg                                                                                                                                                                                                             >  1:66425/1‑250 (MQ=255)                                              agaTAATCTAACAAATATTCTAGCCGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAATAAGTACGATCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGAAGCCGAACAATACCTTCGTATTCCTGATCGCATTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGCCGGccccgc                                                                                                                                                                                                   >  2:190265/1‑248 (MQ=255)                                                      tCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTcc                                                                                                                                                                                             >  1:73999/1‑249 (MQ=255)                                                      tCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTcc                                                                                                                                                                                             >  1:107941/1‑249 (MQ=255)                                                       cacaAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCa                                                                                                                                                                                            <  1:11142/249‑1 (MQ=255)                                                         caAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGCACAATACCTTCGTATTCCTGCTCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTCTGAAGGCCCAGCGACTCCAGACt                                                                                                                                                                                        >  2:108909/1‑251 (MQ=255)                                                                  tGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACAAGCa                                                                                                                                                                                >  2:389780/1‑250 (MQ=255)                                                                      aGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGc                                                                                                                                                                           >  2:54935/1‑251 (MQ=255)                                                                       gCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCt                                                                                                                                                                          <  1:42397/251‑1 (MQ=255)                                                                         ggTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTg                                                                                                                                                                         >  1:230346/1‑250 (MQ=255)                                                                         ggTCGCCAATAAGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGCCTGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATACATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCCTCTGGATGGTATGCCGGCCCCGCGACTCCAGACTTACCTGCAACGGCTgg                                                                                                                                                                        >  2:37808/1‑250 (MQ=255)                                                                              ccAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGAtgctg                                                                                                                                                                   >  1:278639/1‑251 (MQ=255)                                                                                aaTACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGAGAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCATCTGTTCCAGGCTTTCCAGCAACGGCTGATgctgct                                                                                                                                                                 <  2:213811/251‑1 (MQ=255)                                                                                 aTACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGctgctg                                                                                                                                                                >  1:168536/1‑251 (MQ=255)                                                                                         gATAGAACGCGCTCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTAGGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAtt                                                                                                                                                        <  2:121317/251‑1 (MQ=255)                                                                                                cgcgcgCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTcc                                                                                                                                                   >  1:59841/1‑249 (MQ=255)                                                                                                cgcgcgCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAg                                                                                                                                                 <  2:110591/251‑1 (MQ=255)                                                                                                       agacaaGCGCAGTAAGTTCGCTCGTAATGGTACGCGAATCTTCCATCATGTATGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTATTTTATCTTCACGATAAACACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCTAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTa                                                                                                                                          <  2:301913/246‑1 (MQ=255)                                                                                                       aGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGt                                                                                                                                           >  1:432710/1‑250 (MQ=255)                                                                                                       aGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTa                                                                                                                                          <  1:405688/251‑1 (MQ=255)                                                                                                            aGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGc                                                                                                                                      >  1:39516/1‑250 (MQ=255)                                                                                                              cgcAGTAAGTTCGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTTCCAGCAACGGCTTATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGa                                                                                                                                   <  2:389882/251‑1 (MQ=255)                                                                                                               gcAGTCAGTACGCTCGGAATGGTACGCGAATCTTCCATAATGTAGGTCATCAGTTGACGCACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTAAAGGTGCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCACCTGAATGGTATGACGGCCCAGCGACCCCAGACTTAACAGCAACGGCTGCTGCTGCTGGGCGAATTTCACCAGACCAGTACGGc                                                                                                                                      >  2:148413/1‑247 (MQ=255)                                                                                                                          gCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGccc                                                                                                                       >  1:409654/1‑250 (MQ=255)                                                                                                                             cGGAATGGTACGCGAATCTTCCATCATGTAGGGCATCAGTTGACGAACAATACCGTCGTATTCCTGATCGACTTTTTTAGCTTCACGATAAATACGTACCGCTTCGTCAATTTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACg                                                                                                                     <  2:160477/250‑1 (MQ=255)                                                                                                                                                tCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCa                                                                                                     <  1:200274/247‑1 (MQ=255)                                                                                                                                                    tcatGTAGGTCATCAGTTGACGAACAATACCTTCGTATTACTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGc                                                                                              >  1:279701/1‑250 (MQ=255)                                                                                                                                                     catGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCa                                                                                             <  1:358162/250‑1 (MQ=255)                                                                                                                                                        gTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTAACGCTTCGTCAATGTCCATCCGCTCGCAAGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTCCGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATgg                                                                                         >  1:208191/1‑251 (MQ=255)                                                                                                                                                          aGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGCTTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGtt                                                                                       >  1:315055/1‑251 (MQ=255)                                                                                                                                                          aGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGtt                                                                                       >  1:68322/1‑251 (MQ=255)                                                                                                                                                                      gACGAACAATACCTTCGTAGTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAAGGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCGGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCAc                                                                            <  2:66425/250‑1 (MQ=255)                                                                                                                                                                         gAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAAAACGCTCCAGCTCGGCAATGGTTTTACTGATCACCAt                                                                         >  2:131149/1‑250 (MQ=255)                                                                                                                                                                         gAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTATCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATAATCTCCAGCTCGGCAATGGTTTTACttgtcaccct                                                                         >  2:229632/1‑241 (MQ=255)                                                                                                                                                                          aaCAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATa                                                                        >  1:86305/1‑250 (MQ=255)                                                                                                                                                                                   ttCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCACTGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCCCGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGCTCACCATAACCCGTCGc                                                               >  2:96121/1‑250 (MQ=255)                                                                                                                                                                                        aTTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGtcg                                                         >  1:366022/1‑251 (MQ=255)                                                                                                                                                                                             tGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGtcgctcgc                                                    >  1:184599/1‑251 (MQ=255)                                                                                                                                                                                              gATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGcc                                                   >  2:186043/1‑251 (MQ=255)                                                                                                                                                                                                 cGACTTTTTTAGCGTCACGAGAAATACGTACCGCTGCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGACGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCTATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGt                                                 <  2:412424/250‑1 (MQ=255)                                                                                                                                                                                                     ttttttATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCt                                            >  2:337390/1‑251 (MQ=255)                                                                                                                                                                                                       ttttATCTTCACGATAAATACGTAACGCTTCGTCACTGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGCCTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGGTGACAATACGCCCCCGCTCGGCAACGGTTTTACTGCTCACCATAACCAGTCGCAAGGCGCTCGCCGTCGGCTg                                           >  2:244786/1‑250 (MQ=255)                                                                                                                                                                                                               taaCGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGCCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCACTCGCAGGTCCCTCGCCGTCGGCTGAAGGTTGGc                                  >  2:353387/3‑251 (MQ=255)                                                                                                                                                                                                               tCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTGCGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTg                                    <  1:354219/249‑1 (MQ=255)                                                                                                                                                                                                                cACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCg                                 >  1:225551/1‑251 (MQ=255)                                                                                                                                                                                                                cACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCg                                 >  1:138323/1‑251 (MQ=255)                                                                                                                                                                                                                cACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGGTTGGCg                                 <  2:39516/251‑1 (MQ=255)                                                                                                                                                                                                                        aTACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGCCGTTTGGCGATAATgcg                         >  1:176736/1‑251 (MQ=255)                                                                                                                                                                                                                            gTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGcgcacg                     >  1:212528/1‑251 (MQ=255)                                                                                                                                                                                                                                 gCTTCGTCCATGACCATCCGCTCGAAAGCGTCCAGCACGTCGTGCAGCATCTGGAAGGTATGACGGCCCAGCGACTCCAGACTTGCCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGCTTTTGTCCGCCACGGCGCCAATACGATCCAGCTCGGCAATGGGTTTAATTATCACAACAACCAGTAGCAGGTCGCTCACCGTCGGCTGACGCTTGGCGCGACTGCGCACGCAggg                >  2:227859/1‑250 (MQ=255)                                                                                                                                                                                                                                    gcGTCAATGTCCATCCGCTCGAACGCGACCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCtt              <  2:8396/249‑1 (MQ=255)                                                                                                                                                                                                                                      gTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCat           >  1:324151/1‑251 (MQ=255)                                                                                                                                                                                                                                       tCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCatc          >  1:345415/1‑251 (MQ=255)                                                                                                                                                                                                                                        cAATGTCCATCCGCTCGAACGCGGCCAGCACGTCGTGCAGCATCTGGTTGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGATAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCatc          <  1:190265/250‑1 (MQ=255)                                                                                                                                                                                                                                          aTGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCatcgat       >  1:154265/1‑251 (MQ=255)                                                                                                                                                                                                                                               cATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCATCGATCGCCa  <  1:54935/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCATCGATCGCCA  >  NC_000913/3906894‑3907382

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