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	97.6%	127.1 / ‑3.0	42	A147E (GCG→GAG)	phoU	negative regulator of the pho regulon
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (28/12);  total (30/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

CCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCATCGATC  >  NC_000913/3906899‑3907378                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 ccACCGACCTGACGGAAATCCTGCCCCTTCACGTAGTAGACGCTAAACTCACCAATATTCTGGCAGCGGTCGCCAATACGTTCGATCGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGACTCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGAACCGCTTCGTCAATGTCCATCCGCTCGACCGCGTc                                                                                                                                                                                                                                        >  2:128302/1‑251 (MQ=255)   aCCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTcc                                                                                                                                                                                                                                       >  1:174538/1‑250 (MQ=255)      gACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGATTCGTCAATGTCCATCCGCTCGAACTCGTCCAGc                                                                                                                                                                                                                                    >  2:92435/1‑250 (MQ=255)      gACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGCATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGc                                                                                                                                                                                                                                    >  1:195447/1‑250 (MQ=255)          tGACGGAAATCCTGCCCCTTCACGTAGTAGACGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTCAGTACGATCGGAATGGTACGCGGATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCAcgtc                                                                                                                                                                                                                               >  1:199329/1‑251 (MQ=255)                                                   acaAATATTCTGGCAGCGGGCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCAGCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGAc                                                                                                                                                                                      <  2:89840/251‑1 (MQ=255)                                                        tatTCTGGCAGCGGTCGCCCAACCGTTCGATAGAACGCGGGCAGAACAGCGCAGTAGATACGCTCGGAATGGTACGCGCGTCTTCCATACTGCAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAAAGCGGCCAGCACGTCGAGCAGCATCTGGATGGTATGATGCCCCAGCGAAGCCAGACTGAcc                                                                                                                                                                                 >  2:188691/1‑251 (MQ=255)                                                        tatTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTAcc                                                                                                                                                                                 >  2:40420/1‑251 (MQ=255)                                                                       cGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGAtgc                                                                                                                                                                  <  1:15040/251‑1 (MQ=255)                                                                        gCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGAtgct                                                                                                                                                                 >  2:28515/1‑251 (MQ=255)                                                                         ccAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGAtgctg                                                                                                                                                                >  2:132623/1‑251 (MQ=255)                                                                                  tCGATAGAACGCGCGCAGAACAGCGCAGTAAGTCCGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACGCCAGACTTACCAGCAACGGCTGATGCTGCTGGgaga                                                                                                                                                        <  2:31104/250‑1 (MQ=255)                                                                                   cGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAt                                                                                                                                                      >  1:240148/1‑251 (MQ=255)                                                                                      tAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑Ttc                                                                                                                                                  >  1:40095/1‑251 (MQ=255)                                                                                                     aCAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACAGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACgg                                                                                                                                   >  1:21809/1‑251 (MQ=255)                                                                                                         cgcAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGc                                                                                                                                  >  1:166367/1‑248 (MQ=255)                                                                                                                       tCGGAATGGTACGCGAATCTTCCCTCATGTAGGTCATCAGCTGCCGCACAATCCCTTCGAATTCCTGATCGACTTTTTTATCTTCACGCTACATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCCGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGCTTTTGTCCGCCACg                                                                                                                 >  2:167999/1‑251 (MQ=255)                                                                                                                                    cgAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGCTCGACTTTTTTATATTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCCGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCt                                                                                                    >  2:125310/1‑251 (MQ=255)                                                                                                                                    cgAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCt                                                                                                    <  2:166367/251‑1 (MQ=255)                                                                                                                                               tcatGTAGGTCATCAGTTGACGCACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCCCGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGcc                                                                                         >  1:227039/1‑250 (MQ=255)                                                                                                                                                      ggTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCCGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGttt                                                                                  >  2:91632/1‑251 (MQ=255)                                                                                                                                                          atcaGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACAGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACt                                                                              >  1:175524/1‑251 (MQ=255)                                                                                                                                                          atcaGATGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCAATCCGCTCGAACGCGTCCAGCACGTCGAGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTAACCAGCAACGGCGGATGCTGCGGGGA‑AATGTTCACCAGCACAGTAAGGCAAATTTTGTCCCCCACGTCGCCAATACGACCCAGCTCGGCCATGATTTTagt                                                                              >  2:217786/1‑249 (MQ=255)                                                                                                                                                             aGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGAt                                                                           >  1:65246/1‑251 (MQ=255)                                                                                                                                                               ttGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATcc                                                                         >  1:73850/1‑250 (MQ=255)                                                                                                                                                                tGACTAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCa                                                                         <  2:40095/250‑1 (MQ=255)                                                                                                                                                                   cGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAAGGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCAt                                                                     <  2:205974/251‑1 (MQ=255)                                                                                                                                                                    gAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCCGCTCGGCAATGGTTTTACTGATCACCAt                                                                     >  2:77665/1‑250 (MQ=255)                                                                                                                                                                    gAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATa                                                                    <  1:194156/251‑1 (MQ=255)                                                                                                                                                                           aCCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTc                                                             >  1:47436/1‑251 (MQ=255)                                                                                                                                                                                   aTTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGtcg                                                     >  1:239214/1‑251 (MQ=255)                                                                                                                                                                                                ttttttATCTTCACGATAAATACGTACCGCTTCGTCTCTATGCATTCGAACGAAAGCGTCCAGCACGTCGTGCAGCATCTGAATGGAATGACGGCCAAGCGACTCAAGAATTACCAGCAACAGCAGATGCTGCTGGGAAAAT‑TTCTCCAGAGCAGTACGGGAGATTTTGTCCGACACGTCGTCAACACCAGCCAGCTCGGCAATGGTTTTACTGAGCACCATAAACAgaacgacgacactcgccgtaggct                                        >  1:128093/1‑228 (MQ=255)                                                                                                                                                                                                   tttATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGa                                      <  2:164892/250‑1 (MQ=255)                                                                                                                                                                                                    ttATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGCTGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAAACAGTCGCAGGTCGCTCGCCGTCGGCTGACg                                    >  2:167463/1‑251 (MQ=255)                                                                                                                                                                                                        cTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATCACGGCCCCGCGACTCCAGCCTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCCGCGCAGTACGGCAGATTTTGTCCGCCACGGCGCCCATACGCTCCAGCTCGGCACTGGTTTTACTGATAACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGAAGTTTg                                >  2:276049/1‑251 (MQ=255)                                                                                                                                                                                                           cACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCg                             <  1:3374/251‑1 (MQ=255)                                                                                                                                                                                                                 aaaTACGTACCGCTTCGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATg                       <  1:56087/251‑1 (MQ=255)                                                                                                                                                                                                                    gaCGTACCGCTGCGTCAATGTCCATCCGCTCGAACGCGGCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTGTGGCGATAATgcgc                    <  2:84202/250‑1 (MQ=255)                                                                                                                                                                                                                                cGTCAATGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCACTGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGAGATAATGCGCACGCAGGCTTcc        >  2:242090/1‑250 (MQ=255)                                                                                                                                                                                                                                 gtcaatgtcaatCCGCTCTTACGCGTCGAGCACGACTTGCAGCATCTGGATGGTATGACGGGCCAGCGAATCCAGGCTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTACGCCACGTCGCCAATACGCCCCAGCTCGGCAATGGTTTTAATGATCACCATAACCAATCGAAGGTCGCTCGCAGTCCGATGACGTTTGGCGATAATGCGCACGCAGGCTTCat       >  2:34388/1‑251 (MQ=255)                                                                                                                                                                                                                                      tGTCCATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGCTCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGTCGATACTGCGCCCGCAGGCTTCttcgatc  >  2:5896/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggccATCCGCTCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCatcgatc  <  2:73850/248‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 CCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAT‑TTCTCCAGCGCAGTACGGCAGATTTTGTCCGCCACGTCGCCAATACGCTCCAGCTCGGCAATGGTTTTACTGATCACCATAACCAGTCGCAGGTCGCTCGCCGTCGGCTGACGTTTGGCGATAATGCGCACGCAGGCTTCATCGATC  >  NC_000913/3906899‑3907378

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