Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,910,755	G→A	D257D (GAC→GAT)	pstS ←	phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,910,755	0	G	A	98.5%	243.3 / NA	68	D257D (GAC→GAT)	pstS	phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (39/27);  minor base C (1/0);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGTTTTTCCACTCTTCGTTCACTTTC  >  NC_000913/3910511‑3911001                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGc                                                                                                                                                                                                                                                   >  1:72104/1‑250 (MQ=255)        atTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTtg                                                                                                                                                                                                                                           >  1:559513/1‑251 (MQ=255)           ggTCTTCCACGCAGCGCGAACCTGCTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTtgctg                                                                                                                                                                                                                                        >  2:187220/1‑251 (MQ=255)                tCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCAtt                                                                                                                                                                                                                                    >  2:340060/1‑250 (MQ=255)                 ccACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTa                                                                                                                                                                                                                                   <  1:249951/250‑1 (MQ=255)                  cACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGc                                                                                                                                                                                                                                 <  1:703757/251‑1 (MQ=255)                     gCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGaa                                                                                                                                                                                                                              <  1:485153/251‑1 (MQ=255)                      cAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTTGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCGCCAATCTGCACCTTTTGCTGCATTAGCGAAg                                                                                                                                                                                                                             <  2:71573/251‑1 (MQ=255)                       aGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGt                                                                                                                                                                                                                            >  1:203689/1‑251 (MQ=255)                           cgAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGttttt                                                                                                                                                                                                                        >  2:150826/1‑250 (MQ=255)                              aCCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTttcttc                                                                                                                                                                                                                     >  1:571719/1‑251 (MQ=255)                                  gTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGGCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGGTTTCTgcggtc                                                                                                                                                                                                                 <  2:352087/251‑1 (MQ=255)                                    tCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTATTAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcggtcgg                                                                                                                                                                                                               >  2:334228/1‑251 (MQ=255)                                       aCTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACt                                                                                                                                                                                                            >  1:42513/1‑251 (MQ=255)                                                     cAGGCTGGCGTAATCCAGGTCGTTCACCTGTTTAGCCCCGGTTTTGTCCGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGCTCTTTGTGGATCAGAATGAACGTGGTAGAGGAAATAGGCCATGCCTCTTCGCCTTTCTGGTTGGTCAGATCATGAGCGAAGGTTTTGCTCAAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCAt                                                                                                                                                                                              >  1:48360/1‑251 (MQ=255)                                                      aGGCTGGCGTAATCCAGGTCGTTCTCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGTCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATc                                                                                                                                                                                             <  2:388935/251‑1 (MQ=255)                                                          tGGCGTAATCCAGGTCTTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCg                                                                                                                                                                                         >  2:570510/1‑251 (MQ=255)                                                              gTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCgg                                                                                                                                                                                        >  1:130350/1‑248 (MQ=255)                                                                     aGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGttt                                                                                                                                                                              <  1:586808/251‑1 (MQ=255)                                                                      ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTg                                                                                                                                                                             >  2:609387/1‑251 (MQ=255)                                                                       gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTgg                                                                                                                                                                            >  1:475778/1‑251 (MQ=255)                                                                          gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGgtgt                                                                                                                                                                         >  1:26908/1‑251 (MQ=255)                                                                                     tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCTCCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGgtt                                                                                                                                                               <  1:708452/250‑1 (MQ=255)                                                                                             gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCGTTGTGTATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTGTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGGCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTGAGTGGACGCCAGGTTGTTCTGCtt                                                                                                                                                      <  2:619203/251‑1 (MQ=255)                                                                                                    aCGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTaa                                                                                                                                               <  1:556795/251‑1 (MQ=255)                                                                                                          aGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCatat                                                                                                                                         <  1:187220/251‑1 (MQ=255)                                                                                                           gTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACAATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATAtt                                                                                                                                        >  1:285768/1‑251 (MQ=255)                                                                                                              gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCaa                                                                                                                                     <  1:449815/251‑1 (MQ=255)                                                                                                              gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCaa                                                                                                                                     <  2:248654/251‑1 (MQ=255)                                                                                                                agaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACa                                                                                                                                   <  1:340060/251‑1 (MQ=255)                                                                                                                agaaTTTCAGCACTTCTGTGCCTTGTTCCGGTTTCTTCTGATATTTGTAGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTGTGCTGCATTAGCGACGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGGTCAGTTTGTTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAAc                                                                                                                                    >  2:35676/1‑250 (MQ=255)                                                                                                                 gaaTTTCAGCACTGCTGTGCCTTGTTCTGGTTTCTTCTGATCTTGGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGGTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACAt                                                                                                                                  <  1:73023/251‑1 (MQ=255)                                                                                                                     ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCGGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATATGCCATGCATCGTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCGTTTGCTGCATTAGGGAAGTGTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGTTGTACGCTAGGTGGTCCTGCGTCGCGGAAGCATAGTCAACATAAc                                                                                                                               <  2:665932/250‑1 (MQ=255)                                                                                                                         gCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAAtt                                                                                                                          >  1:300310/1‑251 (MQ=255)                                                                                                                           aCTTCTTTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAATCATATTCAACATAACCAATTGc                                                                                                                        <  1:546731/251‑1 (MQ=255)                                                                                                                               cTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCAccc                                                                                                                    >  1:229665/1‑251 (MQ=255)                                                                                                                                 gtgCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCAcc                                                                                                                     <  2:149864/248‑1 (MQ=255)                                                                                                                                    ccTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGc                                                                                                                 >  1:527573/1‑249 (MQ=255)                                                                                                                                    ccTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGc                                                                                                                 >  1:452132/1‑249 (MQ=255)                                                                                                                                             ggTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAGAGGCCATGCATCTTCGCCTTTCTGGTTTGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCGGCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTga                                                                                                        <  1:283938/249‑1 (MQ=255)                                                                                                                                               tttCTTATGATCTTTGTGGATAAGAATGAACGTGGTAGAGGTAATAGGCCCTGCATCTTCGCCTTTCTGGATGGTAAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAAGATAACCAATTGCACCCGGCAGACGCTgaacg                                                                                                     >  2:640575/1‑250 (MQ=255)                                                                                                                                                           tttGTGGATCAGAATGAACGTGGAAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGACAGCGGCGATAc                                                                                        >  2:521690/1‑251 (MQ=255)                                                                                                                                                              gtgGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACcgt                                                                                     >  2:136351/1‑251 (MQ=255)                                                                                                                                                                    cAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTa                                                                                <  1:22200/250‑1 (MQ=255)                                                                                                                                                                                ggTAGCGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGTGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCCTTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGAGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTTCACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTAAATTTAACGccc                                                                    >  2:115960/1‑250 (MQ=255)                                                                                                                                                                                  tAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACACGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGccc                                                                    >  1:404437/1‑248 (MQ=255)                                                                                                                                                                                          aTAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCgg                                                         <  2:130350/251‑1 (MQ=255)                                                                                                                                                                                               ccATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGCCTAACCGGTTTACCATCAGCGGCGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCAtt                                                    >  2:372788/1‑251 (MQ=255)                                                                                                                                                                                               ccATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCAtt                                                    <  1:164936/251‑1 (MQ=255)                                                                                                                                                                                                 aTGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCGTTTGCTGCATTAGCGAAGTTGTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCAtttt                                                  <  2:568985/251‑1 (MQ=255)                                                                                                                                                                                                  tGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTCTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTAGTTACCTTTACCGCCCAGCCCGATCGGCCAtttt                                                  >  1:278297/1‑250 (MQ=255)                                                                                                                                                                                                     aTCTTCGCCTTTCTGGTTGGACAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCAAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCATCCGGCAGACGCTGAACGAACGCGACGATAACGTCGTTACCTTTACAGCCCAGACCGATCGGCCATTTTACgg                                              >  2:275709/1‑251 (MQ=255)                                                                                                                                                                                                           gCCTTTCTGGTTTGTCAGATCTTGAGCGAAGGTTTTGCTCCACTCTGCACCTTTTTCTGCATTAGCGAAGTTTTCTTCGATCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATCCCGTCTTTAACTTTACCGCCCAGACCGATCGGCCCTTTTACGGCAGAGc                                        >  2:213650/1‑251 (MQ=255)                                                                                                                                                                                                                  tggttggtCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTACCCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTAAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGATAGAGCCAGTAcc                                 >  1:392266/1‑251 (MQ=255)                                                                                                                                                                                                                  tggttggtCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTAcc                                 >  1:568989/1‑251 (MQ=255)                                                                                                                                                                                                                   ggttggtCAGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTAcc                                 <  2:571719/250‑1 (MQ=255)                                                                                                                                                                                                                           aGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACgttgtt                        >  1:396719/1‑251 (MQ=255)                                                                                                                                                                                                                           aGATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACgttgtt                        <  1:270230/251‑1 (MQ=255)                                                                                                                                                                                                                           aGATCCTGAGCGAAGATTTTGCTACAATCTGCACCTTTTGCTCCCTTCGAGAAGATTTCTTCGGTCGGAATAACTGGTTTACCATCAGCGGAGATAAGTTTGGAGTACGCCAGGTTGTTCTACTTCGCGTAAGCATATTCAACATAACCAATTGAACCTGGCAGACGCTGAACGAACGAGGCGATACACTCGTAACCTGTACCGCCCAAAACGATAGGCAGTTTTACGGTagtgacctcatcaacattgat                        >  2:376105/1‑232 (MQ=255)                                                                                                                                                                                                                            gATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCCCCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGttt                       <  2:527573/251‑1 (MQ=255)                                                                                                                                                                                                                            gATCCTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGttt                       >  1:552246/1‑251 (MQ=255)                                                                                                                                                                                                                                cTGAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGTTTTTcc                   >  1:499811/1‑251 (MQ=255)                                                                                                                                                                                                                                  gAGCGAAGGTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGGACCAACGTTGTTTTTCCa                  <  2:72104/250‑1 (MQ=255)                                                                                                                                                                                                                                     cGAAGGTTTAGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGTTTTTCCActct              >  1:159023/1‑251 (MQ=255)                                                                                                                                                                                                                                         ggTTTTGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGCTACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGTTTTTCCACTCTTCGt          >  1:59886/1‑251 (MQ=255)                                                                                                                                                                                                                                            tttGCTCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGTTTTTCCACTCTTCGTTc        >  1:385124/1‑250 (MQ=255)                                                                                                                                                                                                                                                 tCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGTTTTTCCACTCTTCGTTCACTTTc  >  2:225762/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tCCAATCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGTTTTTCCACTCTTCGTTCACTTTc  >  1:290579/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACCATCAGCGGAGATCAGTTTGGTGTACGCCAGGTTGTTCTGCTTCGCGTAAGCATATTCAACATAACCAATTGCACCCGGCAGACGCTGAACGAACGCGGCGATACCGTCGTTACCTTTACCGCCCAGACCGATCGGCCATTTTACGGTAGAGCCAGTACCAACGTTGTTTTTCCACTCTTCGTTCACTTTC  >  NC_000913/3910511‑3911001

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