Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,878,077	C→T	intergenic (+459/+492)	iap → / ← cas2	alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,878,077	0	C	T	97.2%	109.5 / ‑4.5	36	intergenic (+459/+492)	iap/cas2	alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (12/23);  total (12/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

AACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTT  >  NC_000913/2877838‑2878318                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  aaCGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGTCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCgg                                                                                                                                                                                                                                        <  1:252702/251‑1 (MQ=255)   cgcgGGGAACTCGTAGTCCTTCATCCCCCCTATGTCTGATCTCCCGGTTTATCTCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCAGGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCTCGGGGAACTCAAGCTGGTTGGCAA‑TCTCTGTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCgggg                                                                                                                                                                                                                                      <  1:410556/251‑1 (MQ=255)          aaCTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTa                                                                                                                                                                                                                               <  1:473107/251‑1 (MQ=255)                 aGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTcc                                                                                                                                                                                                                         <  1:363029/250‑1 (MQ=255)                           tCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCgg                                                                                                                                                                                                               <  1:112804/250‑1 (MQ=255)                                    gTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAg                                                                                                                                                                                                     <  1:97069/251‑1 (MQ=255)                                                   ttATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCCGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTCGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTCTccc                                                                                                                                                                                      >  1:29559/1‑251 (MQ=255)                                                   ttATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATcc                                                                                                                                                                                       <  1:190204/250‑1 (MQ=255)                                                     agccccGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGTTTTATCCCCGCTGGCGCGGGGACCTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGGTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCg                                                                                                                                                                                    <  2:518608/249‑1 (MQ=255)                                                                    gggAACCCCCGGGGGATCATGTTTTCGGTCATGCGCCCCCCGGTTTACCCCCGGTGGCGCGGGGAACTCGGGGCGGCTGGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGGTTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACt                                                                                                                                                                     <  2:501419/251‑1 (MQ=255)                                                                                 ggATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGgc                                                                                                                                                                                                                                            <  2:250730/167‑1 (MQ=255)                                                                                 ggATAATGTTGACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGgc                                                                                                                                                                                                                                            >  1:250730/1‑167 (MQ=255)                                                                                  gATAATGTGTACGGTCATGCGCCCCCCGGTTTACCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGAcg                                                                                                                                                       <  2:540994/251‑1 (MQ=255)                                                                                            aCGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTAt                                                                                                                                             <  1:252231/251‑1 (MQ=255)                                                                                                    gcgcCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGAtt                                                                                                                                     <  1:458732/251‑1 (MQ=255)                                                                                                           ccGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGtt                                                                                                                              >  1:34788/1‑251 (MQ=255)                                                                                                                 gatacccGCTGGCGCGGGGAACTCGGGGCGGCTGGCCTGGCAGCCAGCGCCAGCCGCGGTTTATCCCCGCGGGCGCGGGGAACGCAAGCTGGCTGGCAA‑GCTCCTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCCCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCTGTTTATcc                                                                                                                         <  2:104108/246‑1 (MQ=255)                                                                                                                    ccccGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGc                                                                                                                     >  1:10565/1‑251 (MQ=255)                                                                                                                               cgGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCgggg                                                                                                           >  1:427160/1‑250 (MQ=255)                                                                                                                                                  tgccttgcAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGaa                                                                                        <  1:526500/250‑1 (MQ=255)                                                                                                                                                   gccttgcAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAAt                                                                                      <  1:423136/251‑1 (MQ=255)                                                                                                                                                              agcTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGCCGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACCTGAGACCGCTCAGACATTCCAGCCCCg                                                                            >  2:235378/1‑250 (MQ=255)                                                                                                                                                                cTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGTGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATc                                                                         >  1:55963/1‑251 (MQ=255)                                                                                                                                                                cTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGGCCGCTCAGAAATTCCAGACCCGATc                                                                         >  2:534964/1‑251 (MQ=255)                                                                                                                                                                  ccAGCAGCGGTTGATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAT‑GCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCa                                                                       <  1:315355/251‑1 (MQ=255)                                                                                                                                                                       agcGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAggcg                                                                                                                                                             >  1:94114/1‑160 (MQ=255)                                                                                                                                                                       agcGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAggcg                                                                                                                                                             <  2:94114/160‑1 (MQ=255)                                                                                                                                                                          agTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGCTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCCGGGTATGCGCGATTCGCGGTTTATCCCCCCTGGCGCGGGGAACTAGCTACCGCTCCGAAATTCCAGACCCGCTCCACAAGGttt                                                               >  2:432193/2‑251 (MQ=255)                                                                                                                                                                                cccccGCTGGCGCGGTGAACTCAAGCTGGCTGGCAA‑TCTCTTGCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATcccc                                                         <  2:34788/250‑1 (MQ=255)                                                                                                                                                                                 ccccGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCg                                                        <  2:284138/251‑1 (MQ=255)                                                                                                                                                                                      cTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGgc                                                   >  1:419005/1‑251 (MQ=255)                                                                                                                                                                                           gcgGGGACCCCAAGCTGGCTGGCAAGTC‑CTTTCGGGGTGAGGCCGGTGTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCGGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCggg                                              <  2:29559/251‑1 (MQ=255)                                                                                                                                                                                                gAACCCAAGCTGGCTGGAAA‑TCTCTGTCGGTGTGAGTCCGGTTTATCCTGGCTTGCGCGGTGAACTCTAGTTTCCGTATCGCCGGATTTATAAAGCAGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGAGTCGCGGTTTATCCCCGCTGGCGCGGGGAACCCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCTGGGAAct                                         <  2:237221/251‑1 (MQ=255)                                                                                                                                                                                                                          gtCGGGGTGAGGCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCg                  <  1:136973/248‑1 (MQ=255)                                                                                                                                                                                                                                      ccGGTTTATTCTCGCGGGCGCGGGGAACGCTAGTTGCCGTAGCCCCGGATATATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAAATCGCAGGCGGCGGCGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGt    <  2:152487/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGAGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGACATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACCTTATCAATTACAACCGACAGGGCACCCGGttt  >  2:561598/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  AACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAA‑TCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTT  >  NC_000913/2877838‑2878318

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