Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,521	T→C	L204L (TTG→CTG)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,521	0	T	C	94.8%	161.5 / ‑0.5	59	L204L (TTG→CTG)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (19/36);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.34e-01

CAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCG  >  NC_000913/4626275‑4626763                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    cagATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGac                                                                                                                                                                                                                                                 <  1:393101/251‑1 (MQ=255)      gAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCt                                                                                                                                                                                                                                            <  1:447006/251‑1 (MQ=255)      gAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCAGGCAACGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCt                                                                                                                                                                                                                                            >  1:162638/1‑251 (MQ=255)            gCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCACCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCa                                                                                                                                                                                                                                      <  1:368103/251‑1 (MQ=255)               gTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCaa                                                                                                                                                                                                                                     <  2:269782/249‑1 (MQ=255)                         ggCGAAGAGTCGCTGCACCAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGTATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGgtgt                                                                                                                                                                                                                         <  1:43163/251‑1 (MQ=255)                           cGAAGAGTCGCTTCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTTACGCGCTTCGCCAAAACGCGCGGTGTgg                                                                                                                                                                                                                       <  2:164512/251‑1 (MQ=255)                                   ggCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAAGCTCAATATGGTGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                                                                                                                                                               <  2:256174/250‑1 (MQ=255)                                   cGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                                                                                                                                                               <  1:39220/251‑1 (MQ=255)                                   cGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTc                                                                                                                                                                                                               <  1:205671/251‑1 (MQ=255)                                       ggaaCAGGGGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCGAAATATGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtgg                                                                                                                                                                                                           <  2:390621/249‑1 (MQ=255)                                             ggTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCAtggttggac                                                                                                                                                                                                     >  1:450199/1‑246 (MQ=255)                                                                cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATggct                                                                                                                                                                                  <  1:88115/250‑3 (MQ=255)                                                                cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTAAAGTCATCGCCGGGCAGCGTGGCGCAGGTACGTGAAACGGCGGCTTATCTGACGCGCTTCGCCACAACGCGCGGGGGGGCGATTGTCAtgggtggacacgagatcgggaaagcac                                                                                                                                                                                   >  1:433846/1‑226 (MQ=255)                                                                cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTc                                                                                                                                                                                 <  2:143392/251‑1 (MQ=255)                                                                cGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCAGCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATggct                                                                                                                                                                                  <  1:230845/250‑3 (MQ=255)                                                                cGCCTTGGCCTGCCGACTGACAATCTCAATAGGTTGGCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGAATCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAATCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACg                                                                                                                                                                                                 <  2:433846/235‑1 (MQ=255)                                                                                  gACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGt                                                                                                                                                               <  1:66870/251‑1 (MQ=255)                                                                                  gACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTCGCTGGCTGGCCCGaaa                                                                                                                                                                 <  1:141425/249‑1 (MQ=255)                                                                                   aCAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGG‑CACGTAACCAAAACTGGCTCGCTGGCGGGCCCCAAAggg                                                                                                                                                              >  1:81246/1‑249 (MQ=255)                                                                                          cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGa                                                                                                                                                         <  1:411433/249‑1 (MQ=255)                                                                                          cAATATGTAGTCGGAAACCAGCATCGAGCAGAGCTGCCTGATTGCCGAAGAAGAGCATCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTTCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGTCGATTGTCATGGTTGGA‑CTCGTAACCAAAGATGGCTCGCTGGCGGGTCCTAAAGTGCTg                                                                                                                                                           <  1:56120/247‑1 (MQ=255)                                                                                          cAATATGGTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTg                                                                                                                                                           <  1:418521/247‑1 (MQ=255)                                                                                              atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGCTCCAGGTGATGCATATGGCGGATGTACAGTCCTCGCCTTGCAGCGTGGCGCAGGTGCGTGAAACTGCGGCTTATATGACGCGCTTCGCCCAAACCCGCGGTGTGGCGCTTTTCATTGTTTGG‑CCCCTAACCCAAGATGGCTCCCTGGCGGGCCCTCAAGTGTTGGTACACtt                                                                                                                                                   >  2:38037/1‑250 (MQ=255)                                                                                                  tgtCGGAAACCGGCATCGATCAGATCTGCCTCGTTACCGAAGAAGAGCAACCGAAGCTGATGGTAAGTGACTCTGTCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAATGTGCTGGAACACTGTAt                                                                                                                                                <  2:217997/250‑1 (MQ=255)                                                                                                                aTCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGAAGCGCTTCGCCCAAACGCGCGGTGTGGCGGTTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGTCGGGCCCTAAAGGGCTTGAAAACTTTATAGGCTGTTCGGTGCt                                                                                                                                 >  1:204899/1‑251 (MQ=255)                                                                                                                     gcAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTgg                                                                                                                            <  1:284157/250‑1 (MQ=255)                                                                                                                             gccggATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATTCATATGGCGGAGGTACATTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCg                                                                                                                      <  2:283149/245‑1 (MQ=255)                                                                                                                                      ccGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTcc                                                                                                           <  1:206336/251‑1 (MQ=255)                                                                                                                                        gaagaaGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCa                                                                                                                                                                                                              <  1:390988/151‑1 (MQ=255)                                                                                                                                        gaagaaGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGGGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCa                                                                                                                                                                                                              >  2:390988/1‑151 (MQ=255)                                                                                                                                              gagCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCg                                                                                                   <  1:38037/251‑1 (MQ=255)                                                                                                                                              gagCAACCGAAGCTGATGGGAATTGACTCGAGCCAGGTGATGCATATGGCGGATGTACAGTCAGCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGTCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCGCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGGTCGGTGCGTTTGGATGGCGATGCCGACTCCCGTTTTCg                                                                                                   <  2:75371/251‑1 (MQ=255)                                                                                                                                                   aCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGGTTTTGGATGGCGATGCCGACTCCCGTTTTAGCActt                                                                                              >  1:304387/1‑249 (MQ=255)                                                                                                                                                       aagtggagggTAAGTGACTCGATACAGGTGATGCATATGGCGGAGTTACAGTCAGCGACGGGCAGCGTGGCGGAGGTGCGTGAAGAGGCGGCTTAGCTGACGCGCTTCGCCAAAACGCGAGGTGTGGCGATGGTCATGGGTGGA‑CACGTAACCAAAGATGGCTCGCGGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCTGTTTTCGCAcgttacc                                                                                          <  2:101316/243‑7 (MQ=255)                                                                                                                                                         gCTGATGGTAATTGACTAGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTTTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTTGAACACTGTATCGCCTGTTCGGGGCTTTTTGGTGGCGATGCCGACTCCCGTTTTCGCActtttcgc                                                                                         >  1:356192/1‑243 (MQ=255)                                                                                                                                                               gtgaaTTGACTCGTTCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGTA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACTCTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGTCTCCCGTTTTCGCACTTTACGCAGCCAt                                                                                   <  1:278031/247‑1 (MQ=255)                                                                                                                                                                      gacccGATCCAGGTGAGGCATATGGCGGATGTACCTGCACCGCCTGGCAGCGTGGCGCAGGGGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAcc                                                                            <  2:12537/246‑1 (MQ=255)                                                                                                                                                                         ccGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCtt                                                                        <  1:415845/250‑1 (MQ=255)                                                                                                                                                                          cGATCCAGGTGATGCATAGGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTc                                                                       <  2:86818/251‑1 (MQ=255)                                                                                                                                                                              ccAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCGGACGCGCTTCTCCAAAACGCGCGGTGTGTCGATTGTCATGGTTGGA‑CCCGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGgcg                                                                   <  1:281096/251‑1 (MQ=255)                                                                                                                                                                              ccAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGGGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGGACACTTTATAGACTGTTTGGTGCTTTTTGGTGGGCAGGCCGACCCCCGTTTTCGCACCTTTCGCAGCCATAAAAACCCCTTCgggg                                                                   >  2:263853/1‑249 (MQ=255)                                                                                                                                                                                  gTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCCAAAGTGCTTGAACACTGGATCGACTGTTCGGTGCTTTTGGGTGGCGATGCCGCCCCCCGTTTTCGGCCTTTACGCAGCCATAAAAACCGCTTCGGCGCggg                                                               >  1:70731/1‑250 (MQ=255)                                                                                                                                                                                                ggATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCg                                                 >  1:395397/1‑251 (MQ=255)                                                                                                                                                                                                ggATGTACAGGCATCGCCTGGCCGCGTGGCGCAGGTGCGTGAAACGGCTGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATGGTCATGGTTGGA‑CACGTAACCAAAGATGGCGCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATTCCGACTCCCGTTTTCTCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCg                                                 <  2:414208/251‑1 (MQ=255)                                                                                                                                                                                                           cATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATg                                      >  1:501662/1‑251 (MQ=255)                                                                                                                                                                                                              cGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCTACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAcc                                   >  1:176752/1‑251 (MQ=255)                                                                                                                                                                                                                    gTAGCGTGGCGCAGGTACGTGAAACGGCGGCTTATTTGACACGCTTCGCCACAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAAc                               >  1:242143/1‑249 (MQ=255)                                                                                                                                                                                                                     caGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGAt                                                                                                                    <  2:127043/162‑1 (MQ=255)                                                                                                                                                                                                                     caGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGAt                                                                                                                    >  1:127043/2‑163 (MQ=255)                                                                                                                                                                                                                      aGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAggg                           <  1:192747/251‑1 (MQ=255)                                                                                                                                                                                                                       gCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAgggg                          >  1:227028/1‑251 (MQ=255)                                                                                                                                                                                                                            gcgcAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCg                     <  1:73599/251‑1 (MQ=255)                                                                                                                                                                                                                               cAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTg                   >  1:357240/1‑250 (MQ=255)                                                                                                                                                                                                                                    gCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTc              <  1:14916/250‑1 (MQ=255)                                                                                                                                                                                                                                       tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGGAGGCAGc           >  1:353413/1‑250 (MQ=255)                                                                                                                                                                                                                                       tGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGc           >  1:396639/1‑250 (MQ=255)                                                                                                                                                                                                                                           aCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTTGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAAccc      >  1:128835/1‑251 (MQ=255)                                                                                                                                                                                                                                           aCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTTGGA‑CACGTAACCAAAGATGGCTCGCTGGCGGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTGTACGCAGCCATAAAATCCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAAccc      <  1:485927/251‑1 (MQ=255)                                                                                                                                                                                                                                               gggctgttttGACCCGCTTCGCAAAATCGGGCGTGGTGGGTATTGTCA‑GGTGGGGGCACGTAACCAAAGATGGTCCGCGGGCTGGCCCGAAAGTGCTGGATCCCTGTATCGACTTTTCGGTGCTTTTGGATGGCGATGCCGACCCCCGTTTTCGCCCCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTAGGCGTCTTCGCGATGACCGACCAGGGGCTTCGTGAAGTCAGCAACCCTTCg  <  2:242143/244‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGG‑CACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCG  >  NC_000913/4626275‑4626763

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