Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,441,065	Δ1 bp	coding (611/621 nt)	rpsD ←	30S ribosomal subunit protein S4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,441,065	0	T	.	100.0%	310.9 / NA	67	coding (611/621 nt)	rpsD	30S ribosomal subunit protein S4
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (36/31);  total (36/31)

TGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACGACCGTTTACCATA  >  NC_000913/3440822‑3441307                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   ggCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGGCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTcg                                                                                                                                                                                                                                            <  2:396435/250‑1 (MQ=255)       aCTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGAt                                                                                                                                                                                                                                       <  2:513505/250‑1 (MQ=255)         tCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATc                                                                                                                                                                                                                                      >  2:211491/1‑249 (MQ=255)                  acacCATCAATCTCAACCGCGGTCACCGCTGAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGACACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCTAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGACCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGGGTTCGtt                                                                                                                                                                                                                           <  2:414917/251‑1 (MQ=255)                    acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAACCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTCGAAACTCTGTCACAGAACCCTGCATTGTGGCCTCTCTTTGGAACTAAGCTTTACTTGGAG‑CAAGCTCGACGATCAGGTGTTCGTTa                                                                                                                                                                                                                          >  1:183433/1‑250 (MQ=255)                         cAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAAAACGGCGCAGTGCGCTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGACCCCTGCATTGTGTCCTCTCTTTGGTGCTAAGCTTTACTTGGAG‑AAAGCTAGACGATCAGGTGTTCGTTAATGTcc                                                                                                                                                                                                                    <  1:545338/251‑1 (MQ=255)                                aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcagaca                                                                                                                                                                                                             <  1:442915/251‑1 (MQ=255)                                 acctgggtaacagcgcAACCCGGCATCGATGAGAGAAGAATACGGCGCAGTCAGTTACCCAGAGAATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCACGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGGCAACGATACCTGCAATGTGTCCCCTCTTTGGTACTAATCTATACTTTGAG‑AAAGATCGACACTTATGGGGTCGTTCATGTCCGcagac                                                                                                                                                                                                              >  2:317294/13‑249 (MQ=255)                                     cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTAAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGCTATCAACCAGGCGCGGTATTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAAGTGTTCGTTAATGTCCGCAGACAGATc                                                                                                                                                                                                         >  2:83445/1‑250 (MQ=255)                                         cACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGTCGATCAGGTGTTCGTTAATGTACGCAGACAGATCAGaa                                                                                                                                                                                                     <  2:317955/250‑1 (MQ=255)                                           ccGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCt                                                                                                                                                                                                                                              >  1:189874/1‑207 (MQ=255)                                           ccGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTGTGGTACTAAGCTTTACTTGGAG‑AAAGCt                                                                                                                                                                                                                                              <  2:189874/207‑1 (MQ=255)                                             gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGGATGGCCAAATCCACGCTCTAAAGGCTCAAGGGGCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTTTCCTCTCTTTGGTGCTAAGCGTTACTTGGAG‑AAAGCTCGACTATCAGTTATTCGTTAATGGCCGCAGTCAGATCAGAACGCTc                                                                                                                                                                                                <  2:516690/251‑1 (MQ=255)                                                     cGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTAc                                                                                                                                                                                        >  1:373572/1‑251 (MQ=255)                                                        cATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCt                                                                                                                                                                                     <  1:487834/251‑1 (MQ=255)                                                               gagaGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACg                                                                                                                                                                               <  1:244242/250‑1 (MQ=255)                                                                    cAGAATACGGCGCATTGCGGTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCCCCGTGGCGTGCGTCTAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGGCACAGAACCCTGCATGGTGTCCTCTCTGTGGTACTAAGCTGTACTTGGAG‑AAAGCGCTACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAGCGGTCCGGCTTACGCGTAAACGTACCtt                                                                                                                                                                         <  2:286961/251‑1 (MQ=255)                                                                    cAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCtt                                                                                                                                                                         <  1:84789/251‑1 (MQ=255)                                                                     aGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAACAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTc                                                                                                                                                                        >  1:223391/1‑251 (MQ=255)                                                                        aTACGGCGCAGTGCGTTACCCAGAGGATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCAt                                                                                                                                                                     <  2:533177/251‑1 (MQ=255)                                                                            ggCGCAGTGCGTTGACCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGTCGTGCAACGAACTCACTTGCTCGATATCAAACAGGCGCGGTTTTAGGAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCGTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCCGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTTCGGCTTGCGCTTAAACGTCCCTTCCATcgtg                                                                                                                                                                 <  2:118750/251‑4 (MQ=255)                                                                            ggCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTg                                                                                                                                                                 <  1:428712/251‑1 (MQ=255)                                                                              cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGcc                                                                                                                                                               >  2:117875/1‑251 (MQ=255)                                                                                        tACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAAc                                                                                                                                                      <  1:117875/250‑1 (MQ=255)                                                                                         aCCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACtt                                                                                                                                                    >  2:402846/1‑251 (MQ=255)                                                                                             agagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcc                                                                                                                                                  >  1:56365/1‑249 (MQ=255)                                                                                                        aaaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAAc                                                                                                                                                                                <  2:80078/208‑1 (MQ=255)                                                                                                        aaaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAAc                                                                                                                                                                                >  1:80078/1‑208 (MQ=255)                                                                                                              aCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTTCCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCAc                                                                                                                               <  2:30973/251‑1 (MQ=255)                                                                                                                     aaaGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTTCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTACACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTc                                                                                                                        >  1:93700/1‑251 (MQ=255)                                                                                                                      aaGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACCGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGCG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTAACGCTGCTcc                                                                                                                       >  2:468666/1‑250 (MQ=255)                                                                                                                      aaGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACGCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTca                                                                                                                       <  1:436413/251‑1 (MQ=255)                                                                                                                       aGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcag                                                                                                                      <  1:211520/251‑1 (MQ=255)                                                                                                                                ggTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCa                                                                                                             <  1:255979/251‑1 (MQ=255)                                                                                                                                   cACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCCGCTCAGCCAGCTCCaga                                                                                                           <  2:155899/250‑1 (MQ=255)                                                                                                                                         ggCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCt                                                                                                     <  1:141916/250‑1 (MQ=255)                                                                                                                                                gtcgaagtaaatcacTCGATATCACACAGACGCGGTTTTAGAAACTATCTCAAAGACCCCTGCATCCTCTCCTCTCTTTGGTACTGCGATTTACTTGGCG‑CCAGCTCGACGATCAGGTGTTCGTTCATGTCCGCAGACAGATCAGAACGCCCCGGCATACGCTTCAACCTACCTTCCATCTTACCAGCATCAACTTCCAGACCGGTTGACTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCAcgcg                                                                                             >  2:527596/15‑251 (MQ=255)                                                                                                                                                  cGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCAcgcg                                                                                             >  1:544884/1‑249 (MQ=255)                                                                                                                                                  cGAACTCACTTGCACGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCTCAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCgag                                                                                           <  1:302903/251‑1 (MQ=255)                                                                                                                                                    aaCTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGAc                                                                                         >  2:158710/1‑251 (MQ=255)                                                                                                                                                                cGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGGTGGCTTTTCACGCTGCGCATCCAGCTCCAGAGCGGCTTTTGCGCGAGACTGCTTCTTCg                                                                               <  2:88908/249‑1 (MQ=255)                                                                                                                                                                    atCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTtct                                                                         >  1:435687/1‑251 (MQ=255)                                                                                                                                                                      cAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCa                                                                       <  1:211491/251‑1 (MQ=255)                                                                                                                                                                       aaCCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCAc                                                                      >  1:157975/1‑251 (MQ=255)                                                                                                                                                                       aaCCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCAc                                                                      <  1:367326/251‑1 (MQ=255)                                                                                                                                                                              cgcgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGATCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCCCGAATGCt                                                               >  1:429713/1‑251 (MQ=255)                                                                                                                                                                                  gTTTTAGAAACTATGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑GAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAACCGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTAc                                                            >  1:129411/1‑250 (MQ=255)                                                                                                                                                                                       aGAAACTCTGTACCAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGTTCAGAACGCTCCGGCTTACGCTTAACCGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTc                                                      >  1:406574/1‑251 (MQ=255)                                                                                                                                                                                          caCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCAtt                                                   >  2:132965/2‑251 (MQ=255)                                                                                                                                                                                             tctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCgg                                                >  2:394924/1‑251 (MQ=255)                                                                                                                                                                                               tGTCACAGAACCCTCAATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTAAGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCAAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTAACGAATGCTTACACCGTCATTCGGAc                                              >  2:221805/1‑251 (MQ=255)                                                                                                                                                                                                gTCACAGAACCCTTCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCTGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCATGTTGGCTGTTCACGCTGCTCAGCCAGCTTCAGAGCGGCTTTCACGCGAGACGGCTTCTTCGCTTTCTCACGAATGCTTACAAAGTCATTCGGAc                                              <  2:495075/250‑1 (MQ=255)                                                                                                                                                                                                 tCACAGAACCCTGCATTGAGTCCTCTCTTTGGAACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCAAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGATTCTTCGCTTTCTCAAGAATGCTTACAACGTCATTCGGACt                                             >  1:476390/1‑250 (MQ=255)                                                                                                                                                                                                  cacaGAACCCTGCATTGTGTCCTCTCTTTGGGACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGGCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCAGCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCGTTCTCACGAAGGCTTACAACGTCAGTCGGACt                                             <  2:26451/249‑1 (MQ=255)                                                                                                                                                                                                         cccTGCATTGTGTACTATCTTTGGTACTAAGCTTTACTTGGAG‑CAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGATTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGAt                                     >  2:37445/1‑250 (MQ=255)                                                                                                                                                                                                          ccTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCCACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATa                                    >  2:124807/1‑250 (MQ=255)                                                                                                                                                                                                            tGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATa                                    >  1:289852/1‑248 (MQ=255)                                                                                                                                                                                                             gCATTGTGTCCTCACTTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGAGCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCATCTTCCAGCCATGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATaag                                  <  1:509660/249‑1 (MQ=255)                                                                                                                                                                                                             gCATTGGGTCCTCTCTTTGGTACTAAGCTTTACGTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCAGCTTGCCAGAATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTGCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATaagaa                                <  2:93700/251‑1 (MQ=255)                                                                                                                                                                                                                       ctctctTTGGCACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTa                       >  2:458142/1‑250 (MQ=255)                                                                                                                                                                                                                         ctctTTGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTaac                     >  1:309824/1‑250 (MQ=255)                                                                                                                                                                                                                             ttGGTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAacac                >  1:495565/1‑251 (MQ=255)                                                                                                                                                                                                                               ggTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACg               >  1:419376/1‑250 (MQ=255)                                                                                                                                                                                                                               ggTACTAAGCTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCAGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCACCTTACAGACAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACg               >  1:290808/1‑250 (MQ=255)                                                                                                                                                                                                                                        cTTTACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGCTCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACGACCGTTTAcc     >  1:271940/1‑251 (MQ=255)                                                                                                                                                                                                                                           tACTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTACTGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGATTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACGACCGTTTACCAt   >  2:72116/1‑250 (MQ=255)                                                                                                                                                                                                                                            aCTTGGAG‑AAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACGACCGTTTACCATa  >  1:548605/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   TGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTTTCACGCGAGACTGCTTCTTCGCTTTCTCACGAATGCTTACAACGTCATTCGGACTAACCTGATAAGAAGCGATGTTAACAACACGACCGTTTACCATA  >  NC_000913/3440822‑3441307

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