Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,487	T→C	G192G (GGT→GGC)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,487	0	T	C	100.0%	142.7 / NA	48	G192G (GGT→GGC)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/23);  total (25/23)

AAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGA  >  NC_000913/4626252‑4626732                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        aaaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCCTCGCCCGGCGGCGTGGCGCGGGTg                                                                                                                                                                                                                                          >  2:148646/1‑251 (MQ=255)            aaaCTGGCCCAGCAGATGAAAACGATGTATGTCACCGGCGAAGAGTCGCTGCAACAGGAGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGGGCGTGGAAcggc                                                                                                                                                                                                                               >  1:392566/1‑251 (MQ=255)                        cagATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGac                                                                                                                                                                                                                   <  1:229379/251‑1 (MQ=255)                         agATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGGAACGGCGGCttttcttacg                                                                                                                                                                                                                  >  2:40660/1‑243 (MQ=255)                         agATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGacg                                                                                                                                                                                                                  >  1:331747/1‑250 (MQ=255)                              aaaaCGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCtt                                                                                                                                                                                                             <  1:177822/251‑1 (MQ=255)                                         tGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGCCGCGCCTTGCCAAAAcgc                                                                                                                                                                                                  >  2:357774/1‑251 (MQ=255)                                             cccGGAGAAGAGGCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAACCCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCTCgg                                                                                                                                                                                              <  2:173097/250‑1 (MQ=255)                                                ggCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCGGACGCGCTTCTCCAAAACGCGCGgtgt                                                                                                                                                                                           <  1:103090/251‑1 (MQ=255)                                                  cGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGCTGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCCGCGTGGCGCAGGTGCGGGAAACGGCGGCTTATCTGACGCGCCTCGCCCAAACGCGCGGGGTgg                                                                                                                                                                                         >  2:111345/1‑251 (MQ=255)                                                  cGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAATAATCGCCTGGCAGCGGTGCGCCGGTGCCTGAAACGGCGGCTTTTCTGACCCGCTTCGcccacacagccggtgtgg                                                                                                                                                                                         >  2:9771/1‑235 (MQ=255)                                                       agTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGAtt                                                                                                                                                                                    <  1:176943/251‑1 (MQ=255)                                                            cTGCAACAGGTGGCGATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGGGCCGGTGCGTGAAACGGAGGCTTATCTGGCGCCCTTCGCCCAAACGCGCGGCGTgggggttttcat                                                                                                                                                                               >  2:120539/1‑241 (MQ=255)                                                                    ggTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAAGTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggttggac                                                                                                                                                                       <  1:191822/251‑6 (MQ=255)                                                                         cAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCCCTTCGCCAAAACGCGGGGGGTGGCGGTTGTCTTGGTTGGGCACGtt                                                                                                                                                                  >  2:404747/1‑250 (MQ=255)                                                                           aTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTATCATGGTTGGACACGTAAc                                                                                                                                                                <  1:321269/251‑1 (MQ=255)                                                                              cgcgcTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCa                                                                                                                                                              >  1:129521/1‑250 (MQ=255)                                                                                       cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGCGCAGATCTGCCTGATTGCCGAAGACGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTAAAGTCATCGCCTGGCAGCGTGGCGCAGGTGAGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAAAGCGCGGTTGGGCGATTGTCATGGGTAGACACGTAACCCAAGAttgctc                                                                                                                                                    >  1:460729/1‑246 (MQ=255)                                                                                                   ccGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGGTGGACACGTAACCAAAGATGGCTCGCTGGCGGGccc                                                                                                                                        >  1:359096/1‑251 (MQ=255)                                                                                                                 cccgtccgatatggAAACCAGCATCGTGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGAGGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAac                                                                                                                          <  2:200376/239‑1 (MQ=255)                                                                                                                          gtCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATc                                                                                                                 <  1:111345/250‑1 (MQ=255)                                                                                                                             ggAAACCAGCATCGCGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTTAAGGGCTGGAacagaacgggag                                                                                                                >  1:200376/1‑240 (MQ=255)                                                                                                                             ggAAACCAGCATCGAGCACATCTGCCTGTGTGCAGAAGAAGAGCAACCGAAGCTGAGGGTAATTGACGCGACCCAGTGGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTTAAACTGCGGCTTATCTTACGCGCTTCCCCAAAACGCGCGGTGTGGCGAGTGTCATGGTTGGACTCGTAACCAAAGATGGCTCGCTGGCTGGCCCTAAAGTGCTGGAACACTG‑TATCgcct                                                                                                             <  2:425107/251‑4 (MQ=255)                                                                                                                                     gCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTg                                                                                                     <  1:338391/251‑1 (MQ=255)                                                                                                                                       aTCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCAAATGGCGGATGGACAGTCATAGCCTGGCAGCGTGGCGCAGGGGCGGGAAACGGAGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGAGGGGCGATTGTCATGGGTGGAAACGTAACCAAAGAGGGCTCGCTGGCGGGCCCTAAAGAGCTGGAAAACTG‑TATCGACAGTTCGGGGCt                                                                                                   >  1:176099/1‑251 (MQ=255)                                                                                                                                               gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑Ta                                                                                                                   <  2:235181/227‑1 (MQ=255)                                                                                                                                               gATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑Ta                                                                                                                   >  1:235181/1‑227 (MQ=255)                                                                                                                                                    gCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGCCTGTTCGGTGCTTTTGGA‑TGGCg                                                                                       >  1:76635/1‑249 (MQ=255)                                                                                                                                                       agATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGcc                                                                                  <  1:40660/250‑1 (MQ=255)                                                                                                                                                         aTTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TAACGCATGTTCGGTGCTTTTGGA‑TGGCGATGCCgg                                                                                >  1:347776/1‑250 (MQ=255)                                                                                                                                                                    agagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTc                                                                     <  1:148646/251‑1 (MQ=255)                                                                                                                                                                     gagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCg                                                                    <  1:447822/251‑1 (MQ=255)                                                                                                                                                                       ccaccgtaaGCTGATGGTAATTGACGCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTACCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCa                                                                  <  1:127369/244‑1 (MQ=255)                                                                                                                                                                           ccGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGGTTTTGGA‑TGGCGATGCCGAATCCCGGTTTCGCActtt                                                              >  1:423422/1‑248 (MQ=255)                                                                                                                                                                             gAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTTGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTGCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCActttac                                                            <  1:104364/251‑6 (MQ=255)                                                                                                                                                                               aGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCCGGCGGGGCCCAAAATGCTGGAAAACTG‑TATCGCCTTTTTGGTGCTTTTTGG‑TGGCGATGCCGACTCCCCTTTTCGCActtttcgc                                                          >  2:144475/1‑244 (MQ=255)                                                                                                                                                                                   gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGcc                                                      <  1:380788/251‑1 (MQ=255)                                                                                                                                                                                        tAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGGTTGTCATGGTTGGCCACGTAACCAAAGATGGCTCGCGGGCGGGCCCTAAAGTGCTGGAACACTG‑TATAGACTGTTTGGGGGGTTTGGG‑TGGGGATGCCGACTCCCCtttttgccctttaccccccccaca                                                  >  2:9160/1‑230 (MQ=255)                                                                                                                                                                                             gACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAcc                                             <  2:305367/250‑1 (MQ=255)                                                                                                                                                                                                         gTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGt                                <  1:298075/251‑1 (MQ=255)                                                                                                                                                                                                         gTGATGCATATGGCGGATGTACAGACATCGCCTGGCAGCGGGGAGCAGGTGCGGGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCCCGCTGGCGGGCCCTTAAGTGCTGGGACACTGGTA‑CGCATGTTCGGTGCCTTTGGAAGGGCGAT‑CAGCCTCCCGTTTTCCCAACTTACGCCGCCATAAAAACCCCttaggggaggg                                >  2:396544/1‑242 (MQ=255)                                                                                                                                                                                                            aTGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGaa                             <  1:35723/251‑1 (MQ=255)                                                                                                                                                                                                              gCATAGGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATg                           <  2:372209/251‑1 (MQ=255)                                                                                                                                                                                                                atatGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACCCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑GATCGACTGTTCGGGGCTATTGAA‑TGGCGTTGCCGACTCCCGGTTTCGCACTTTACGCAGCCATAAAAACCGCTTCAGGGCGGTGAATg                           >  1:457861/1‑249 (MQ=255)                                                                                                                                                                                                                  atGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTTGGAAACGTAACCAAAGATGGCTCGCTGGGGGGCCCTAAAGAGGGGGAGCACTG‑TTTCGACTGGTCGGGGCGTTTGGA‑TGGGGAGGCCGGCTCCCGTTTTTGCAATTTACGCCGCCCTAAAAACCCCTTCgggccgggggatgggat                       >  2:292824/1‑236 (MQ=255)                                                                                                                                                                                                                      cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACCCTG‑GATCGACTGTTCGGTGCTTTTGGG‑GGGCGGTGCCGACCCCCGTTTTCGCACTTTACGCAGCCATAAAAACCggctcgggcgggggaatgggctgggg                   >  2:45882/1‑226 (MQ=255)                                                                                                                                                                                                                             accgTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTAACTGACTCGCTTCTCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGTTCGCTGGCGGGCCCTAAAGTTCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATGAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTcg            <  2:176099/248‑1 (MQ=255)                                                                                                                                                                                                                                        cTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGACACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGa  >  1:51192/1‑250 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        AAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTG‑TATCGACTGTTCGGTGCTTTTGGA‑TGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGA  >  NC_000913/4626252‑4626732

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