Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,798,171	C→T	D258N (GAC→AAC)	pheS ←	phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,798,171	0	C	T	95.6%	131.1 / ‑1.5	46	D258N (GAC→AAC)	pheS	phenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (22/21);  total (23/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 = 7.05e-01

GAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGATAAACACGGCCAGGCGCGATGATACGAATCGGTGGCTGCTGGGCTTTCATGGTGCG  >  NC_000913/1797936‑1798417                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       gAATTTCATAATCTATTCCTGCCTTATTTAAACGGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGaa                                                                                                                                                                                                                                           <  1:21445/249‑1 (MQ=255)          aaTCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAAATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTGTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTaaa                                                                                                                                                                                                                                <  1:180976/251‑1 (MQ=255)           aTCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAAc                                                                                                                                                                                                                               >  1:258038/1‑251 (MQ=255)               aTTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGa                                                                                                                                                                                                                           <  1:88150/251‑1 (MQ=255)                              aaaCTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACgg                                                                                                                                                                                                             <  1:167084/250‑1 (MQ=255)                                       gagtaaaCGCAGATCGTTTTCGAAGAAAGAACGCATGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGTTGCCAACGTTACGCAAC‑TCGTTCTGATGCACCACCCCGCAGCCCAGCGC‑TGCCAGCCATGTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTGCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCt                                                                                                                                                                                                   <  2:187944/247‑1 (MQ=255)                                             aCGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCAGCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAATCGAATCTGCAAAt                                                                                                                                                                                             <  1:79196/251‑1 (MQ=255)                                              cGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTACGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCACAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATc                                                                                                                                                                                            >  1:65167/1‑251 (MQ=255)                                                cAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGCTGCACCATCCCGCCGCCCAGCAC‑TTCCAGACATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCtt                                                                                                                                                                                          >  2:545624/1‑251 (MQ=255)                                                     cGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTc                                                                                                                                                                                      >  1:422967/1‑250 (MQ=255)                                                           cGAAGAACAAACGCAGGTCGGTGACGCCGTAACGCAACATAGTTAGACTCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACATCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGAGGCACCATCCCGCAGTCCAGCTC‑TTCCAGCCAGTTACCGTATTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTTCAAATCTTCCTCAaagaa                                                                                                                                                                                <  2:75770/250‑1 (MQ=255)                                                                     aCGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAgg                                                                                                                                                                      >  1:361400/1‑250 (MQ=255)                                                                           gTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGt                                                                                                                                                               <  1:46109/251‑1 (MQ=255)                                                                           gTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGt                                                                                                                                                               <  1:296824/251‑1 (MQ=255)                                                                             cGGTGACGCGGTAACGCAACAGAAGCAGACGGTCCATCCCCATCCCGAAAGCGAAACCAGAGTAAACTTCCGGGTCGAGGCCAACGTTAACCAAC‑TCGTTCGGATGCACCATCCTGCAGCCCAGCAC‑TTCCAGCCATTTACCGTGTTTACCCATGACGTGCACTTCTGCAGAAGGTGCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCCCAAAGAAGTTACGCAGGAAGTCGTg                                                                                                                                                              <  2:475529/250‑1 (MQ=255)                                                                              ggTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGc                                                                                                                                                             >  1:120609/1‑250 (MQ=255)                                                                                   cGCCGAAACGCAACATAGTCAGACGCACCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGGCGATGACAACGTTACGCAAC‑ACGTTCGGCTGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCGCTTCTGCAGAAGGTTCGGTAAACGGGAATTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCATGAATTCGTGCAGCGTg                                                                                                                                                       <  2:58773/251‑1 (MQ=255)                                                                                           cGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTc                                                                                                                                                 >  1:456120/1‑249 (MQ=255)                                                                                           agCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGACTCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTc                                                                                                                                                 >  1:456111/2‑249 (MQ=255)                                                                                                   aGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTa                                                                                                                                        >  1:425282/1‑250 (MQ=255)                                                                                                      cAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGc                                                                                                                                    >  1:304685/1‑251 (MQ=255)                                                                                                        gACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTg                                                                                                                                  >  1:410218/1‑251 (MQ=255)                                                                                                        gACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATACACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTg                                                                                                                                  >  1:327379/1‑251 (MQ=255)                                                                                                         ccGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTGCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGa                                                                                                                                 <  2:456120/250‑1 (MQ=255)                                                                                                         aCGCCCCACCCCCATCCCGAAGTCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTAACCATGACGTTCACTTCTGCAGAAGGTGCTGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGTAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGa                                                                                                                                 <  2:456111/251‑1 (MQ=255)                                                                                                          cgcgccatccccatccCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTGACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGAt                                                                                                                                <  1:448394/247‑1 (MQ=255)                                                                                                                 tccccatcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGt                                                                                                                          >  1:387445/1‑250 (MQ=255)                                                                                                                      atcccGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATcaac                                                                                                                    <  1:494852/251‑1 (MQ=255)                                                                                                                           gAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATc                                                                                                                >  2:108161/1‑250 (MQ=255)                                                                                                                                 gAAACCAGAGTAAACTTCCGGGTCGATGCAAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAACCGGGAAGTAGGACGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGCAGTCGTGCAGCGAGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATAAACAATCAGACCtt                                                                                                         >  2:556446/1‑251 (MQ=255)                                                                                                                                          gTAAACTTCCGGGTCAATGCCAACGTTACGCAACCAC‑TTCAGATGCACCATCCCGCAGCTCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTATGCAGCAGGCTCGGTAAACGGGAAGTAGGAAGGACGGAAGCCAATCTGCAAATATTCCTCAACGAAGTTACGCAGGATGTCGTGCAGCGTGCCTTTCAGGTTCGTAAAGCGAATGTTGGTATCAACAATCAGACCTTCCATCTGAt                                                                                                >  2:221755/1‑251 (MQ=255)                                                                                                                                                       tCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCgt                                                                                   >  1:162131/1‑251 (MQ=255)                                                                                                                                                               aaCGTTAAACAAC‑ACGAACGGATGCACCATCACGCCGCCCA‑CACCTTCCAGCCCTTTACCGTTTTAACCCATGACGTTAACTTCTGCAGAAGGTTCTGTAAACGGGAAGTAGGCAGGACGGACACGACTCGGCAAATCTTCCTCACAGACGTTCCGCAGGAAGGCGTGCAGCGTGCATTCCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCCTCTGAAGGAACATCGGCgattgagcct                                                                           >  2:567829/1‑242 (MQ=255)                                                                                                                                                                aCGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTg                                                                          <  1:430384/251‑1 (MQ=255)                                                                                                                                                                        cAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑GTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAAGCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAATCTGATGTTGGTAGCAACAATCAGACCTGCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTa                                                                    <  2:258038/249‑1 (MQ=255)                                                                                                                                                                                 tCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACg                                                           >  1:79528/1‑250 (MQ=255)                                                                                                                                                                                  cgtattctccaaccctcAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCAATTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGGCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGAt                                                         <  2:304685/235‑1 (MQ=255)                                                                                                                                                                                   ggATGCACCATCCCGCAGGCCAGTAA‑TTCCAGCCATTTACCGTTATTACCCATGACGTTCACGTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTGCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCAGCTGATGGAACATCGTCGTGTGAGGCTGGTCGTAGTCGTTACGATa                                                        <  2:79528/251‑1 (MQ=255)                                                                                                                                                                                         aCCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGGAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGATAAACACg                                                  <  2:120609/251‑1 (MQ=255)                                                                                                                                                                                           cATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTATGGTCGTAGTCGTTACGATAAACACGGc                                                >  1:107279/1‑251 (MQ=255)                                                                                                                                                                                                cGCAGACCAGCGC‑TTCCAGCCTGTTACAGTTTTTTCCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGTGAAGGAGGAAGGACGGAAGCGAATCTGCAAAGCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGGAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGATAAACACGGCCAGgc                                           <  2:152604/251‑1 (MQ=255)                                                                                                                                                                                                                    ccATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGATAAACACGGCCAGGCGCGATGATACGAATCGGTg                        >  1:220778/1‑251 (MQ=255)                                                                                                                                                                                                                     cATTTACCGTTTTTACCCATGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGATAAACACGGCCAGGCGCGATGATACGAATCGGTgg                       >  1:227017/1‑251 (MQ=255)                                                                                                                                                                                                                        ttACCGTTTTTACCCAGGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGATAAACACGGCCAGGCGCGATGATACGAATCGGTGgct                     <  2:220778/250‑1 (MQ=255)                                                                                                                                                                                                                                       aTGACGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGATAAACACGGCCAGGCGCGATGATACGAATCGGTGGCTGCTGGGCTTTCATGGt     <  2:422967/251‑1 (MQ=255)                                                                                                                                                                                                                                          aCGTTCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGCAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATAAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGGTACGACAAACACGGCCAGGCGCGATGATACGAATCGGTGGCTGCTGGGCTTTCATgggacg  >  1:501833/1‑247 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       GAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAAC‑ACGTTCGGATGCACCATCCCGCAGCCCAGCAC‑TTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAACAATCAGACCTTCCATCTGATGGAACATCGGCGTGTGAGTCTGGTCGTAGTCGTTACGATAAACACGGCCAGGCGCGATGATACGAATCGGTGGCTGCTGGGCTTTCATGGTGCG  >  NC_000913/1797936‑1798417

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