Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,162	T→A	S417S (TCT→TCA)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,162	0	T	A	100.0%	172.5 / NA	53	S417S (TCT→TCA)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (23/30);  total (23/30)

GGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAA  >  NC_000913/4626925‑4627400                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                ggaaattctgctgTCTGTCTGGCACCGGCACGGGGGGCTGCAAATGGCCGATCAGGATGTGTTGGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCGTTGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCACAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAac                                                                                                                                                                                                                                    <  2:272108/244‑1 (MQ=255)    aagccTGCTGGCTTTCTTGCCCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑g                                                                                                                                                                                                                                <  2:386556/248‑1 (MQ=255)           cgggctttcttGCACCGTCTCGGTGTTCTGCAAATGGCCGATCTGGATGTGGTTGTTAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTGCCAGCCGGCGTGACATACCGCTGCCATAGGAGCTGGTGGTGTTTGGTGAAGTCTGGCTGGCAGGGGAGATCCGCCGGTTGCCCAGCGTTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCg                                                                                                                                                                                                                         <  1:384155/242‑1 (MQ=255)                 gtctttcccgggagtgggggtCTGCAAATGGCCGACCTGGATTTGTTTGTGAACGTGGTCGGGGGCGGGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCGGCGTTAGAGACCGCTGCCTCAGGATCTGGTGGTGTTTGGTGAAGCGGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACTAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCggg                                                                                                                                                                                                                    <  2:505363/232‑1 (MQ=255)                       cACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGTGTGAAGGTAACCGAAACCAGCGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGt                                                                                                                                                                                                             <  1:273977/251‑1 (MQ=255)                         ccGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTc                                                                                                                                                                                                           >  1:609796/1‑251 (MQ=255)                             cacggtgggctgcgagcggcggagcAGGATGTGGTGGTGATCGGGGTCGGTGGCGTGAGGGCAGCCGAAACCAGGGCCGACTTAGCGTTATTGCTGGCGTTGGTTTCCAGCCTGCGTGAGAGCCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCTGTGCCCAGCGGTCATGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGTCGGGCGATTGTTCcgg                                                                                                                                                                                                        <  2:265148/226‑1 (MQ=255)                             cACGGTAGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGAAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGCCAGGAACGAATTTCAGAAGCGGCGGAACAA‑GGATTTCGCCGGGCGTTTGTTCcggc                                                                                                                                                                                                       >  1:632105/1‑251 (MQ=255)                             aaCGGGGGTCTGAAAATGGCCTACCAGGATGTGTTGGGGAAAGTGGTCGGTGGCGTGAAGGTAACCGAATCCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCcgg                                                                                                                                                                                                        <  2:107367/249‑1 (MQ=255)                               cGGTGGTCGGCAAATGACCGATCAGGATGTGTTTGGGAACGTGGTCGGCGGCGTGAAGGAAACCGACACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGAGCACAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAA‑ACAGGTTTTCGACGGGAGCTTGTTCcggcgg                                                                                                                                                                                                     >  2:525456/1‑251 (MQ=255)                                        gaaaggggCCGATCGGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTGCCAGCCGGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGtac                                                                                                                                                                                            <  2:81565/245‑3 (MQ=255)                                               gCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGAGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTACAGTACCGaaaa                                                                                                                                                                                      >  2:491896/1‑250 (MQ=255)                                                     cAGGATGTGTTGTTGAAAGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGGGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgcgc                                                                                                                                                                                <  1:528094/250‑1 (MQ=255)                                                           gtgtTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAgg                                                                                                                                                                         <  1:236611/251‑1 (MQ=255)                                                                tgtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATg                                                                                                                                                                     <  1:491896/250‑1 (MQ=255)                                                                 gtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAACAAGCGCCGGAAGGGATGc                                                                                                                                                                    >  2:232563/1‑250 (MQ=255)                                                                                   gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcgg                                                                                                                                                                                                                                          >  1:463124/1‑163 (MQ=255)                                                                                     gAAGGGAACCGAAACCAGCCCCGCCTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGTCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAc                                                                                                                                               <  2:83922/251‑1 (MQ=255)                                                                                      aaGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAct                                                                                                                                              >  1:149109/1‑251 (MQ=255)                                                                                       aGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCt                                                                                                                                                                                                   >  1:315603/1‑197 (MQ=255)                                                                                           aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTcc                                                                                                                                          <  1:285628/250‑1 (MQ=255)                                                                                           aaCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCAGGCGGCTAACGTAACGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTCAAAAACTCTCCg                                                                                                                                         >  1:73171/1‑251 (MQ=255)                                                                                                   ccAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt                                                                                                                                 <  1:523349/251‑1 (MQ=255)                                                                                                   ccAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCtt                                                                                                                                 >  1:576170/1‑251 (MQ=255)                                                                                                    cagcgCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTa                                                                                                                                <  1:446095/247‑1 (MQ=255)                                                                                                     agcgCCGACTTAGCGTTACTGCTGGAGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGGTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGGTTGTTCCGGCGGCTAACGTACCGAAAAAAGCGACGGAAGGGATGAAGCATTTTGGCGTTAAAAAGctcctcaccccgctt                                                                                                                                 >  2:163077/4‑237 (MQ=255)                                                                                                           gACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGGGtt                                                                                                                         >  1:166709/1‑251 (MQ=255)                                                                                                           gACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAA‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAGAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAAGCGCTTAGCGTGtt                                                                                                                         >  1:521929/1‑251 (MQ=255)                                                                                                            agttAGCGTTACTGCTGGCGAGGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGTGCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCGCAGCGGTCAGGGACGAAGTGCAGAAGCGTCGAAACAC‑GGTTTTCGCCGGGCGATGGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTGTGGCGTTAAAAAACTCTCCGACGCGGGTAGTGTGTTc                                                                                                                        <  2:79731/249‑1 (MQ=255)                                                                                                                 ggtgtACTGCTGGCGGGGGTTTCCAGCCTGCGTGACAGATCGCTGCCACAGTATCTGGTGGTGTTTGGTGAAGGCGGGCTGGCAGGGGAGATCCTCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTGTGCGTGTTcgacg                                                                                                                    <  2:290367/246‑1 (MQ=255)                                                                                                                    ttACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACtt                                                                                                                >  1:303629/1‑251 (MQ=255)                                                                                                                         gctgggtttGGTTTCCAGCCTTCGTGAAAGATCTCGGCCACAGGATCTGGGGGTGTTTGGGGAAGTCGGGCTGGCAGTGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATGTCAGAAGCGGCGAAACAC‑GGGTTTCGCCTGGCGATTGTTCCGGTGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGTCGCGCTTAGCGTGTTCGACGATTTATAAt                                                                                                           <  2:539755/243‑1 (MQ=255)                                                                                                                         gctgGCGATGGTGTCCAGCCGGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATaa                                                                                                            <  1:609692/250‑1 (MQ=255)                                                                                                                          ctgGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATg                                                                                                          >  1:636584/1‑251 (MQ=255)                                                                                                                             gCGATGGTTTCAAGCCTGCTGGACAGACCGCTGCCACAGGAGCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTAGAATgagg                                                                                                       <  2:609796/251‑2 (MQ=255)                                                                                                                                      tCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACgga                                                                                               <  1:525456/250‑1 (MQ=255)                                                                                                                                                   aCAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGGTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATgtcgt                                                                                 <  2:166709/251‑1 (MQ=255)                                                                                                                                                               cacAGGATCTGGTGGTGTTTGGAGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGGCTTATAATGAGGAATACGGGGGGAGATATGTCGTCTTTTGATTAcc                                                                     >  2:108868/3‑251 (MQ=255)                                                                                                                                                                acAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATgag                                                                                                        <  1:159622/214‑1 (MQ=255)                                                                                                                                                                         ggtggtGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAAc                                                              <  1:564928/248‑1 (MQ=255)                                                                                                                                                                          gtggtgTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGcc                                                          <  1:123332/251‑1 (MQ=255)                                                                                                                                                                               gtTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCa                                                      <  1:33711/250‑1 (MQ=255)                                                                                                                                                                               gtTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCa                                                      >  1:673933/1‑250 (MQ=255)                                                                                                                                                                                  tGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTCGCGTGTTCGACGACtttt                                                                                                              >  2:23174/1‑189 (MQ=255)                                                                                                                                                                                  tGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTtat                                                                                                              <  1:23174/191‑1 (MQ=255)                                                                                                                                                                                            gggCTGGCAGTGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGc                                        >  1:541912/1‑251 (MQ=255)                                                                                                                                                                                                   cAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTAc                                 >  1:415020/1‑251 (MQ=255)                                                                                                                                                                                                      gggAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTAcagc                              >  1:622946/1‑251 (MQ=255)                                                                                                                                                                                                                          aGCGGTCAGGAACGAATGTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGt          <  1:499498/251‑1 (MQ=255)                                                                                                                                                                                                                             ggTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACa                                               <  2:34286/211‑1 (MQ=255)                                                                                                                                                                                                                             ggTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACa                                               >  1:34286/1‑211 (MQ=255)                                                                                                                                                                                                                              gTCAGGAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATg       >  1:227851/1‑250 (MQ=255)                                                                                                                                                                                                                                  ggAACGAATTTCAGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCaa  <  1:465289/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                GGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACAC‑GGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAA  >  NC_000913/4626925‑4627400

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