Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,532,409	T→G	intergenic (‑758/+28)	sgcX ← / ← yjhP	putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,532,409	0	T	G	100.0%	268.0 / NA	71	intergenic (‑758/+28)	sgcX/yjhP	putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (47/24);  total (47/24)

CTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAG  >  NC_000913/4532166‑4532655                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGAAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGccagcc                                                                                                                                                                                                                                                  >  1:267096/1‑250 (MQ=255)   aaaCCTTATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTTGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTgccagcca                                                                                                                                                                                                                                                 <  2:99918/249‑1 (MQ=255)   aaaCCTGATACCGCTTCTTCCCCGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCAATGGCACTGTGCGCATCTGGTGCAGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATCCATCAGATCAAAGCAAAGTTAAAGCCTCGATATTGccagcc                                                                                                                                                                                                                                                  >  1:43007/1‑248 (MQ=255)    aaCCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGa                                                                                                                                                                                                                                               >  1:180042/1‑250 (MQ=255)     aCCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGaa                                                                                                                                                                                                                                              >  1:431649/1‑250 (MQ=255)            aCCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCggg                                                                                                                                                                                                                                      >  2:17029/1‑251 (MQ=255)                  tcttcCCAGGGTTTGAGATTTCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTTGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATACGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAAcc                                                                                                                                                                                                                                  <  2:243695/249‑1 (MQ=255)                      cccAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCtgttgt                                                                                                                                                                                                                            <  1:187627/251‑1 (MQ=255)                       ccAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTgttgtt                                                                                                                                                                                                                           >  1:152239/1‑251 (MQ=255)                          gggTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTAATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑ACCTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTAcc                                                                                                                                                                                                                        <  2:356145/251‑1 (MQ=255)                                gagaTTGCCGGTGCTGTCTCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCg                                                                                                                                                                                                                   <  2:77155/250‑1 (MQ=255)                                  gaTTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATg                                                                                                                                                                                                                >  1:342702/1‑251 (MQ=255)                                     tGCCGTTGCTGTCGCTGCATTACGCAAAATATTGCTCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGCCAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGACTATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGAAAAGCCTCGATATTGCCAGCCCGCAATCCCGGGAACCTGTTGTTACCGGGCGATgcgg                                                                                                                                                                                                             >  2:448757/1‑248 (MQ=255)                                         gTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCg                                                                                                                                                                                                          >  2:48953/1‑250 (MQ=255)                                          ttGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGa                                                                                                                                                                                                         >  1:325949/1‑250 (MQ=255)                                           tGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAac                                                                                                                                                                                                       >  1:364604/1‑251 (MQ=255)                                                 cGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTAGTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGcc                                                                                                                                                                                                  <  2:258369/250‑1 (MQ=255)                                                  gatGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACATGTTGTTACCGGGCGATGAGTGCGAACACGcccc                                                                                                                                                                                                >  2:298896/3‑251 (MQ=255)                                                     gCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATCTTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGTGATGAGTGCGAACACGcccccgc                                                                                                                                                                                             >  2:334807/1‑248 (MQ=255)                                                     gCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGc                                                                                                                                                                                             >  1:328142/1‑251 (MQ=255)                                                          aCGCACAATATTGATCACGGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACAACAGATCAAAGCAAAGGTAAAGCCTCGATATGGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAg                                                                                                                                                                                        <  2:425914/251‑1 (MQ=255)                                                           cGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGc                                                                                                                                                                                       >  2:147615/1‑251 (MQ=255)                                                           cGCACAATACTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGGTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCATATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCTGGAACCTGTTGTTGCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGc                                                                                                                                                                                       <  1:247340/251‑1 (MQ=255)                                                              acaATCTTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGACATTCCGCCAGCCTGACAG‑AACTGTTCATGTTGTGCCATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACATGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGAAtt                                                                                                                                                                                    >  2:242779/1‑251 (MQ=255)                                                                      gATCACTGGCACTGTGTGCATCTGGTGCCGACATCCGCATTCCGCCAGCCTGCCAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCcgcgcg                                                                                                                                                                            >  2:231391/1‑251 (MQ=255)                                                                         cACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGt                                                                                                                                                                           >  1:330772/1‑249 (MQ=255)                                                                                  tgtgtgCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATg                                                                                                                                                                                                                >  1:37261/1‑203 (MQ=255)                                                                                  tgtgtgCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATg                                                                                                                                                                                                                <  2:37261/203‑1 (MQ=255)                                                                                        catctgtggCCGACATCGGCATTCCGCAAGCCTTACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGTCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTc                                                                                                                                                          <  2:25002/243‑1 (MQ=255)                                                                                        catatGGTGCCGACATAGGAATTCCGCAAGCCTG‑CAGGCAATGATCATGTTGTGAAATCTGATGCTCCTTATGGAGACATCTTGATTACTAAACAAACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGCTATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGACAACGCCCCAGCCAAAGCATTCCCGCGCGTAGATCACGTAGCGTTTc                                                                                                                                                          >  2:286407/5‑251 (MQ=255)                                                                                        cATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAATGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTc                                                                                                                                                          >  1:116342/1‑251 (MQ=255)                                                                                               tGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGATCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGcgcg                                                                                                                                                    >  2:90665/1‑250 (MQ=255)                                                                                                  cGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGt                                                                                                                                                >  1:144997/1‑251 (MQ=255)                                                                                                     cATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTa                                                                                                                                             <  1:90665/251‑1 (MQ=255)                                                                                                                cGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTg                                                                                                                                  >  1:1280/1‑251 (MQ=255)                                                                                                                 gCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCgg                                                                                                                                                        <  1:263739/228‑1 (MQ=255)                                                                                                                 gCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCgg                                                                                                                                                        >  2:263739/1‑228 (MQ=255)                                                                                                                                      cATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCg                                                                                                             <  2:161799/251‑1 (MQ=255)                                                                                                                                         gttgtGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCa                                                                                                          <  1:48953/251‑1 (MQ=255)                                                                                                                                            gtgACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAgg                                                                                                        >  2:135701/1‑250 (MQ=255)                                                                                                                                             tgACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGAt                                                                                                      >  1:145334/1‑251 (MQ=255)                                                                                                                                             tgACATATGATGCCCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCACAGGTAAAGCCTAGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGAt                                                                                                      >  2:33653/1‑251 (MQ=255)                                                                                                                                                     gatgCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCa                                                                                                          >  1:318604/1‑239 (MQ=255)                                                                                                                                                     gatgCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCAGCGTCa                                                                                                          <  2:318604/239‑1 (MQ=255)                                                                                                                                                     gatgCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCGGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTc                                                                                               <  2:325949/250‑1 (MQ=255)                                                                                                                                                      atgCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCCAAGCAAAGGTAAAGCCTTGATATTGCCAGCCAGAAATCCCGGGAACCTGTTTTTCCCGGGCTATGAGTTCCAACACGCCCCAGCCAAAGCATTCCAGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGCCGAAGGCATCGTCAGGATTCGCTTcca                                                                                             <  1:266777/251‑1 (MQ=255)                                                                                                                                                           ccTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCGGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGCCTTCGGCGGCCAAGTCATCGTCAGGATTCGCTTccagcc                                                                                          >  2:333432/1‑249 (MQ=255)                                                                                                                                                           ccTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccagccag                                                                                        <  1:147615/251‑1 (MQ=255)                                                                                                                                                                  gTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGgcgc                                                                                  >  1:40510/1‑250 (MQ=255)                                                                                                                                                                         cTTGATTATTAAACATACTGAATATATATTGTGCAGGTGATACAACAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTCAGTCACCAAGCGTTTCGGCGCGATATTTAGTTCAGCCCTGACTTCGGCGGCGAAGTCATAGTCAGGATTCGCTTACAGCCAGCGGAGCaacgtc                                                                            >  1:447243/1‑244 (MQ=255)                                                                                                                                                                         cTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATAGTCAGGATTCGCTTCCAGCCCGCGGCGCATGGTc                                                                            >  2:361282/1‑249 (MQ=255)                                                                                                                                                                              ttAATAAACATACTGAATATGGATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTAACGGGCGATGAGTGCGAACACGCCCCAGCCAAAGAATTCCCGCGCGTAGGTAACGTAGCGTTTAGGCGCGATGTTTAGTTCGGCCCTGCCTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGACGCCTGTGCAGCCAtt                                                                     >  2:206413/1‑251 (MQ=255)                                                                                                                                                                                aataaACATACTGAATATGTATTGGGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAtttt                                                                   <  2:364604/251‑1 (MQ=255)                                                                                                                                                                                  aaaaCATACTGAATATGTATTGTGCAGGTGATACATCAGAGCAAAGCAAAGGTAAAGCCTCGATAGTGCCAGCCAGAAATCCCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATGCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGGCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGc                                                                 <  2:40510/250‑1 (MQ=255)                                                                                                                                                                                    aaCATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCgg                                                               >  1:332934/1‑251 (MQ=255)                                                                                                                                                                                     aCATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATAGTGCCAGCCAGAACTACCGGGAACATGTTGTTACCGGGCGATGCGTGCGCAAACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGAGGCGAAGTAATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCCGCCCTTTTGCGGc                                                              >  2:200094/1‑251 (MQ=255)                                                                                                                                                                                      cATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTAAGCCATTTTGCGGCt                                                             >  2:243003/1‑251 (MQ=255)                                                                                                                                                                                        aaaggAATATTTATTGTGCAGGTGATACATCAGATCAAAGAAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAATCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCATCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCTGCTTc                                                           <  2:382467/247‑1 (MQ=255)                                                                                                                                                                                                      tgtgCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGCTGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGACTGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGc                                             >  1:302331/1‑251 (MQ=255)                                                                                                                                                                                                       gagCAGGTGATACATCAGATGAAAGCAAAGGTAAAGCCTCGATAGTGCCAGCAAGAAATACCGGGAACTTGGTGTTACCGGGCGATGAGTGAGAACACGCCCGAGCAAAAGCATGCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTGTAGTTCGTCCCTGATTCCGGCGGCGAAGTCAGGGACAGGTTTCGCTTCCAGGCAGGGGGGCATGGTCAGCGATTTTGCGGATGCATTCCGGTCCCAgt                                             <  2:192602/248‑2 (MQ=255)                                                                                                                                                                                                                     tCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGcc                               >  1:315008/1‑250 (MQ=255)                                                                                                                                                                                                                       aGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGt                            >  1:112705/1‑251 (MQ=255)                                                                                                                                                                                                                           caaagcaaagGTAAAGCCTCGATATTGCCAGCCAGAAAGACCGGGAACATGTTGTTCCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGCCATCGTCAGGATTCGCTTCCACCCAGCGGCGCAAGGTCAGCCATTTTGCGGCTTCATACCGATCCCAGACTTCCTGGTCTGCCAGTAcc                         >  2:284387/1‑250 (MQ=255)                                                                                                                                                                                                                             aagcaaagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCAAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAtt                      >  1:324617/1‑251 (MQ=255)                                                                                                                                                                                                                                aaaagGTAAAGCCTCGATATTGCCAGCCAGAAAAACCGGGAACCTGTTGTGACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAttacc                   <  2:267096/250‑4 (MQ=255)                                                                                                                                                                                                                                 aaagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTcca                  <  1:135701/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTccacc                >  2:351830/1‑250 (MQ=255)                                                                                                                                                                                                                                        aagacTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTCCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCAGGGTCAGCCATTTTGCGGCTTTATACCGGGCCCAGCCTTCCTGGTCTGCCATTCACATTTCCACCACGTc           >  2:14053/5‑251 (MQ=255)                                                                                                                                                                                                                                          gCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCCGCCTTCCTGGTCTGCCAGTACCCTTTCCACCCCGTCAt         >  2:251949/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTGCGGCGGCGAAGTCAGCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATa        <  1:16311/251‑1 (MQ=255)                                                                                                                                                                                                                                               gATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCCTAGCCg    >  2:381816/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tatTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCgcg  >  1:412452/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAG‑AACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCCGAG  >  NC_000913/4532166‑4532655

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: ‑