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	97.9%	165.5 / NA	49	A520V (GCT→GTT)	aspS	aspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/23);  minor base C (1/0);  total (26/23)
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.99e-01

AACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCA  >  NC_000913/1948719‑1949203                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                aaCGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCg                                                                                                                                                                                                                                            <  1:201498/251‑1 (MQ=255) aaCGCGAATGATTTAGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTCGCAAAGCTCGGTGCTTCAGTCAACCGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAAAGTCACGGATATTGTCGGCTTCGGTCAGTAGAATGGACAGCCGGTCCAGACCGTATGCCAG‑AACTGCGTGCGGCGGAGACCCGTA‑TTTAAGACCGt                                                                                                                                                                                                                                              >  2:91000/1‑249 (MQ=255)                            cATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGACGGGATAGCAACGCTCGGTGCTTCCGTCATCAGCCACGCCGCCGCCGCGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCAAGGTCAGCCGGTCAAGACCGAATGCCAG‑ACATGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGCAGCCGAATTTCTCGCGCTGt                                                                                                                                                                                                                 >  2:424652/1‑251 (MQ=255)                               atCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTcc                                                                                                                                                                                                              >  1:54905/1‑251 (MQ=255)                                     tATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCAAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCACTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATACCGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGACAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCg                                                                                                                                                                                                        >  2:207096/1‑251 (MQ=255)                                      aTTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGt                                                                                                                                                                                                       >  1:141506/1‑251 (MQ=255)                                          tcAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGTGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGCGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGAt                                                                                                                                                                                                   >  2:12287/1‑251 (MQ=255)                                           cAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCTAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGa                                                                                                                                                                                                    <  2:354282/249‑1 (MQ=255)                                           cAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTg                                                                                                                                                                                                     <  1:149487/248‑1 (MQ=255)                                             gCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATAcc                                                                                                                                                                                                >  2:342598/1‑251 (MQ=255)                                                ttcttcACAACCTGAATGCTCAGCTGAGACAGTGCAGTCGGGTTAGCAAAGATCGGCGCTTGAGTCAACAGACACGAAGCCGCCGTGGGTTTCGGGCAGGCGATAACGACACGGATAATGTCGGTGCCGTTACGCAGCATCGTCCGCAGGTAAAGACCGAATGGCAGGCCCTG‑GTGCGGCGGAGTACCGTACTTGCGGA‑CGTCGGGCCGGACACCGAATTTCTCGCGCTGTTCCTCTTCGTTGATAccc                                                                                                                                                                                               >  2:258832/1‑249 (MQ=255)                                                     cacaACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATAc                                                                                                                                                                                        >  1:414925/1‑251 (MQ=255)                                                       caACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCCAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATAccc                                                                                                                                                                                      >  2:165748/1‑250 (MQ=255)                                                            tGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAacac                                                                                                                                                                                 >  1:345798/1‑251 (MQ=255)                                                             gAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACAcc                                                                                                                                                                                >  1:173534/1‑251 (MQ=255)                                                                         cagcCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCata                                                                                                                                                                    <  1:194096/251‑1 (MQ=255)                                                                         cagcCAGTGCAGTCGGGTTAGCAAAGATCGGGGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGGGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGGTGATACCCAGAATACCAAACACCGTCTGCTGCat                                                                                                                                                                     >  2:50874/1‑250 (MQ=255)                                                                          agcCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCatat                                                                                                                                                                   <  2:363480/251‑1 (MQ=255)                                                                          agcCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCatat                                                                                                                                                                   >  1:281203/1‑251 (MQ=255)                                                                                     tCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATgg                                                                                                                                                        >  2:68337/1‑251 (MQ=255)                                                                                                      ggTGCGTCAGTCAGCAGACACGCTGCCGCCGTGGTTGTCGTGAAGGCGATAACGTCACGGATATTGGCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑GCCTGCGGGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCTAATTTCGCGCGCTGTTCCTGTTCTTTGATACTCAGAATACCAAACACCGTCGGCTGCATATCACCATTATGGTTACGTACTGAgccacc                                                                                                                                       <  2:439782/251‑1 (MQ=255)                                                                                                           ttCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATg                                                                                                                                                         >  2:34926/1‑228 (MQ=255)                                                                                                           ttCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATg                                                                                                                                                         <  1:34926/228‑1 (MQ=255)                                                                                                            tCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGccc                                                                                                                                   <  1:248739/249‑1 (MQ=255)                                                                                                                tcatcaGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGCGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTAATGAACCACCGCCCACTTCg                                                                                                                             >  2:291711/1‑251 (MQ=255)                                                                                                                   ccaGACACGCTGCCGCCGTGGTTTTCGGGAAGGTGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACTCCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTaa                                                                                                                          <  1:103884/250‑1 (MQ=255)                                                                                                                               ccgccgTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGAcc                                                                                                              <  2:368023/251‑1 (MQ=255)                                                                                                                                 gccgTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCat                                                                                                            <  1:173619/251‑1 (MQ=255)                                                                                                                                             gggAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGCTGACCATATCGTAAGCGtt                                                                                                >  2:234904/1‑251 (MQ=255)                                                                                                                                              ggAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTc                                                                                               <  1:305027/251‑1 (MQ=255)                                                                                                                                              ggAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTc                                                                                               <  1:243957/251‑1 (MQ=255)                                                                                                                                                 aGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGcc                                                                                            >  1:210362/1‑251 (MQ=255)                                                                                                                                                     gATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCg                                                                                        >  2:19042/1‑251 (MQ=255)                                                                                                                                                     gATAACGTCACGGATATTGTCGGTGCCGGTCAGAAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCGGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCg                                                                                        >  1:36416/1‑251 (MQ=255)                                                                                                                                                          cGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAacac                                                                                                                                                                                 >  1:341178/1‑157 (MQ=255)                                                                                                                                                          cGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAacac                                                                                                                                                                                 <  2:341178/157‑1 (MQ=255)                                                                                                                                                          cGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCAttt                                                                                   <  2:144062/251‑1 (MQ=255)                                                                                                                                                          cGTCACGGATATTGTAGGTGCCGGTCAGCATCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGACGAACATAACGTCAGCGGTCTCCCCCGCAttt                                                                                   >  2:21072/1‑251 (MQ=255)                                                                                                                                                                ggATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGt                                                                             <  1:131905/251‑1 (MQ=255)                                                                                                                                                                   tatTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGc                                                                           >  1:399334/1‑250 (MQ=255)                                                                                                                                                                              ccGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑GCCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTc                                                               <  2:173534/251‑1 (MQ=255)                                                                                                                                                                                                   cGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCgg                                          <  1:171914/251‑1 (MQ=255)                                                                                                                                                                                                         aGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCggtgaggt                                    >  2:164735/1‑251 (MQ=255)                                                                                                                                                                                                         aGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCggtgaggt                                    <  1:144417/251‑1 (MQ=255)                                                                                                                                                                                                         aGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTGTCggtgaggt                                    <  2:414925/251‑1 (MQ=255)                                                                                                                                                                                                              gAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACg                               >  2:289604/1‑251 (MQ=255)                                                                                                                                                                                                                            gcgtgcgGCGGAGTACCGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCAtt                    <  1:170129/249‑1 (MQ=255)                                                                                                                                                                                                                              gtgcgGCGGAGTACCGTA‑CTGCAGAACGGCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGGTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTACCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCGTTCGGTGAGGTGAACGGATGGTGCATTGcct                <  2:249847/251‑2 (MQ=255)                                                                                                                                                                                                                                            cGTA‑TTTCAGAACGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCCCCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTGCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCa  <  2:281203/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                AACGCGAATGATTTCGTGTATTTGAGTCATATCAGTTATTCTCAGCCTTCTTCACAACCTGAATGCTCAGCTCAGCCAGTGCAGTCGGGTTAGCAAAGCTCGGTGCTTCAGTCATCAGACACGCTGCCGCCGTGGTTTTCGGGAAGGCGATAACGTCACGGATATTGTCGGTGCCGGTCAGCAGCATGGTCAGACGGTCAAGACCGAATGCCAG‑ACCTGCGTGCGGCGGAGTACCGTA‑TTTCAGAGCGTCGAGCAGGAAGCCGAATTTCTCGCGCTGTTCCTCTTCGTTGATACCCAGAATACCAAACACCGTCTGCTGCATATCACCATTATGGATACGTACTGAACCACCGCCCACTTCGTAACCATTGATGACCATATCGTAAGCGTTCGCCACCGCATTTTCCGGTGCAGCTTTCAGTTCTGCAGCCGTCATATCTTTCGGTGAGGTGAACGGATGGTGCATTGCCGTCAGGCCGCCTTCA  >  NC_000913/1948719‑1949203

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