Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,316	C→T	A78A (GCG→GCA)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,316	0	C	T	100.0%	233.0 / NA	68	A78A (GCG→GCA)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/37);  total (31/37)

GGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAA  >  NC_000913/3816070‑3816556                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 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)                               cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACgcg                                                                                                                                                                                                                >  1:176700/1‑249 (MQ=255)                                aGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACgcgg                                                                                                                                                                                                               >  1:127309/1‑248 (MQ=255)                                 gcgtGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCTGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTTAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCGCCCTGCTTACCGTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGatcgcggc                                                                                                                                                                                                              <  2:210769/248‑8 (MQ=255)                                    tGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGgcagc                                                                                                                                                                                                           >  1:190173/1‑251 (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)                                                                                               gaagcagTACGCGTGCCACCATCCGCCTGACGCACGTCGCAGTCCAGCGTACTGGTGCACTCACCCAGCGCTTTCAAATATACTGCCGCGCGAATTGCACGGGCGATCAGACGCTGGATTTACCTTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGAGGCAGCATACCGGCCTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGCAGCGAGGCCCGcc                                                                                                                                                >  2:186850/7‑251 (MQ=255)                                                                                                   cgGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGCGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCcttct                                                                                                                                            >  1:29974/1‑251 (MQ=255)                                                                                                             gCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATagag                                                                                                                                   <  1:292008/250‑1 (MQ=255)                                                                                                               caccaTCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCg                                                                                                                                <  1:428132/251‑1 (MQ=255)                                                                                                                           gAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACacaatac                                                                                                                     <  2:254615/250‑4 (MQ=255)                                                                                                                                    cGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCg                                                                                                           <  1:403903/251‑1 (MQ=255)                                                                                                                                       aGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGcc                                                                                                         >  1:306125/1‑250 (MQ=255)                                                                                                                                           cAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAAACAATACTTTGTTATCGCCAAAtt                                                                                                    >  2:240380/1‑251 (MQ=255)                                                                                                                                               gTAATGGTGAAATCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTATTCAATAGAGGCGGTACACAATACTGTGGGATCGCCAAATTctac                                                                                                >  1:88851/1‑248 (MQ=255)                                                                                                                                                 aaTGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGAcc                                                                                               <  1:200923/250‑1 (MQ=255)                                                                                                                                                     gTGACCTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAAACAATACTTTGGTATCGCCAAAttttcccagccc                                                                                          >  2:35267/1‑241 (MQ=255)                                                                                                                                                         aCTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGTCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTGTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAg                                                                                      <  2:136764/251‑1 (MQ=255)                                                                                                                                                                       ttCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtt                                                                         <  1:222902/250‑1 (MQ=255)                                                                                                                                                                         cAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATGTCCATTGTGCGGCCACCCGGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtttt                                                                       <  2:139371/250‑1 (MQ=255)                                                                                                                                                                                aCTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTTCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTCCTGCATGTTTTGTATAg                                                                 <  2:325160/249‑1 (MQ=255)                                                                                                                                                                                aCTGCCGAGCGAAGTGCACGGGCGACCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACAGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATTCTTTGGTATCGCCAAATTCGACCATCACCGAGCCTTCTGCATGTTTTGTATAg                                                                 <  2:7044/249‑1 (MQ=255)                                                                                                                                                                                  tGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACCCAATACTTTGGTATCGCCAAATTCGACCCACCCCGCGCCTTCTGCCTGTTTTTTATAGttcc                                                             >  2:30754/1‑249 (MQ=255)                                                                                                                                                                                   gCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTGCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGTTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTGGTATAGTTAc                                                             <  2:164877/250‑1 (MQ=255)                                                                                                                                                                                          gAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAACCAATACTTTGTTATCTCCAAAATCGACCCGCACCCAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAgg                                                     >  2:80745/1‑251 (MQ=255)                                                                                                                                                                                                cACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCTCCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACg                                               <  2:243691/251‑1 (MQ=255)                                                                                                                                                                                                 aCGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTTGACCCGCACCGAGCCCTCTGCATGTTTTGGATTGTTTCGAGTAAGGGTAACgg                                              >  2:120068/1‑251 (MQ=255)                                                                                                                                                                                                  cGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACggg                                             >  1:299716/1‑251 (MQ=255)                                                                                                                                                                                                             aCGCTGGATTTCCATTGTGCGTCCACCCTGCTTCCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGat                                  <  2:289358/251‑1 (MQ=255)                                                                                                                                                                                                             aCGCTGGATTTCCATTGTGCGTCCACCCTGCTGACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGat                                  <  1:402223/251‑1 (MQ=255)                                                                                                                                                                                                                     tttCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCATGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACGGGACGCACCTGATTATTGCt                           >  1:182864/1‑250 (MQ=255)                                                                                                                                                                                                                      ttCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTGTGTATAGTTACGAGTCAGGGTAACGGGACGCGCCTGATTAGTGCTa                          <  2:190173/250‑1 (MQ=255)                                                                                                                                                                                                                      ttCCATTGTGCGGCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACGCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCATCACCGAGCCTTCTTCATGTTTTGTATAGTTACGAGTCAGGGTAACGGTACGCATCTGATGAGTGCTa                          <  2:399401/250‑1 (MQ=255)                                                                                                                                                                                                                       tCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAggg                                                    <  1:133278/223‑1 (MQ=255)                                                                                                                                                                                                                        ccATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACgg                       <  1:90046/251‑1 (MQ=255)                                                                                                                                                                                                                            tgtgCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCct                    <  1:49218/250‑1 (MQ=255)                                                                                                                                                                                                                            tgtgCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCct                    <  1:301654/250‑1 (MQ=255)                                                                                                                                                                                                                            tgtgCGTCCACCCGGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCGACGGCct                    <  2:42037/250‑1 (MQ=255)                                                                                                                                                                                                                                 gTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGGATCGCCAAATTCGGCCAGCACCGGGCCTTCTGCATGTTTTGTATAGTTACGAGTAAGGGTAAAGGGAAGCACCTGATTATTGCTACGGCCTGCTggg              >  1:129834/1‑250 (MQ=255)                                                                                                                                                                                                                                   ccACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGGATCGCCAAATTTGACCAGCACCCAGCCTTCTGCCTGTTTTGTATAGTTACGGGTCCGGGTAAAGGGCCGCACCTGATTATTGCTACGGCctgcctggcg            >  1:17819/1‑245 (MQ=255)                                                                                                                                                                                                                                       ccTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCata        <  2:176700/251‑1 (MQ=255)                                                                                                                                                                                                                                       ccTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCata        <  1:75994/251‑1 (MQ=255)                                                                                                                                                                                                                                          gCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCTCCTGATTATTGCTACGGCCTGCTGGACGCATAtt      <  1:372273/250‑1 (MQ=255)                                                                                                                                                                                                                                          gCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCatat       <  1:35267/249‑1 (MQ=255)                                                                                                                                                                                                                                              aCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGaaa  <  1:134325/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 GGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAA  >  NC_000913/3816070‑3816556

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