Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,893	A→C	A219A (GCT→GCG)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,893	0	A	C	100.0%	203.8 / NA	62	A219A (GCT→GCG)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (39/23);  total (39/23)

TTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTA  >  NC_000913/3815647‑3816137                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     ttAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGcc                                                                                                                                                                                                                                                   >  1:102553/1‑251 (MQ=255) ttAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGcc                                                                                                                                                                                                                                                   >  2:144854/1‑251 (MQ=255)  tAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGcc                                                                                                                                                                                                                                                   >  2:421522/1‑250 (MQ=255)            ggCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGCTTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAggtg                                                                                                                                                                                                                                          >  1:288956/1‑246 (MQ=255)             gCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACTCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                                       <  1:224317/251‑1 (MQ=255)                     gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGACAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTc                                                                                                                                                                                                                               >  1:395653/1‑251 (MQ=255)                          aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGCAGCTCTTCATggg                                                                                                                                                                                                                          >  2:349891/1‑251 (MQ=255)                           ggggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGGAGTAGCTCTTCATgggg                                                                                                                                                                                                                         >  1:175350/1‑250 (MQ=255)                             ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGa                                                                                                                                                                                                                       <  1:296712/251‑1 (MQ=255)                                ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTCAGCGCAACTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCCGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                                    >  1:417956/1‑251 (MQ=255)                                 tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAAc                                                                                                                                                                                                                     <  2:102553/249‑1 (MQ=255)                                      gcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGc                                                                                                                                                                                                              >  1:384951/1‑251 (MQ=255)                                         cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGNAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                                                                                                                                                           >  2:246107/1‑251 (MQ=255)                                          ggATTTCAGCGTAAACGCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                                                                                                                                                                                                          <  2:262865/251‑1 (MQ=255)                                           aaggtccggcgaaaCTCGACAAACTTTAACAGCCGCTCGCTGAGCGGAAAATCAATAATCTTGCGCTGATAGGGTGTCATGCCTTCGCTGCTTATCTTACCTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGGCTTAAAAATCTGTTTGCCAGCGCCGCCTGCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCGCGCGCCAGCGCCAACAAGGTGGTTAGCTCTTCATGGGTGGACGTCGCGCCtact                                                                                                                                                                                                         <  2:385157/240‑4 (MQ=255)                                            aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGCGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                                                        >  2:120919/1‑251 (MQ=255)                                            aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg                                                                                                                                                                                                        <  1:48753/251‑1 (MQ=255)                                                                aaCTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGTTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAatg                                                                                                                                                                                    <  2:175350/251‑1 (MQ=255)                                                                                tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                                                      <  2:186771/249‑1 (MQ=255)                                                                                tgctAAGCGCAAAGTCAATAAACTGGCGCTGATATGGTTTCATGCCGTCGCTCCTCATCTTACTTTTCAACAGCCAAAAAAAAGGCGACTAATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCTGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCGGCCGTCCCCTGCACTTCGATGATGCGCCCGTCTTCg                                                                                                                                                                     <  2:184425/250‑1 (MQ=255)                                                                                    aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                                                                                   <  1:445219/248‑1 (MQ=255)                                                                                       cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                                             <  1:349891/251‑1 (MQ=255)                                                                                       cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                                             <  1:421522/251‑1 (MQ=255)                                                                                        gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACt                                                                                                                                                            >  1:199131/1‑251 (MQ=255)                                                                                          aaaTTCACTAAACAGGCGCTGATATGGTTTCATGCCTTCGCTCCTAATCTTACTTTTCTACAGACAAAAACAAGGCGCATCCTCAGTCGCCTTAAAAATAAGTTTGCCAGCGCAGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAAACAGGTGAGAAGCCCTTAATGGGTGAACGGCCCGCCTTCAGCCGTCCCATGCACTCAAATGATGCCACCGTCTTCGGtcatccccac                                                                                                                                                          >  2:436232/1‑246 (MQ=255)                                                                                               cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                                                                                      >  1:407867/1‑250 (MQ=255)                                                                                                    aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAACAAAA‑GGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGAGTCGCTACACTGGATTCGA‑T‑T‑CCCCCTCGGCCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCATTCTGCCGTCCCCTGCACTTCAATGAAGCGCCCGTCTTACGTCATCCCTACGTTCATGGCgg                                                                                                                                               >  2:463087/1‑251 (MQ=255)                                                                                                             tGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                                        <  1:482581/250‑1 (MQ=255)                                                                                                              gATATAGTTTCTTGCCGTCTCTCCACATCTTGTTTTTCTACAGTCTAAAAAAAGGCGACTCCTCAGTCTCCTTAAAAATCCGTTTGCCCGCGCCGCCTGCTGTGTCCCTCCAATTCATTCGA‑GGC‑CCCC‑TCGCGCCAGCGCCAACTAGGTGAGTAGCTCTTCATGGGTGAACCGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctcttt                                                                                                                                        <  2:350526/249‑3 (MQ=255)                                                                                                                 atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGACGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTACTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCgtctg                                                                                                                                   >  2:289857/1‑247 (MQ=255)                                                                                                                 atGGTTACATGCCTTCGCTTCTCAGCTTCCCTTTCTCCAGACAAAAAAAAGGCGATTCATCAGTCGCTTTAAAAACCAGTTTGCCCGCGCCGCCTTCTGCGTCGCTAAAATGTATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCGTTATGGGTGAACGGCTCTCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGCATGCGGTCATCCCTACTCTCATGTCGGTCTCTTCGGct                                                                                                                                    <  2:381285/250‑2 (MQ=255)                                                                                                                           gccctCGCTCCTCATGTTACTTTTCTACAAACAAAAAAAAGGCGACACATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTTCAATGGATTCGC‑T‑T‑CCCCCTCGGGCCAGCGTCAACAAGGTGCTTAGTGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCaat                                                                                                                         <  2:54483/247‑2 (MQ=255)                                                                                                                                 gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTATCAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCCACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATTCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCTTGGTCTTCAACGTCTTc                                                                                                                   >  2:156377/1‑251 (MQ=255)                                                                                                                                    cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTcc                                                                                                                  >  1:366602/1‑249 (MQ=255)                                                                                                                                       cATCTGACTTCTCTAACGACAAAAAAAAGGCGACACAGCAGTCGCCTTAAAAATCACTTTGCCGGCGCCGCCTTCTGCGTCGCTACAATGGATTCGACT‑C‑CCCC‑GCTGGCCAGCGCCAACCAGTTGTGTAGCTCTTCATTGTTGAACGGCTCGCCCTCTGCCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCTGCGTCTGCGTCTTCAACGTATTCCAGagc                                                                                                             <  2:20196/251‑3 (MQ=255)                                                                                                                                            tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCaa                                                                                                         >  2:252826/1‑250 (MQ=255)                                                                                                                                             aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaa                                                                                                        <  1:137367/250‑1 (MQ=255)                                                                                                                                                tttCTACAGACAAAAAAAAGGCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCaaattgc                                                                                                    >  2:161382/1‑247 (MQ=255)                                                                                                                                                           aaataaaaGGCGACTAATCAGTCGCCTTTCACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑TTCACA‑TCAGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTCCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTTCGTTCATGTCGGTCTCGGCGGCGGCGTCTGCCACGTATCCCAGATCGCaatttgcttcaacgttt                                                                                          >  2:28341/5‑235 (MQ=255)                                                                                                                                                              aaaaaGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCTGTAATCACTACGTTCATGTCGGTCTCTGCGCCTGAGTCTTCAACGTATTCCAGCTCGctaatttcttcaccgttttct                                                                                       >  2:256339/1‑230 (MQ=255)                                                                                                                                                                 atggCGATTCATCAGTAACCTTAAAGATCAGTTTGACAGCCACGAATTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGCCGGTCTCTGCGGCTGAGTCTTCAACGTATTTCAGCTCGCAAATTGCTTCCCCGTTTACAAtac                                                                                   >  2:273046/3‑249 (MQ=255)                                                                                                                                                                    gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                 >  1:350878/1‑250 (MQ=255)                                                                                                                                                                    gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTAAGCAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCCTCACTACGTTCATTTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                 >  1:245359/1‑250 (MQ=255)                                                                                                                                                                          catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGa                                                                            <  2:377798/249‑1 (MQ=255)                                                                                                                                                                             caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTg                                                                       >  1:355915/1‑251 (MQ=255)                                                                                                                                                                                    cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGGGAGAAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTTCGGCTAc                                                                >  1:413564/1‑251 (MQ=255)                                                                                                                                                                                        aaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATc                                                            >  1:284774/1‑251 (MQ=255)                                                                                                                                                                                             cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCCCCGTTTACAATACCCACAGAAACTGCGGCTACCATCCCttt                                                       >  2:41603/1‑251 (MQ=255)                                                                                                                                                                                                   gccacccccgccTTCTGCGTCTCTACAATGGATTCGAGT‑C‑CCCC‑TCGGGCCAGCGCCAACAAGGTGGGGAGCTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGGCCGGCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCTGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTGCACCGTTTACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                                  <  1:273046/243‑1 (MQ=255)                                                                                                                                                                                                     cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                                                                                                                                                                                         <  1:344277/97‑1 (MQ=255)                                                                                                                                                                                                     cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                                                                                                                                                                                         >  2:344277/1‑97 (MQ=255)                                                                                                                                                                                                        cgccgcCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCACGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTgg                                            >  1:102560/1‑251 (MQ=255)                                                                                                                                                                                                         gccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                           >  1:451433/1‑251 (MQ=255)                                                                                                                                                                                                         gccgccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                           <  1:3097/251‑1 (MQ=255)                                                                                                                                                                                                         accgccTTCTGCGTCGATACAATGGATTCGA‑T‑T‑CCCCCTCGGTCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATTATGCGCCCATCTTCGGTCATCACTACGTTCATGTCGATCTCTTCGGCTGCGCCTTCAACGTATTCCCGCTCGCACATTGCTTCACCTTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCCTCggttttgt                                           >  2:63/2‑245 (MQ=255)                                                                                                                                                                                                            gccTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGtttt                                        >  1:58825/1‑251 (MQ=255)                                                                                                                                                                                                                tcggCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCGCGGGCCAGCGTCAACAAGGTGAGTAGTTCGTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAg                                     <  2:333823/247‑1 (MQ=255)                                                                                                                                                                                                                 cTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCACTTTCATCGGGTTGGTTTTAAGCt                                   >  2:484868/1‑251 (MQ=255)                                                                                                                                                                                                                   gCGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATAGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTg                                 >  2:54780/1‑251 (MQ=255)                                                                                                                                                                                                                    cGTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCCACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACTGCGGCTACCATCCCCTTCATCGGATTGGTGTTAAGCTTGc                                >  2:441607/1‑251 (MQ=255)                                                                                                                                                                                                                     gTCGCTACAATGGATTCGA‑T‑T‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTAAGCtttcc                               >  2:336014/1‑248 (MQ=255)                                                                                                                                                                                                                                                     ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGCGTCTTCAACGTATTCCAGATCGCCAATTGCTTCACCGTTTACAATACCGGCCGAAACTGCGGCTACCATCCCTTTCCTCGGGTTGGTTTTTAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtg  >  2:3875/1‑250 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     TTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑T‑T‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTA  >  NC_000913/3815647‑3816137

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