Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,823,288	C→T	G161D (GGC→GAC)	recA ←	DNA recombination/repair protein RecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,823,288	0	C	T	94.2%	187.3 / ‑1.3	70	G161D (GGC→GAC)	recA	DNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (42/23);  total (44/26)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCA  >  NC_000913/2823050‑2823531                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     aCCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCtt                                                                                                                                                                                                                                          >  1:538802/1‑251 (MQ=255)       cacGTTTTCGCCCTCTTTCCCCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGAttt                                                                                                                                                                                                                                    <  1:466704/251‑1 (MQ=255)       cacGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGAttt                                                                                                                                                                                                                                    >  1:147243/1‑251 (MQ=255)           ttttCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGc                                                                                                                                                                                                                                >  1:113197/1‑251 (MQ=255)              tCGCCCGCTTTCACCGGGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTGCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCCCACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTTTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCtt                                                                                                                                                                                                                              <  2:113197/250‑1 (MQ=255)                  ccTCTTTCCCCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGc                                                                                                                                                                                                                         <  1:398297/251‑1 (MQ=255)                    tctTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGt                                                                                                                                                                                                                       >  1:568256/1‑251 (MQ=255)                            cgcgCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGCTCTGGTTGATGCAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGCTTTCGCCTTCGATTTCCCCTTTCGGCGTCAGCgccc                                                                                                                                                                                                               >  2:73656/1‑250 (MQ=255)                            cgcgCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTgccg                                                                                                                                                                                                               >  1:230232/1‑251 (MQ=255)                            cgcgCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTgccg                                                                                                                                                                                                               >  1:544014/1‑251 (MQ=255)                                   aTACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTCCCGAACATCACACCAATTTTCATACGGATCTGGATGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATCCGTGCCGCAAGGCCCCTTTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAggcccgcaccgg                                                                                                                                                                                                         >  2:295692/1‑239 (MQ=255)                                         cgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAac                                                                                                                                                                                                  <  1:306318/251‑1 (MQ=255)                                              gTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAacgacg                                                                                                                                                                                              >  1:26625/1‑250 (MQ=255)                                                      gAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCt                                                                                                                                                                                     >  1:455810/1‑251 (MQ=255)                                                         cAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACt                                                                                                                                                                                  >  1:284865/1‑251 (MQ=255)                                                               cGTAGAATTTCAGCGCTTTACGACCGGTAGTGGTTTCCGGGTTACCGACCATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCTTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCa                                                                                                                                                                            <  1:458371/251‑1 (MQ=255)                                                                gTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGcc                                                                                                                                                                             >  1:37746/1‑249 (MQ=255)                                                                     aTTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAAcgc                                                                                                                                                                      <  1:493804/251‑1 (MQ=255)                                                                         cAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGcc                                                                                                                                                                   >  1:383080/1‑250 (MQ=255)                                                                             gcgTTTCCACGGGTAGGGGTTTCGGGGTTACCGAACATCACACAAATTTTCATAGGGATCTGGTTGATGAAGAGCAGCAGAGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGTTTACGCATCGCTTGGCTCATCATGCGTGCCTGAAGGCCCATGCGAGAGTCGTTGATTTCGCCTACGATTTCCTCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGc                                                                                                                                                              <  1:56040/251‑1 (MQ=255)                                                                              cgTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCg                                                                                                                                                             >  1:375289/1‑251 (MQ=255)                                                                                     accGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGACGTCGATTTCGCCTTCGCTTTCCGCTTTCGGCGTCAGGGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCAGAACGCGACAGGGAGTcaccat                                                                                                                                                      >  2:47040/1‑248 (MQ=255)                                                                                                        ttACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGcc                                                                                                                                   >  1:456444/1‑251 (MQ=255)                                                                                                              aaCATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTc                                                                                                                             >  1:314651/1‑251 (MQ=255)                                                                                                               aCATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGcc                                                                                                                                                                                                                                            >  1:487271/1‑139 (MQ=255)                                                                                                               aCATCACACCAATTTTCATACGGATCGGGTTGATGAAGATCAGCAGCGTTTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGcc                                                                                                                                                                                                                                            <  2:487271/139‑1 (MQ=255)                                                                                                                cATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCg                                                                                                                           >  1:254028/1‑251 (MQ=255)                                                                                                                cATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCg                                                                                                                           >  1:7668/1‑251 (MQ=255)                                                                                                                 aTCACACCACTTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGACGTCGATTTCGCATTCGATTTCCGCTTTCGGCATCAGTGCCGCCACGGAGTCATCGACGATAACGTCTACTGCGCCAGAACGCGCCAGGACGCAACAGATTTCACGTGCCTGATCGCCGGTGTCCgg                                                                                                                          >  2:553942/1‑251 (MQ=255)                                                                                                                 aTCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCgg                                                                                                                          >  1:98621/1‑251 (MQ=255)                                                                                                                   caccccAATTTTCATGCGGAGTGGGTGGATGAAGACCAGCGGC‑TGTTTGGACTTCTTCGGGTTCCCC‑GCCAGCTTACGCATCGCCTGGCCCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCt                                                                                                                        <  2:98621/247‑1 (MQ=255)                                                                                                                        cAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCCGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCAAGTGCCTGCTCGCCGGTGTCCGGCTggg                                                                                                                     >  1:72311/1‑249 (MQ=255)                                                                                                                        cAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCAGCACGCGCCAGGGCGTCACAGATTTCCAGTACCTGCTCGACGGTGTCCGGCTggg                                                                                                                     >  2:479593/1‑249 (MQ=255)                                                                                                                           ttttCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGcac                                                                                                                >  1:126849/1‑251 (MQ=255)                                                                                                                                   cGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGtt                                                                                                        >  1:122121/1‑251 (MQ=255)                                                                                                                                          ggTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCAGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGGCG‑at                                                                                                 >  1:553447/1‑250 (MQ=255)                                                                                                                                           gTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑atat                                                                                               <  1:528779/251‑1 (MQ=255)                                                                                                                                            ttGATGAAGATAAACAGCGTGTA‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCAACGCCTGGCTCATCATACGTGCCGCAAGGCCCCTGTGAGAGACGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGAGCCGCCACGGAGTCCACGACGATACCGTCTACTGCGCCAGAACGCGCCAGGTCGTCACAGATTTACATAGCCTGCTCGCCGGTTTCCGTCTGGTAGCCCAGAAGGATTTCG‑atcac                                                                                              >  2:446040/1‑248 (MQ=255)                                                                                                                                             tgatgaAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCg                                                                                             >  1:450236/1‑251 (MQ=255)                                                                                                                                                tgaAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACg                                                                                          >  1:289077/1‑251 (MQ=255)                                                                                                                                                  aaGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCCGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGcc                                                                                        >  2:225258/1‑251 (MQ=255)                                                                                                                                                        agcagcGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTAAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCCCAGCAGGTTGTCG‑AAATCGCCGCCCAgatt                                                                                  >  2:365765/1‑248 (MQ=255)                                                                                                                                                             cGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTg                                                                             >  1:509041/1‑251 (MQ=255)                                                                                                                                                              gtgtT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGc                                                                            >  1:231468/1‑251 (MQ=255)                                                                                                                                                              gtgtT‑GGACGGCATAAGGTTACCC‑GCCAGCTTACGCGTCGCCTGGCTCATCATACGTGCCGCAAGGCCCCTGTGACAGTCATCGATTTCGCCTTCGCTTTCCGCTTTCGGCGTCAGAGACGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCAGGACGAGCCAGGGCGTCACAGATGTCCAGTGCCTGCTCGCCTGTGTCCGGCTGGGTGCAAAGCCGTTGGTCG‑ATATCGACCCCCAGGTTACGTGc                                                                            >  2:535184/1‑251 (MQ=255)                                                                                                                                                              cgatcgggACTGCTGCAG‑TTACCCCGCCATCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCTATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTg                                                                             <  1:40910/245‑1 (MQ=255)                                                                                                                                                                    ggACTGCTTCAGGTTACAC‑GCCAGCTTACGCTTCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGCGAGTCGACGCTTTCGCCTTCGATTTCCGCTTTCGGCGTCCGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCCCCGGGGGCCGGCTGGGAGCACAGGCGTTTGTCGCAT‑TCGCCGCCCAGTTTACGTg                                                                             >  2:40910/1‑245 (MQ=255)                                                                                                                                                                      aCTGCTTCAGGTTACCC‑GCCAGCGTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGAtt                                                                     <  1:290474/251‑1 (MQ=255)                                                                                                                                                                       ccatagagtgttcccc‑GCCAGCTTACGCATCTCCTGGCTCATCATTAGTGCAGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTg                                                                    <  2:438121/238‑1 (MQ=255)                                                                                                                                                                         gctgaagtttcccc‑GCCAGCTTTCGCATCGCCTGGCTCATCATACGGGCCGCAAGGCCCATGGGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTggg                                                                  <  2:460143/240‑1 (MQ=255)                                                                                                                                                                         gcgtCAGGTTGCCC‑GCCAGCTTACGCATCGCCTGGCTCATCAGACGTGCCGCAATGCCCATGTGAGAGGCGTCGATGTCGCCGTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCGAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGGCT‑ATATCGACGCCCAGTTTACGTGCGTAGAtt                                                                     <  1:67591/245‑1 (MQ=255)                                                                                                                                                                         gCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTggg                                                                  >  1:2632/1‑251 (MQ=255)                                                                                                                                                                          cTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑CTATCGACGCCCAGTTTACGTGCGTAGATTGGGt                                                                 >  1:367421/1‑251 (MQ=255)                                                                                                                                                                            tcggcgtACCC‑GCCGTCTTACGCAGCGCCCGGCTCATCATACGTGCCGCCAGGCCCATGTGTGGGTCTTCCATTTCGCCTCTGATTTCCCCTTTCGGAGTCAGTGCCGCCACGGCGTCAACGACGATTACGTCTACTGCGCCAGTACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGTTGGGAGCACAGCAGGTTTTCG‑ATATCGACGCCCAGTGTTCGTGCGTAGATTGGGTc                                                                <  2:14977/244‑1 (MQ=255)                                                                                                                                                                             cAGGTTCCCC‑GCCAGCTTACGCATCGCCTGGCGCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGGCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACTTGCGTAGATTGGGttc                                                               <  2:230232/250‑3 (MQ=255)                                                                                                                                                                             cAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGCGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCAGCCACGGAGTCAACGACGATATCGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGCGCACCGCAGGTTGTCG‑ATATCGCCGCCCAGTTTACGTGCGTAGATTGGGTccc                                                              >  2:467802/1‑250 (MQ=255)                                                                                                                                                                               ggTTCCCC‑GCCAGCTTGCGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTTCGTAGATTGGGTCCAgc                                                            <  2:69913/251‑1 (MQ=255)                                                                                                                                                                                 ttACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAgcgc                                                          >  1:188517/1‑251 (MQ=255)                                                                                                                                                                                  tACCC‑GCCAGCTTACGCATAGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGcgcg                                                         >  1:306311/1‑251 (MQ=255)                                                                                                                                                                                          cAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAg                                                  <  1:436414/251‑1 (MQ=255)                                                                                                                                                                                           ggCTTACGCATCGCCTGGCTCAGCATACGTGCCGCCAGGCCCAGGGGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATGTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATACCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGc                                                 <  2:7668/250‑1 (MQ=255)                                                                                                                                                                                                       gCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGc                                     >  2:72752/1‑251 (MQ=255)                                                                                                                                                                                                               catcatACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGt                               <  1:486679/249‑1 (MQ=255)                                                                                                                                                                                                               cagcatACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTAAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTTCGTAGATTGTGTCCATCGCGTGTTCAGCATCGATAATCGCACAGGttt                             <  1:84568/248‑1 (MQ=255)                                                                                                                                                                                                                   atACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTAcc                         >  1:187953/1‑251 (MQ=255)                                                                                                                                                                                                                            gCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctg                <  1:472966/251‑1 (MQ=255)                                                                                                                                                                                                                             cAAGGCCCATGGTAGAGTCGTCGATTCCGCCTTCGATTTCCGCTTACGGCGTCAGTGCCGCCACGGAGGCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTGCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGAGTGGTTCCAGCGCGTGTTCAGCATCGATAATCGCACAGGTTTTACCTTCAcgcttc               <  2:39180/251‑3 (MQ=255)                                                                                                                                                                                                                                   ccgggggagaGTCGTCGATTTCTGCGTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGGCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGGCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑TTATCGACGCCCAGTTTTCGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTTACGctgcgctgcg         <  2:568256/245‑1 (MQ=255)                                                                                                                                                                                                                                       gtgAGAGTCGTCGTATTCGACTTCGATTTCCGCTTGCGGCGGCAGTGCCGCCAAGGAGGCAACGACGATTGCGTCTATTGCGCCATATCGCGCCAGGGCGTCACAGAGTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGGTGGTTGTCG‑AGATCCACACCCAGTCTACGGGCGTTGATTGGGTCCAGCGTGTTTTCAGGAGCGATGAACGCAGAGTttttgactgcactttacgctgcggc       <  2:340819/249‑22 (MQ=255)                                                                                                                                                                                                                                          agagTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTTCGTAGATTGGGTCCAGCGCGCGTTCACCATCGATAAACGAACAGGTTTTACCTTCACGCGGCGCTGTGGGGATcc  >  2:121256/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     ACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTT‑GGACTGCTTCAGGTTACCC‑GCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCG‑ATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCA  >  NC_000913/2823050‑2823531

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