Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,042	T→C	M44T (ATG→ACG)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,042	0	T	C	100.0%	140.4 / NA	46	M44T (ATG→ACG)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/18);  total (28/18)

GTCACCATCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAA  >  NC_000913/4625795‑4626282                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 gTCACCATCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAAcggt                                                                                                                                                                                                                                               >  1:210881/1‑251 (MQ=255)     ccATCCGTCACGCTGACCTGATGGGGGAATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGACATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGATCCAACACGCGCCTTTGTTTGTAATGACTGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGgcg                                                                                                                                                                                                                                           >  1:263997/1‑251 (MQ=255)      cATCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGgcgc                                                                                                                                                                                                                                          >  1:48437/1‑251 (MQ=255)             cACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGTATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAACCGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGTGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGAGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACg                                                                                                                                                                                                                                    >  2:51225/1‑250 (MQ=255)                     ccTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCACCGCGGGCGCGTACCGAGCGtctc                                                                                                                                                                                                                            >  2:130793/1‑250 (MQ=255)                              ggTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGTCACAGCATTAACGAAGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGATCACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATg                                                                                                                                                                                                                  <  1:147198/251‑1 (MQ=255)                               gTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGc                                                                                                                                                                                                                 >  1:204385/1‑251 (MQ=255)                                 aTTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGcc                                                                                                                                                                                                                <  1:46249/250‑1 (MQ=255)                                       cATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGTGCAGTGCAGTGCCTGTCATGCCTGGAACACCAGCACCGAGGTGCGTCGTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCGCATCGGCTATGCCGGTAgcg                                                                                                                                                                                                         <  2:22810/251‑1 (MQ=255)                                        atcctaccagccAGCCTGCATCAAAAATAATTGCTCGCCAGCCCTCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGAGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACCCCATCACCCAGGCGCGTCTTGCTGCGTCGCCAACGGTGGCGCGCAACGAGCGTCTCGGCGGCCATGCCGGTAgcgc                                                                                                                                                                                                        >  2:43592/12‑251 (MQ=255)                                               cAGGCAGCCTGAATCAGAAGTAATTGCGCGCCCGCCAGCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGt                                                                                                                                                                                                  <  1:78164/250‑1 (MQ=255)                                                  gCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGcc                                                                                                                                                                                              >  2:206040/1‑250 (MQ=255)                                                    aGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCGGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCaaa                                                                                                                                                                                            <  2:210881/251‑1 (MQ=255)                                                                aaGTAATAGCGCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAAct                                                                                                                                                                                >  1:125395/1‑251 (MQ=255)                                                                  gTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGGGGCAAAAGTCCAGCAActct                                                                                                                                                                              >  2:204216/1‑251 (MQ=255)                                                                     aTTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGAGGGGCCGATTATCCGCGCTGGCAGGGGAAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGGGGCAAAAGTCCAGAAACTCTCCg                                                                                                                                                                           >  1:126226/1‑251 (MQ=255)                                                                           cgcccgccATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGatat                                                                                                                                                                       <  1:151267/249‑1 (MQ=255)                                                                                     ccTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAAGGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCAAGAAACTCTCCGATATCAGCCTTGAag                                                                                                                                                            <  1:126272/250‑1 (MQ=255)                                                                                         cgggcggCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGTGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACCCCATCACCGAGGTGCGTCGTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGTCAAAAGGCCAGAAAGTCTCCGATATCAGCCTTGAATAGCTg                                                                                                                                                       <  2:243694/251‑1 (MQ=255)                                                                                               gcACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGt                                                                                                                                                 <  1:35836/251‑1 (MQ=255)                                                                                                                    ccGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAGGGCCAGAAACTCTCCGCTATCAGCCTTGAAGAGCTGCCGCGGTTTTCCACCGGATTTAAag                                                                                                                              >  2:984/1‑249 (MQ=255)                                                                                                                       tGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTc                                                                                                                         >  1:13496/1‑251 (MQ=255)                                                                                                                         gCAAAAGCTCCAAAACGCGCCTTTGTTTGTACTGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGAGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTCGCGCCGGGGTGGCAAAAGTCCAGAAACtctc                                                                                                                                                                             >  2:262191/1‑197 (MQ=255)                                                                                                                         gCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACtctc                                                                                                                                                                             <  1:262191/197‑1 (MQ=255)                                                                                                                                  ccAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGGGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACt                                                                                                              >  1:73837/1‑251 (MQ=255)                                                                                                                                               gtggtggtAATGAATGCGGGGCCGATTATCCGCGCGGGCAGGGGCAGTGCAGTGCCTGTCATGCCGGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGGAGCGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCGTCCGCGTTTTTCCACCGGATTTAAAGATTTCGACCGCGTGCTAGgcggcggcg                                                                                                   <  2:263997/243‑1 (MQ=255)                                                                                                                                                    tgtAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGc                                                                                             <  1:42815/250‑1 (MQ=255)                                                                                                                                                     gtAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCCTGCCTGGACCACCATCACCGAGGTGCGTCTTGCTGCGTCGCCACCGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGCGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGccc                                                                                           >  2:26848/1‑250 (MQ=255)                                                                                                                                                          gaatgCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCACCGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGCTTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGa                                                                                        >  2:155505/1‑249 (MQ=255)                                                                                                                                                                  ggCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAATTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTc                                                                              <  1:156624/251‑1 (MQ=255)                                                                                                                                                                    ccGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTg                                                                            >  2:66724/1‑251 (MQ=255)                                                                                                                                                                        ttATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTg                                                                        >  1:22222/1‑251 (MQ=255)                                                                                                                                                                        ttATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTg                                                                        >  1:165891/1‑251 (MQ=255)                                                                                                                                                                            ccGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCgg                                                                    >  2:2122/1‑251 (MQ=255)                                                                                                                                                                            ccGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGAGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTCAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGCAGTGCCATTCTGCTTGGCgg                                                                    >  2:126233/1‑251 (MQ=255)                                                                                                                                                                                gcTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAc                                                                >  1:41375/1‑251 (MQ=255)                                                                                                                                                                                gcTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAc                                                                <  1:67237/251‑1 (MQ=255)                                                                                                                                                                                  tGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGGGCCAGGAAGTGCCATTCTGATTGGCGGTAAccc                                                              >  2:54471/1‑251 (MQ=255)                                                                                                                                                                                      aGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTggg                                                          >  1:232743/1‑250 (MQ=255)                                                                                                                                                                                           cagtgcagtgCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCggg                                                     <  2:18931/251‑1 (MQ=255)                                                                                                                                                                                                  gtgcCTGTCATGCCTGGAACACCAGCTCCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATcc                                              <  2:165891/251‑1 (MQ=255)                                                                                                                                                                                                    gcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCa                                             <  1:39872/250‑1 (MQ=255)                                                                                                                                                                                                              ccTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgttgctg                                  >  2:188390/1‑245 (MQ=255)                                                                                                                                                                                                                           cACCGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGc                      >  1:146083/1‑250 (MQ=255)                                                                                                                                                                                                                              cGAGGTGCGTCTTGCTGCGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAAc                  <  1:26848/251‑1 (MQ=255)                                                                                                                                                                                                                                aGGTGCGTCTTGCTGCGTCGCCACCGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAAACGCTGTGCAAACTg                >  1:59959/1‑251 (MQ=255)                                                                                                                                                                                                                                              gcGTCGCCAACGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGCTTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGAACCCCTGGTGCGGGGAAATCCACGCTGTTGCTGCAACCGCTGGGCACACTGGCCCAGCAGATggc  >  2:212352/1‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 GTCACCATCCGTCACGCTGACCTGATGGGGGTATTCTGCATCCTCTCAGGCAGCCTGAATCAGAAGTAATTGCTCGCCCGCCATCCTGCGGGCGGCACAGCATTAACGAGGTACACCGTGGCAAAAGCTCCAAAACGCGCCTTTGTTTGTAATGAATGCGGGGCCGATTATCCGCGCTGGCAGGGGCAGTGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAA  >  NC_000913/4625795‑4626282

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