Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,948,964	G→A	A520V (GCT→GTT)	aspS ←	aspartate‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,948,964	0	G	A	94.9%	164.3 / ‑1.2	60	A520V (GCT→GTT)	aspS	aspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (23/33);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGT  >  NC_000913/1948717‑1949210                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gCAACGCTAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTGCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGt                                                                                                                                                                                                                                                     <  1:167375/251‑1 (MQ=255) gCAAAGCGAATGATTTCGGTTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGt                                                                                                                                                                                                                                                     <  2:516192/251‑1 (MQ=255)     cgcgAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAg                                                                                                                                                                                                                                                 <  1:55393/251‑1 (MQ=255)     cctgAATGATGTCGTGTATTGGAGTAAGAACAGTTCTTCCCATCCTCCTTCACAACCTGAATTATCATCTCAGCCAGCGCAGTCGGGTTAGCAAAGCTCGGTGCTTAAGTCATCAGACACGTTGCCGCCGTGGTTTTCGGGTAGTCGATAACGTCACGGATATTGTCGTTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGGACCGTATCTCATAACGTCGat                                                                                                                                                                                                                                                 <  2:362043/248‑2 (MQ=255)      gcgAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGc                                                                                                                                                                                                                                                >  1:89715/1‑251 (MQ=255)        gAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGc                                                                                                                                                                                                                                                >  1:111051/1‑249 (MQ=255)            gATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGACAGACGGTCAAGACCGAAAGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGaa                                                                                                                                                                                                                                           >  1:345198/1‑250 (MQ=255)               ttCGTGTAGTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCGGAAGGCTCAGCGCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACTCGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGcc                                                                                                                                                                                                                                        <  2:75776/250‑1 (MQ=255)                 cGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGa                                                                                                                                                                                                                                      >  1:430843/1‑250 (MQ=255)                  gtgtATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAAt                                                                                                                                                                                                                                    >  1:219959/1‑251 (MQ=255)                  gtgtATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGAGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGCCAAGACCGGATGCCAGACCTGCGTGCGGCGGGGCACCGTATTTCAGAACGTCGAGCAGGCAGCCGAAt                                                                                                                                                                                                                                    >  1:396264/1‑251 (MQ=255)                           aatCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCGCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCgcgc                                                                                                                                                                                                                           <  2:12671/249‑1 (MQ=255)                              cATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGGCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTAACGAATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGt                                                                                                                                                                                                                        >  1:375105/1‑251 (MQ=255)                                     gTTATTCTCAGCCTTCTTCACAACCGGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCtt                                                                                                                                                                                                                 <  2:399221/251‑1 (MQ=255)                                        aTTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGt                                                                                                                                                                                                              >  1:225664/1‑251 (MQ=255)                                             cAGCCTTTGTCACACCCAGAATGCTCGTCTCAGCCAGTGCCGGCGGGTTAGGCAAGCTCGTTGCTTCATTCATCATACACGCTGCCGCAGTGGTTTTCGGGAAGGCGATAACTGCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATa                                                                                                                                                                                                         <  2:418922/251‑1 (MQ=255)                                             cAGCCTTCTCCACAACCTGAATGCGCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGAt                                                                                                                                                                                                          <  1:346561/250‑1 (MQ=255)                                                 cttcttCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGAGGAGTACCGTATTTCAGAACGTCGAGCAg                                                                                                                                                                                                                                              <  1:445576/210‑1 (MQ=255)                                                   gcttcACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCGCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCa                                                                                                                                                                                                     <  1:290095/248‑1 (MQ=255)                                                        ccaACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTAGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATAcc                                                                                                                                                                                              <  1:115251/250‑1 (MQ=255)                                                         caACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATAcc                                                                                                                                                                                              >  1:74193/1‑250 (MQ=255)                                                            ccTGAAGGCTAAGCTCACCCAGTGCAGGCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAac                                                                                                                                                                                          <  2:541550/251‑1 (MQ=255)                                                               gattgctagtctcggcCAGTGCAGGCGGGTCAGCAAAGCTCGGTGCTTCAGTCACCAGACACGCGGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGTAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACAcc                                                                                                                                                                                       <  2:291913/237‑1 (MQ=255)                                                                    ctcagctcagCCAGTGCCGTCGGGTTAGCAAAGATCGGTGCTTAAGTCATCAGACAAGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGCAATTGTCGGTGTCGGTCAGCAGCATTGTCAGACGGCCAAGACCGGATGCCAGACATGCGTGCGGCGGCGTACCGTATTTACGAACGTCGAGCAGGAAACCGAATTTCTCGCGCTGTTCCTCTTCGTTGGTAccccgcatcccaaccaccctcgg                                                                                                                                                                                  >  2:538899/1‑231 (MQ=255)                                                                      cagctcagcCAGCGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGAACGAATGCCAGACCTGCGCGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTATCGCGCTGTTCCTCATCGTTGATACCCAGataatcaaaacccgtctgct                                                                                                                                                                                >  1:331221/1‑232 (MQ=255)                                                                       agctcagcCAGTGCAGTCGGGTTAGCAAAGCTCGGGGCTTCAGTCAGCAGACACGCTGCCGCCGTGGTTTGCGGGAAGGCGATAACTACACGGATATTGGCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTctgctg                                                                                                                                                                               <  2:26510/251‑1 (MQ=255)                                                                          tcagcCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGGATTTCTCGCGCTTTTCCTCTTCGTTGATACCCAGagaacgggagagcggcgtgtaggg                                                                                                                                                                            >  2:442402/1‑228 (MQ=255)                                                                          tcagcCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAg                                                                                                                                                                                                    <  1:442402/227‑1 (MQ=255)                                                                               cagcgcagggtggTTAGCAAAGCTCGTTGCCTCAGTCATCAGACCCGCTGCCGCCGTGGTTTTCGGGTAGGCGATAACGTCACGGATATTGTCGGTGCCTGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCATACCTGCGTGCGGCGGATTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCATCTTCGTTGATACCCAGAATACCAAACCCCGTCTGCTGCCTATCAc                                                                                                                                                                       <  2:340943/240‑1 (MQ=255)                                                                               cAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCAc                                                                                                                                                                       >  1:244600/1‑251 (MQ=255)                                                                                  tGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCAt                                                                                                                                                                    >  1:309074/1‑251 (MQ=255)                                                                                      gTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATg                                                                                                                                                                <  1:28547/251‑1 (MQ=255)                                                                                           gtgaGCAAACCTCGGAGCCGCAGTCATCAGACAAGCTGCAGACTTGGTTTGCGGGAATGCGAGAACGTCACGGATATTGTCGGGGACGGTCAGTAGCATGGTCAGACGTTCAAGACCGAAGTCCAGTCATGCGTGCGGCGGATTACCGTAGTTCAGAACGTCACTCAGGAAGCCGAATTGCTCGCGCTGTGCCTCTTCGGTGATACCCAGAATACCAAACACCGTCTGCTGCATATTCCCATTATGGATAc                                                                                                                                                           <  2:58192/248‑1 (MQ=255)                                                                                                cAAAGCTCGGTGCTTCAGTCATCAGAAACGCTGCCGCCGTGGTTTTCGGGGAGGCGATAACGTCCCGGATATTGTCGGAGCCGGTCTGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTGTTCGTTGTTACCCAGAATACCAAAATCCGTATGCTGCATATCACCATTATGGATActgact                                                                                                                                                      >  2:12621/1‑246 (MQ=255)                                                                                                                       aGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAAcc                                                                                                                               <  1:282212/251‑1 (MQ=255)                                                                                                                          cacGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCAtt                                                                                                                            <  1:119941/251‑1 (MQ=255)                                                                                                                             gCTGCCGCCGTGGTTTTCGGGAAGGCGTTAACGTCACGGATATTGGCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATtgat                                                                                                                         <  2:377414/251‑1 (MQ=255)                                                                                                                                    ccgTGGTTTTCGGGAAGGCGATAACGGCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCata                                                                                                                  <  2:424227/251‑1 (MQ=255)                                                                                                                                      gttggtttCGGGAAGGCGATTACGTTACGGATATGGTCGGTGCCGGTCAGCAGCATGGTCAGACTGTCAATACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGATCGTCTAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATTCCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATc                                                                                                                <  2:490548/246‑1 (MQ=255)                                                                                                                                      gTGGTTTTCGGGAAGGCGATAACGTCACGGACATTCTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACAGATTGCCAGACCTACGTGCGGCGGAGCACCGTATTTCAGAACGTCGAGCAGCAAGACGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTAAACCAACGCCCACTTCGTAACCATTGATGACCATATc                                                                                                                >  1:459293/1‑251 (MQ=255)                                                                                                                                       tGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCTCCCACTTCGTAACCATTGATGACCATATCg                                                                                                               >  1:103002/1‑251 (MQ=255)                                                                                                                                        ggTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCg                                                                                                               <  1:98899/250‑1 (MQ=255)                                                                                                                                              cGGGAAGGCGATAACGTAACGGAGATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGGCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGt                                                                                                        <  2:89715/251‑1 (MQ=255)                                                                                                                                                        aTAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGc                                                                                              >  1:363995/1‑251 (MQ=255)                                                                                                                                                          aaCGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCAt                                                                                            <  1:456575/251‑1 (MQ=255)                                                                                                                                                             gTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCAtttt                                                                                         <  2:430843/251‑1 (MQ=255)                                                                                                                                                                  ggATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTTAGCGTTCGCCACCGCATTTTCCGGt                                                                                    >  2:398650/1‑251 (MQ=255)                                                                                                                                                                                ccGGTCATCAGCATGGTCAGACGGTCAAGACCGAAGTCCAGATCGGCGTGCGGCGGAGTTCCGTAGTTCAGAACGTCGAGCAGGAAGCCTAATTTCTCGCGCTGTGCCTCTTCGTTGATACCCAGAATACCAAACACCTTCTGTTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGAGGACCATATCGTAAGCGTTCGCCCCCGAATTTTCCGGTGCAGCTTTCAGTTc                                                                      <  2:402147/251‑1 (MQ=255)                                                                                                                                                                                ccGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCGTTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTc                                                                      <  1:53215/251‑1 (MQ=255)                                                                                                                                                                                   gTCAGCAGCATGGCAAGACGGCCAACACCGAATGCCAGACCTGCGCGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGCTCCTCTTCGTTGATACCCAGAATTCCAAACACCGTCTGCTGCATCTCACCCTTATGGATACGTACTGAACCACCGCCCACTTCGTAACAATTGATGGCCATATCGTACGCGTTCGCCACCCAATTTTACGGTGCAGCTTTCCGTTCTGc                                                                   >  2:226656/1‑251 (MQ=255)                                                                                                                                                                                          gcaTTGTCAGACGGTCAAGACCGAAAGCCAGACCGGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCCCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCTCCATGATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTc                                                            <  2:34563/251‑1 (MQ=255)                                                                                                                                                                                                 cAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCtt                                                     >  1:122390/1‑251 (MQ=255)                                                                                                                                                                                                               cGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCAACATTATGGATACGTCCTGAACCACCGACCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCCCCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATACCTTTCggtgaggtgatt                                       >  1:362174/1‑249 (MQ=255)                                                                                                                                                                                                                           cggcgtgcgGCGGAGTACCGTATTTCAGAACGTCGAGCATGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGGAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTGCTGCAGCCGTCATAACTTTCGGTGAGGTGACCGGATGGTGCAtt                           <  2:30869/249‑1 (MQ=255)                                                                                                                                                                                                                                 gcggcgGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTc                     <  1:278216/251‑1 (MQ=255)                                                                                                                                                                                                                                 gcggcgGAGTACCGTATTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACAACGTCTGCTGCATATCACCATTATGGATACGTACTGACCCACCGCCCACTTCGTAACCATTGATTACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCCGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGGGCCTTGCCGTc                     >  2:60389/1‑251 (MQ=255)                                                                                                                                                                                                                                  cggcggAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTCCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCa                    <  2:74193/251‑1 (MQ=255)                                                                                                                                                                                                                                      ggAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGGCATATCTTTCGGTTAGGTGAGCGGATGGTGCATTGCCGTCAGgcc                >  2:137457/1‑251 (MQ=255)                                                                                                                                                                                                                                              gTATTTCAGAACGTCGAGCAGGAAGCCGAATTGCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGGTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCAc        <  1:142711/251‑1 (MQ=255)                                                                                                                                                                                                                                                    cAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACcgtcgt  >  2:277817/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GCAACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAGACCTGCGTGCGGCGGAGTACCGTATTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCACCGTCGT  >  NC_000913/1948717‑1949210

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