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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,166	0	C	A	100.0%	228.4 / NA	69	A128A (GCG→GCT)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (42/27);  total (42/27)

GGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCA  >  NC_000913/3815920‑3816410                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     ggTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCgg                                                                                                                                                                                                                                                  >  2:354385/1‑251 (MQ=255)                                 cTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCAAAGCGGTACGCGTGCCACCATCAGCCTGACGCACGTCg                                                                                                                                                                                                                  >  2:207745/1‑251 (MQ=255)                                        gatgCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCAGCAGCCTGAAGCACGTCGCAGTCCa                                                                                                                                                                                                           <  2:438059/251‑1 (MQ=255)                                          tgCGCCCTTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCTTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAg                                                                                                                                                                                                          <  1:29982/250‑1 (MQ=255)                                             gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                                                                                                                                                                                       >  2:90046/1‑250 (MQ=255)                                             gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTa                                                                                                                                                                                                      <  1:180941/251‑1 (MQ=255)                                                   cTTCGGTCATCCCTAAGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTAACATCCCTTTCATTGGATTGGTTTTCCGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCCCCATCAGCCCGACGCACGTCCCCGTCCAGCGTAATGGTg                                                                                                                                                                                                >  1:359794/1‑251 (MQ=255)                                                      cGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAAc                                                                                                                                                                                             <  1:32142/251‑1 (MQ=255)                                                            tcaCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGCTCGCAAATTGCTTCACCGTTTACAATAACGACAGAAACCGCGGCTACCATCCCTTTCCTTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                                                                                                                                                       >  1:431867/1‑251 (MQ=255)                                                              aCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                                                                                                                                                       >  1:98349/1‑249 (MQ=255)                                                                   gTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCTCGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCtt                                                                                                                                                                                <  2:39345/251‑1 (MQ=255)                                                                           cGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCTACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTTCCACCATCAGCCTGAAGCACGTCGCAGGCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATc                                                                                                                                                                         <  2:261527/250‑1 (MQ=255)                                                                                   gCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACGCACCCAGCGCTTTCAAATCTACTGCcg                                                                                                                                                                 <  2:43264/250‑1 (MQ=255)                                                                                        tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTCATGGTGAACTCACCCAGCGCTTTCAAATCTCCTGCCGCGCGa                                                                                                                                                            >  2:75994/2‑250 (MQ=255)                                                                                        tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGcca                                                                                                                                                                                                                                       >  1:41135/2‑175 (MQ=255)                                                                                        tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGcca                                                                                                                                                                                                                                       <  2:41135/174‑1 (MQ=255)                                                                                        tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAACGCTTTCAAATCTACTGCCGCGCGa                                                                                                                                                            >  1:261153/2‑250 (MQ=255)                                                                                           gTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAagtg                                                                                                                                                        >  2:85361/1‑249 (MQ=255)                                                                                                  aCGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGCAGTGCACGGGc                                                                                                                                                 >  2:200923/1‑251 (MQ=255)                                                                                                     tATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                                                                              >  1:243691/1‑251 (MQ=255)                                                                                                           aGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCCGCGGCTACCATCCCTTTCATTGGATTGGTTGTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGatg                                                                                                                                        <  2:302192/251‑3 (MQ=255)                                                                                                             aTCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCt                                                                                                                                      <  1:396126/251‑1 (MQ=255)                                                                                                                       gCTTCACCGTTTACAATACAGACAGAAACCCCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCAGTTTTCCAACAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCAt                                                                                                                            >  1:159969/1‑251 (MQ=255)                                                                                                                        cTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCAtt                                                                                                                           >  1:151756/1‑251 (MQ=255)                                                                                                                           caccGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATtgtg                                                                                                                        <  1:184819/249‑1 (MQ=255)                                                                                                                            accGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTATAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGAGGTACGCGTGCCACCATCGGCCTGACGCACGTCGCAGTCCAGAGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCCAAGTGCACGGTCGAACAGCCGCTGGTTTTCCAttacgc                                                                                                                       >  1:430867/2‑247 (MQ=255)                                                                                                                                  tacaATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCCGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAc                                                                                                                 >  2:352252/2‑251 (MQ=255)                                                                                                                                   acaATACCGACAGAAACCGCGGTTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCATACGCTGGATTTCCATTGTGCGTCCAcc                                                                                                                >  1:139371/1‑251 (MQ=255)                                                                                                                                   acaATACCGACAGAAACCGAGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGCGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAcc                                                                                                                >  2:49218/1‑251 (MQ=255)                                                                                                                                    caATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCCCGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAccc                                                                                                               >  2:193784/1‑251 (MQ=255)                                                                                                                                      ataccGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGGTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTg                                                                                                             <  2:90908/248‑1 (MQ=255)                                                                                                                                           gACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTAAGCTTGCCGTTTTCCACCAGCTTATGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAAATCACACAGCGCTTTCAAATCTACTGCCGCGAGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCttc                                                                                                         >  1:279953/1‑249 (MQ=255)                                                                                                                                            aCAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTAcc                                                                                                       >  1:51288/1‑251 (MQ=255)                                                                                                                                            aCAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTAcc                                                                                                       <  1:85361/251‑1 (MQ=255)                                                                                                                                              aGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCTCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCtt                                                                                                     <  2:167034/251‑1 (MQ=255)                                                                                                                                               gAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCttt                                                                                                    >  2:301654/1‑251 (MQ=255)                                                                                                                                                       ggCTACCATCCCTTTCATTGGATGGGTTTTCAGCTTTCCGTTTTCCACCAGCATCTGTAGCCCATATGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCTAGCGTAATGGAGAACTCACCCAGCGCTTTCAAATCTACCGCTGCACGAGGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTACACCCTGCTTACCTTTCGCTGCtt                                                                                            >  2:130871/1‑251 (MQ=255)                                                                                                                                                        gCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCACGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTc                                                                                           >  1:399401/1‑251 (MQ=255)                                                                                                                                                           aCCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTTTGTAGCGCATCTGCCAGCGCCACTCAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGCAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACg                                                                                        >  2:59546/1‑251 (MQ=255)                                                                                                                                                           aCCATCCCTTTCATTGGATTGGTGTTCAGCGTGCCGTTTTCCTCCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGGGCCACCATCAGCCTGAAGCACGTCGCAGTCCGGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCCCGAAGGGCACGTGCGATCAGACGCTGGATGGCCATTGTGCGTCCACCCTGCTTTCCTTTCGCTGCATCACg                                                                                        <  2:431867/251‑1 (MQ=255)                                                                                                                                                             cATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACg                                                                                        >  2:134325/1‑249 (MQ=255)                                                                                                                                                              aTCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGc                                                                                     >  1:208247/1‑251 (MQ=255)                                                                                                                                                                    atcattgcaTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACCCGCGCCACCATAAGCCTGAAGCACGTCGCAGTCCCGCGTATTGGTGAACTCAACCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCGTTCTGCGTCCACACTGCTTACCTTCCGCTGCTCCACGCGCGTTACg                                                                               >  2:275140/9‑251 (MQ=255)                                                                                                                                                                      cattggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACgg                                                                              <  1:356388/246‑1 (MQ=255)                                                                                                                                                                      cattggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGCATCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACTTCGCAGTCCCGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCTCACGTAGTGCACGGGCGCTCAGACGCTGGCTTTCCATTGGGCGTCCACCCTGCTTACCTTTCGCTGCTTCCCGATCGTTACggg                                                                             >  2:232991/5‑251 (MQ=255)                                                                                                                                                                       attggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGAATCTGCCAGCGCCAAACATGCACACGTAATCGAAGCAGTACGCGTGCCACCATCCGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACATTTCGCTGCTACACGAGCGTTACGGgt                                                                            >  2:182279/4‑251 (MQ=255)                                                                                                                                                                          ggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGtgtg                                                                         <  1:67286/251‑1 (MQ=255)                                                                                                                                                                           gATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTgg                                                                        >  1:289358/1‑251 (MQ=255)                                                                                                                                                                             ttGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGt                                                                      >  1:42037/1‑251 (MQ=255)                                                                                                                                                                               ggTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCCGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGt                                                                      >  2:278577/1‑249 (MQ=255)                                                                                                                                                                                   ttCAGCTTGCCGTTTTCCACCAGCTGCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTACTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAAc                                                                 <  2:419077/250‑1 (MQ=255)                                                                                                                                                                                     cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                                                                              <  2:188296/171‑1 (MQ=255)                                                                                                                                                                                     cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                                                                              >  1:188296/1‑171 (MQ=255)                                                                                                                                                                                     cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACgcg                                                              >  1:176700/1‑249 (MQ=255)                                                                                                                                                                                      aGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACgcgg                                                             >  1:127309/1‑248 (MQ=255)                                                                                                                                                                                       gcgtGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCTGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTTAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCGCCCTGCTTACCGTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGatcgcggc                                                            <  2:210769/248‑8 (MQ=255)                                                                                                                                                                                          tGCCGTTTTCCACCATCTTCGGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAcc                                                                                                                                                                                        <  1:261137/124‑1 (MQ=255)                                                                                                                                                                                          tGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGgcagc                                                         >  1:190173/1‑251 (MQ=255)                                                                                                                                                                                          tGCCGTTTTCCACCAGATTCTGTAGCGCATCTGCCAGCGCCACGACGGCACCCGTCATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAcc                                                                                                                                                                                        >  2:261137/1‑124 (MQ=255)                                                                                                                                                                                                                                gcCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGAcc                   >  2:428132/1‑251 (MQ=255)                                                                                                                                                                                                                                  cacacAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCtt                 >  1:79537/5‑251 (MQ=255)                                                                                                                                                                                                                                   acacAGGCACCCGTAATCGAAGCGGTACGCGGTCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCGGGCCCTGACCttt                <  2:208247/248‑1 (MQ=255)                                                                                                                                                                                                                                       aGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAgg            >  1:140711/1‑251 (MQ=255)                                                                                                                                                                                                                                        ggCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGa           <  1:354385/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGaa          <  1:360591/251‑1 (MQ=255)                                                                                                                                                                                                                                            cccGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCATTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCg                                          <  2:251066/216‑1 (MQ=255)                                                                                                                                                                                                                                            cccGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCg                                          >  1:251066/1‑216 (MQ=255)                                                                                                                                                                                                                                            cccGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAgcg       <  2:151756/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aaTCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCa  <  1:86489/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCA  >  NC_000913/3815920‑3816410

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