Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,851,614	C→T	A252V (GCC→GTC)	ansA →	L‑asparaginase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,851,614	0	C	T	95.2%	125.5 / ‑1.2	42	A252V (GCC→GTC)	ansA	L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (22/18);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

TGCGTTTGCCTCTCCAAACCTTCCTCCGTTACTGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGA  >  NC_000913/1851369‑1851857                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    tGCGTTTGCCTTTCCAAACCTTCCTCCGTTACGGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTc                                                                                                                                                                                                                                                 <  2:185451/250‑1 (MQ=255)       tGCCTCTCCAAACCTTCCTCCGTTACTGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAg                                                                                                                                                                                                                                           <  2:323401/250‑1 (MQ=255)           tctcCAAACCTTCCTCCGTTACTGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAAtt                                                                                                                                                                                                                                      >  2:178563/1‑251 (MQ=255)                  aCCTTCCTCCGTTACTGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACaagaa                                                                                                                                                                                                                               <  1:9054/251‑1 (MQ=255)                            gTTACTGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATc                                                                                                                                                                                                                     <  2:11859/251‑1 (MQ=255)                              tACTGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATcgc                                                                                                                                                                                                                   >  2:303393/1‑251 (MQ=255)                                cTGGAGGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCAGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCgg                                                                                                                                                                                                                 >  1:44940/1‑251 (MQ=255)                                 tGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGt                                                                                                                                                                                                                >  2:130481/1‑251 (MQ=255)                                  gaaaGCAGGTATCAATATTCGGCGTTTAAATACGCCAGCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGGCGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATGGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATGACAAGAAGTCAGCGATCGCGGTa                                                                                                                                                                                                               <  2:59567/249‑1 (MQ=255)                                    aaGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTAtt                                                                                                                                                                                                             >  2:309357/1‑251 (MQ=255)                                              caattatcgtcgTTGGAATACGCCTGCCGCACCGCACGGTGAAGGGGAATTAATCGTTCAACCAATCTCCCCACAACCAATTGGCGTAGTGACGAGTTATCCATGGATTTCTGGTGAAGTCGTGCGCCATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACATAACAAATTCTTCCTGCAGGAATTACAAGAAGGCAGCGAGCGCTGTATTGTGGTGGTCa                                                                                                                                                                                                   <  2:98156/245‑1 (MQ=255)                                                 tatTCGACGGTTGAATACGCCACTCGTACCGCACGCTGAAGGGGAATTAAAAGTACATCCAATCACCCCACAACCAATTAGAGTAGTAACGATTAATCCAGGGATTTATGCTGACGTCGTGCGCAATTTTATGCCCCAACCGGTGACAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTATTCCTGCAGGAATTCGAACACGACAGCGATCGCAGTATTGTGGTGGTCNAcc                                                                                                                                                                                                >  2:374394/1‑251 (MQ=255)                                                        cgtTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGacac                                                                                                                                                                                          >  2:162332/1‑250 (MQ=255)                                                                    ccACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCgg                                                                                                                                                                             >  2:200979/1‑251 (MQ=255)                                                                     cACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGt                                                                                                                                                                            <  1:309357/251‑1 (MQ=255)                                                                                cACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTTAACAt                                                                                                                                                                 <  2:22723/251‑1 (MQ=255)                                                                                 aCGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATg                                                                                                                                                                >  2:230778/1‑251 (MQ=255)                                                                                   ggTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTACCGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAAACGGTGAAAGAATTGATTCTGCGCTCATATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATggg                                                                                                                                                              >  1:260430/1‑251 (MQ=255)                                                                                     tGAAGGGGAATTAATCGTTCATCCAATCACCACACAACCAATTGGCGTAGTGACGCTTTATCCAGGGATTTCTGCTGACGTCGTGCGCAGTTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAAAATGggtg                                                                                                                                                            >  2:230906/1‑251 (MQ=255)                                                                                     tGAAGAGGAATTAATCGCTCATCCAATCAACCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGAGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCATGCAGGCATTACAAGAAGACAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCTGGTAAAGCGAACATGggtg                                                                                                                                                            >  2:349094/1‑251 (MQ=255)                                                                                                     gTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTa                                                                                                                                             >  1:121714/1‑250 (MQ=255)                                                                                                          tCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTAt                                                                                                                                                       >  1:384939/1‑235 (MQ=255)                                                                                                          tCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTAt                                                                                                                                                       <  2:384939/235‑1 (MQ=255)                                                                                                                cACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGccc                                                                                                                                 <  2:152813/251‑1 (MQ=255)                                                                                                                        aaCCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGccc                                                                                                                                 <  2:214108/243‑1 (MQ=255)                                                                                                                        aaCCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGccc                                                                                                                                 >  1:214108/1‑243 (MQ=255)                                                                                                                                   gTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCgg                                                                                                              <  1:130481/251‑1 (MQ=255)                                                                                                                                      gTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATTTATGTCCGGTAAAGTGAACATTGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGc                                                                                                           <  1:235283/251‑1 (MQ=255)                                                                                                                                                   cAGGGATTTCTTCTTTCGTCGTGCGCAATTTTCTGCGGCTACCGGTTAAAGCATTGATTCTGCGCTTCTATGTCTGGGGTAATGCGCCACAAAACAAAGTCTTCCGGCTGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCGGACTCAATGTATGTCCGTTAAAGTGAACATTGGTGGTTATGTCTCCGGTAACGTCCTCGCCCAGTCCGGCTTAATTTGGGGTGCAGATATGACTGTa                                                                                              <  2:260430/251‑1 (MQ=255)                                                                                                                                                      ggATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGaa                                                                                           <  1:131981/251‑1 (MQ=255)                                                                                                                                                       gATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGACCATGCCGTCGTAATTGGCGGTGCAGATTTGACTGTAGa                                                                                            >  2:147570/1‑249 (MQ=255)                                                                                                                                                                       gtGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACAt                                                                                                                                                                 <  2:87585/164‑1 (MQ=255)                                                                                                                                                                       gtGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACAt                                                                                                                                                                 >  1:87585/1‑164 (MQ=255)                                                                                                                                                                                      cgcCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCtt                                                             >  1:380239/1‑249 (MQ=255)                                                                                                                                                                                                      cATTGATTCTGCGCTCCTAGGTCGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACCAGAAGCCTGCGATCGCGGTATTGTGGTGGTCACCCTGACACAATGTATGTCCGGTAAAGTGAACATGTGTTGTTATGCCACCGGTACCGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACGATAGAAGCCACGCTATCCAACCTGCATTACCTGCTTAGCCCTGAACTGGATACt                                           <  2:300000/251‑1 (MQ=255)                                                                                                                                                                                                                      cTATGGCGTGGGTAATGCGCCACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGAACTGACTGTAGAAGCAACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAgg                           >  1:160557/1‑251 (MQ=255)                                                                                                                                                                                                                        ctGGCGTGGGTAGTGCGACACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGAGGTTATGCCACCGGTAACGCCCTCGACCATGCCGGCGTAATTGGCGGTGCAGATATGACCGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCAAGGAACTGGATACTGAAAACATTCGCAAggac                         >  2:67482/2‑249 (MQ=255)                                                                                                                                                                                                                            cGTGGGTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATg                      >  1:27917/1‑250 (MQ=255)                                                                                                                                                                                                                                ggTAATGCGCCACAAAACAAAGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGcc                  >  1:243132/1‑250 (MQ=255)                                                                                                                                                                                                                                    aTGCGGCACAAAACAAAGTGTTCCTGCAGGAATTACAAGGAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTGGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGGCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAAc             <  2:229220/251‑1 (MQ=255)                                                                                                                                                                                                                                              taacaaaGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGAATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGACGGCGTAATTGGCGGTGCAGATATGAATGTAGAAGACACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCTGAAGGCCATGAGCCAAAACCTGCgcggcc   >  2:90246/2‑250 (MQ=255)                                                                                                                                                                                                                                                acaaaGTCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGa  >  2:361048/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    TGCGTTTGCCTCTCCAAACCTTCCTCCGTTACTGGAAGCAGGTATCCATATTCGTCGTTTGAATACGCCACCCGCCCCGCACGGTGAAGGGGAATTAATCGTTCATCCAATCACCCCACAACCAATTGGCGTAGTGACGATTTATCCAGGGATTTCTGCTGACGTCGTGCGCAATTTTCTGCGCCAACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGA  >  NC_000913/1851369‑1851857

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