Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,789,654	T→C	D14D (GAT→GAC)	ydiE →	PF10636 family protein YdiE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,789,654	0	T	C	100.0%	219.7 / NA	61	D14D (GAT→GAC)	ydiE	PF10636 family protein YdiE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/24);  total (37/24)

ACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGC  >  NC_000913/1789408‑1789900                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       aCCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGcc                                                                                                                                                                                                                                                     >  1:152229/1‑250 (MQ=255)      tgtCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATcac                                                                                                                                                                                                                                               >  2:613449/1‑251 (MQ=255)          tgggctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATcaca                                                                                                                                                                                                                                              >  2:542530/1‑248 (MQ=255)           gggctgggAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGAc                                                                                                                                                                                                                                          <  2:266025/251‑1 (MQ=255)             gctgggAGGATACCGAACGCCTGGTCGAAAAACGCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCAGTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCg                                                                                                                                                                                                                                        <  1:273344/251‑1 (MQ=255)                 agaggaTACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCGTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATACGACTCATTGTTATATCAATTATTATTAATTTTTATGCGTTCTACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCATTAACAAGACCGCTTc                                                                                                                                                                                                                                    >  2:664449/2‑251 (MQ=255)                          cGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGcc                                                                                                                                                                                                                           >  2:86447/1‑251 (MQ=255)                                  tGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGc                                                                                                                                                                                                                   >  1:38011/1‑251 (MQ=255)                                      cGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAAt                                                                                                                                                                                                               >  1:649535/1‑251 (MQ=255)                                                         gTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGtt                                                                                                                                                                                             >  1:477831/1‑250 (MQ=255)                                                                         gAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCaa                                                                                                                                                                              <  2:149957/249‑1 (MQ=255)                                                                           aTGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAACCTg                                                                                                                                                                          >  1:351956/1‑251 (MQ=255)                                                                              cGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGAtt                                                                                                                                                                       >  2:200077/1‑251 (MQ=255)                                                                                 gCCCATTCCTGACGGAAAGGGCATTTCTGCGCAACTTTTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCCCGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATc                                                                                                                                                                    <  2:595696/251‑1 (MQ=255)                                                                                    cATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAGTCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGAt                                                                                                                                                                 <  2:219009/251‑1 (MQ=255)                                                                                    cATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGAt                                                                                                                                                                 <  1:542530/251‑1 (MQ=255)                                                                                     aTTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATc                                                                                                                                                                <  2:21895/251‑1 (MQ=255)                                                                                           gACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGa                                                                                                                                                            <  1:630209/249‑1 (MQ=255)                                                                                           gACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACg                                                                                                                                                          <  2:356490/251‑1 (MQ=255)                                                                                           gACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTGCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGAc                                                                                                                                                           <  2:200515/250‑1 (MQ=255)                                                                                                aaTGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAg                                                                                                                                                                                           >  1:288061/1‑213 (MQ=255)                                                                                                aaTGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAg                                                                                                                                                                                           <  2:288061/213‑1 (MQ=255)                                                                                                             gcgcAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTc                                                                                                                                         <  1:613049/250‑1 (MQ=255)                                                                                                                    ttgttgTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAAc                                                                                                                                 >  1:21862/1‑251 (MQ=255)                                                                                                                         gtCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAg                                                                                                                             <  1:228514/250‑1 (MQ=255)                                                                                                                                caacaaATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCaa                                                                                                                     <  1:613449/251‑1 (MQ=255)                                                                                                                                     aaTTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAgctgc                                                                                                                >  1:439251/1‑251 (MQ=255)                                                                                                                                       ttACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTCCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGctgctg                                                                                                              >  1:479404/1‑251 (MQ=255)                                                                                                                                        tACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGt                                                                                                             >  1:582916/1‑251 (MQ=255)                                                                                                                                                      ntCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCt                                                                                               >  2:325039/2‑251 (MQ=255)                                                                                                                                                        cAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCttt                                                                                             >  1:142489/1‑251 (MQ=255)                                                                                                                                                         aGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCttt                                                                                             >  2:613919/1‑250 (MQ=255)                                                                                                                                                             ataataTTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTc                                                                                        >  1:421041/1‑251 (MQ=255)                                                                                                                                                                     gATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCt                                                                                <  2:152229/251‑1 (MQ=255)                                                                                                                                                                            aGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTcc                                                                         >  2:474867/1‑251 (MQ=255)                                                                                                                                                                                   ttGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGtt                                                                  >  1:641532/1‑251 (MQ=255)                                                                                                                                                                                        atatCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAccc                                                             >  1:673337/1‑251 (MQ=255)                                                                                                                                                                                        atatCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTAccc                                                             <  2:430681/251‑1 (MQ=255)                                                                                                                                                                                               ttattattaATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTc                                                        >  1:690645/1‑249 (MQ=255)                                                                                                                                                                                               ttattattaATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATAACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTc                                                        >  1:414201/1‑249 (MQ=255)                                                                                                                                                                                                          ttttATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGTACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGTAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAATCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACt                                            <  1:76590/250‑1 (MQ=255)                                                                                                                                                                                                             tctGCGTTATACGGATAGCAGAATACTCACGCCTGCAACGGACGCCACTCACACGACCGCTTCCCCGCAGCCTATTCGACGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGAGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGAAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAACCGTTTACCCCCGTCAGGTCGACGACTggc                                         >  2:520361/3‑248 (MQ=255)                                                                                                                                                                                                                  nttATACGGATAGAAGAAAACTCAAGCGTGAAACGGACGCCAATCATAAGAAAGATTCACCGCAGCCTATTCGGAGAATTTACAGCAAGACACTGTTAGGTCAGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATATGCTCCGTAAAAATCAGGCTGGCAAGCTGCTGTTGACAAAGTAGCCTTTAACTAGAGCAGCTGACTTCAAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGt                                    >  2:311000/2‑250 (MQ=255)                                                                                                                                                                                                                  gTTATACGGATAGCAGAAACCTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGtc                                   <  2:649535/251‑1 (MQ=255)                                                                                                                                                                                                                  gTTATACGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGtc                                   >  1:21881/1‑251 (MQ=255)                                                                                                                                                                                                                       aCGGATAGCAGACAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTCTTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGCCCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTcatcat                              >  2:354676/1‑251 (MQ=255)                                                                                                                                                                                                                       aCGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTcatcat                              >  2:345963/1‑251 (MQ=255)                                                                                                                                                                                                                       aCGGATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTcatcat                              <  2:582916/251‑1 (MQ=255)                                                                                                                                                                                                                          gATAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCAtctc                           >  2:569741/1‑251 (MQ=255)                                                                                                                                                                                                                           aTAGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATctct                          >  1:77297/1‑251 (MQ=255)                                                                                                                                                                                                                             aGCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATT‑TCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTc                       >  1:173151/1‑251 (MQ=255)                                                                                                                                                                                                                              gCAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTc                       <  2:351956/251‑1 (MQ=255)                                                                                                                                                                                                                               cAGAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCg                      >  1:577160/1‑251 (MQ=255)                                                                                                                                                                                                                                 gAAAACTCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCt                    >  2:504508/1‑251 (MQ=255)                                                                                                                                                                                                                                       tCACGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAgg              >  2:588824/1‑251 (MQ=255)                                                                                                                                                                                                                                         ccGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGAt            <  2:477831/250‑1 (MQ=255)                                                                                                                                                                                                                                         aCGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGAt            <  1:569741/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGCCTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGAtt           <  1:86447/251‑1 (MQ=255)                                                                                                                                                                                                                                            ccTGAAACGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTc         <  1:200077/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aaCGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAAc    >  2:75618/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cGGACGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGc  >  1:310718/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ACCCGTGTCTGGGCTGGGAGGATACCGAACGCCTGGTCGAAAAACTCGCCTCTGCGGTAGATACCCGCTTCTGAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGCAAGAATATCTGCTCCGTAAAACTCAGGCTGGCAAGCTGCTGTTGACCAAGTAGCCTTTAACTCGAGCAGCTGACTTCCAGCCGTTTACCCCAGTCAGGTGGACGACTGACATAGTCATCATCTCTGTCGCTGAAAGGATTTCGCAACGC  >  NC_000913/1789408‑1789900

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