Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	673,386	C→A	P466P (CCG→CCT)	leuS ←	leucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	673,386	0	C	A	98.6%	246.8 / NA	70	P466P (CCG→CCT)	leuS	leucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (36/33);  total (37/33)
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.93e-01

CATCAGTTTGTGGAAGAAGCGGAAGTAGAGCAGGTGCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATCGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGTCAAGACCGTTGAACTCGC  >  NC_000913/673140‑673630                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     cATCAGTTTGTGGAAGAAGCGGAAGTAGAGCAGGTGCATAATGGCGAGTTCAATACCACCCATGTAGATAACCACCGGCAGCCAGTAGTTAGCCGCTTCAGAATCCAGCATAACTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGc                                                                                                                                                                                                                                                  >  1:312679/1‑251 (MQ=255)    cAGTTTGTGGAAGAAGCGGAAGTAGAGCAGGTGCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGGCGGTTTCACGCAGTGCTGGCATACCGTTAACGGGAGGTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTg                                                                                                                                                                                                                                                <  2:83071/250‑1 (MQ=255)    cAGTTTGTGGAAGAAGCGGAAGTAGAGCAGGTGCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAACCCGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCCAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCCCCCACTCCCGCTCTGCTTTAATAGGGCTgg                                                                                                                                                                                                                                               >  1:25360/1‑251 (MQ=255)        ttGTGGAAGAAGCGGAAGTAGAGCAGGTGCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCCGATCTGCTTTAATAGGGCTGGTAAt                                                                                                                                                                                                                                           >  1:119124/1‑251 (MQ=255)         tgtgGAAGAAGCGGAAGTAGAGCAGGTGCATAATGTCGTGTTCAATACCACCAATGTAGATAGCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAAt                                                                                                                                                                                                                                           <  1:11410/250‑1 (MQ=255)         tgtgGAAGAAGCGGAAGTAGAGCAGGTGCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAAt                                                                                                                                                                                                                                           <  1:332718/250‑1 (MQ=255)                  agCGGAAGTAGAGCAGGTGCATAATGGCGTGTTCAATACCCACAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGGTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCAtt                                                                                                                                                                                                                                 >  1:161373/1‑251 (MQ=255)                           agagCAGGTGCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCt                                                                                                                                                                                                                        <  2:222979/251‑1 (MQ=255)                             agCAGGTGCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCAGTACCACATCTTc                                                                                                                                                                                                                      <  2:221402/251‑1 (MQ=255)                                  gTGCATAATGGCGTGTTCAATACCACCCATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGc                                                                                                                                                                                                                  >  1:409850/1‑250 (MQ=255)                                    gCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCCTAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGCGTCGGTTTCACGCAGTGCTGGCATAACGTTAACGGTAGGTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAgg                                                                                                                                                                                                               >  1:224913/1‑251 (MQ=255)                                        aaTGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCTGCAGGATCa                                                                                                                                                                                                           <  1:348622/251‑1 (MQ=255)                                        aaTGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATc                                                                                                                                                                                                            <  1:102894/250‑1 (MQ=255)                                        aaTGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATc                                                                                                                                                                                                            <  2:344696/250‑1 (MQ=255)                                           ggCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCg                                                                                                                                                                                                        >  1:425054/1‑251 (MQ=255)                                                 ttCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCt                                                                                                                                                                                                  >  2:162501/1‑251 (MQ=255)                                                     aTACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGgtc                                                                                                                                                                                              >  2:315152/1‑251 (MQ=255)                                                       accaccAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGgtcg                                                                                                                                                                                             <  2:26264/250‑1 (MQ=255)                                                        ccaccaATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGgtcgtc                                                                                                                                                                                           <  1:325934/251‑1 (MQ=255)                                                         caccaATGTAGATATCAACCGGCAGACCGGAGTTAGACGATTCGGCATCCAGCATACCTTCTTTGTAATGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGCCTCCATAACGGTGTCGACAGTGTCGGTTTCACGCAGTGCTGGCATTCCGTTAACGGTAGTTTTCGCCCACTCCGGATATGTTTTAATAGGGCTGGTAATGTCGTCCATTACCACATCTTCCGGCAGGTTCACCggacagtggccgccc                                                                                                                                                                                          >  2:297546/1‑238 (MQ=255)                                                                 tAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCgg                                                                                                                                                                                  >  1:92571/1‑251 (MQ=255)                                                                  aGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACCTCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGc                                                                                                                                                                                 >  1:283534/1‑251 (MQ=255)                                                                  aGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGc                                                                                                                                                                                 >  2:347925/1‑251 (MQ=255)                                                                        ccACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCCTACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCCGCTGGTCGTCCGGGGTCGGCATTa                                                                                                                                                                             >  2:186436/1‑249 (MQ=255)                                                                                   cAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGgtcgtc                                                                                                                                                                                           <  1:136174/224‑1 (MQ=255)                                                                                   cAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGgtcgtc                                                                                                                                                                                           >  2:136174/1‑224 (MQ=255)                                                                                   cAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTc                                                                                                                                                                >  1:156250/1‑251 (MQ=255)                                                                                   cAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTc                                                                                                                                                                <  1:93809/251‑1 (MQ=255)                                                                                                  tCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCgg                                                                                                                                                 >  2:281627/1‑251 (MQ=255)                                                                                                      aaTCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATc                                                                                                                                             >  1:366090/1‑251 (MQ=255)                                                                                                           aGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGCCTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGGCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCCTCGGAATCGgcgc                                                                                                                                        >  2:390554/1‑251 (MQ=255)                                                                                                                aCCTTCTTTGTACTGCGCGCAAGTGAAGCGCGCATAGCCCCAGGAGGACTCCATACAGGTGTCGAAAGTGTCGGTTTCACGGCGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGGTCTGCTTTAATCGGGCTGGTACTGCCGTCCATTACCCCATCTTCCGGCCGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGGACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGcccc                                                                                                                                   >  2:332017/1‑251 (MQ=255)                                                                                                                       ttGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTaa                                                                                                                              <  1:167963/249‑1 (MQ=255)                                                                                                                         gTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAACGGTGTCGCAAGCGTCGGTTTCACGCAGTGCTGGCCTACAGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATATTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCCGCGTCACCATCGGAATCGGCGCGCCCCAGTCACGcc                                                                                                                          >  2:433332/1‑250 (MQ=255)                                                                                                                          tACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCCCCGGCAGCTGGTCGTCCGGGGTCGGCCTTACGGTACCGTCTTCCAGCGTCCCCATCGGACTCGGCGCGCCCCAGTAACGCTg                                                                                                                         >  2:336140/1‑251 (MQ=255)                                                                                                                           aCTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTg                                                                                                                         <  1:34396/250‑1 (MQ=255)                                                                                                                            cTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGAc                                                                                                                       >  2:392570/1‑251 (MQ=255)                                                                                                                               cGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACggg                                                                                                                    >  1:243140/1‑251 (MQ=255)                                                                                                                                 ggCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCTCCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGACCCAGTAACGCTGACGGGaa                                                                                                                  <  1:365491/251‑1 (MQ=255)                                                                                                                                   cAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAac                                                                                                                <  2:160318/251‑1 (MQ=255)                                                                                                                                    aaGTGTAGCGCGCATAGTACCAGGATTACTCCATAAAGGTGTAGAAAGTGTCGGTTTCACGCAGTGTTGGCATACCGTTAACGGTCGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCCGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTACCCATCGGAATCTGCGCGCCCCAGTAACGCTGACGGGAAacc                                                                                                               >  1:462504/1‑250 (MQ=255)                                                                                                                                      gtgtAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACTCTGACTGGAAACAcc                                                                                                             <  2:161373/251‑1 (MQ=255)                                                                                                                                      gtgtAGCGCGCATAGTACCAGGAGGACTCCATAAAGGGGTCGAAAGTGTCGGTTTCGCGCAGTGCTGGCATACCGTTAGCGGTAGTTTTCGCCCACCCCGGATCTGCTTTAATAGGGCTGGTAATTCCTCCCATTACCGCATCTTCCGTCAGCATCACCGGCAGCGGGTCGTCCTGGGGCGGCATTACGGTACCGTCTTCCAGCGCCACCATCCGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACAcc                                                                                                             <  1:212707/251‑1 (MQ=255)                                                                                                                                        gtAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGCAAGTGTCGGTTTCACGCAGTGCTGGCCTACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACAcccc                                                                                                           >  1:305445/1‑251 (MQ=255)                                                                                                                                        gtAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAACGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGAAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCCCCATCGGAATCGGCGCGCCCCAGTAACGCTGACCGGAAACAcccc                                                                                                           >  2:53349/1‑251 (MQ=255)                                                                                                                                        gtAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACAcccc                                                                                                           >  1:269206/1‑251 (MQ=255)                                                                                                                                                  tAGTACCAGGCGGACTCCATAAAGGCGTCGCAAGTGTCGGTTTCACGCAGAGCTGGCACACCGTTAACGGTAGTTTTCGCCCACTCCGGCTCTGCTTTAATAGAGCTGGTAATGCAGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGCCTTCCAGCGTCACCACCGGAATCGGCGCGCCCCCGCAACGCTGACGGGAAACCCCCCAGTCgcgc                                                                                                   >  2:198134/1‑249 (MQ=255)                                                                                                                                                         aggaggACTCCATAAAGGTGTCGAAAGGGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGCTCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCCGCGTCACCATCGGAATCGGCGCGCCCCAGTCACGCTGACGGGCAACACCCCAGTCGCGCAGGCggggg                                                                                          >  2:663/1‑248 (MQ=255)                                                                                                                                                           gaggaCTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGtt                                                                                        <  2:425054/251‑1 (MQ=255)                                                                                                                                                              gaCTCCATAAAGGTGTCGAAAGTGTCGGTTTCTCGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCAc                                                                                     <  2:265153/251‑1 (MQ=255)                                                                                                                                                               aCTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGGAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACCCCCCAGTCGCGCAGGCGTTAGTTCACt                                                                                    <  2:456204/251‑1 (MQ=255)                                                                                                                                                                    aTAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATGCCGTTAACGGTAGTTTTCGCCCACTCCGGAGCTGCTTTAATATGGCTGGTAATGCCGTCCATTACCTCATCTTCCGGCAGGATCACCGGCAGCTGGTCGGCCGGTGTCGGCAGTACGGTGCCGTCTTCCAGCGTCACCATCGGAATGTGCGCGCCCGAGTTACGTTTACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACg                                                                               <  1:215948/251‑1 (MQ=255)                                                                                                                                                                    aTAAAGGTGTCGAAAGAGTCGGTTTAACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGCATCGGCGCGCCCCAGTAACGCTGAAGGGAAACACCCCCGTCGCGCAGGCgggagctccatttacc                                                                               >  2:258753/1‑237 (MQ=255)                                                                                                                                                                      aaaGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGGTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCGGACGTGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGc                                                                              <  2:240243/250‑1 (MQ=255)                                                                                                                                                                       aaGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTc                                                                            <  1:186436/251‑1 (MQ=255)                                                                                                                                                                              cGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACg                                                                     <  1:129848/251‑1 (MQ=255)                                                                                                                                                                              cGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACg                                                                     <  2:243140/251‑1 (MQ=255)                                                                                                                                                                                     gtCGGTTTCACGCAGTTCTGGCATACCGTTAACGGTAGTTGTCGCCCACTCCGGATCTGCTTTAAGAGGGCTAGTAAGGCCGTCCATTACCACATCTGCCGGCAGGATCACCGGCAGCTGGTCGTCCGGAGTCGGCATTACGGTACCGTCTTCCAGCTTCACCATCGGAATCGGCGCTCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATc                                                               <  1:487831/250‑1 (MQ=255)                                                                                                                                                                                     gtCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGGTTTCGCCCACTCCGGATCTGCTTGAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTTACGGGAAGCACCCCAGTCGCGCAGGCGGTAGTTCACTTTGCGCTCGCCAACGCCCATc                                                               <  2:61868/250‑1 (MQ=255)                                                                                                                                                                                            tCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGcagtcn                                                        <  2:366090/250‑2 (MQ=255)                                                                                                                                                                                             cACGCAGTGCTGGCATACCGTTAACGGTCGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTAAGCCCGCCCACGCCCATCGcagtcagt                                                      >  2:39487/1‑251 (MQ=255)                                                                                                                                                                                                  aGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGGCAGGTTATc                                                 >  1:8663/1‑251 (MQ=255)                                                                                                                                                                                                           aTACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACCCCCATCGCCGTCAGTTTATCGGCGATggg                                        >  2:390651/1‑250 (MQ=255)                                                                                                                                                                                                           aTACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCGGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGc                                        <  2:409850/251‑1 (MQ=255)                                                                                                                                                                                                                                   cACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGTCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGTCAAg                <  1:435118/251‑1 (MQ=255)                                                                                                                                                                                                                                   cACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGgc                                                                                                                                          <  2:250163/129‑1 (MQ=255)                                                                                                                                                                                                                                   cACTCCGGATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGgc                                                                                                                                          >  1:250163/1‑129 (MQ=255)                                                                                                                                                                                                                                         ggATCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCCGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCCGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGTCACGACCGtt          >  1:297325/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCTGCTTTAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTTTCGGCGATGGCGTTGAAGGCCGCTTCATGGTCAAGACCGTTg         >  1:5992/1‑250 (MQ=255)                                                                                                                                                                                                                                                  gtAATAGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGGAGGGCACTTTACGCGCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGTCAAGACCGTTGAACTCGc  <  1:198134/249‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     CATCAGTTTGTGGAAGAAGCGGAAGTAGAGCAGGTGCATAATGGCGTGTTCAATACCACCAATGTAGATATCCACCGGCAGCCAGTAGTTAGCCGCTTCGGAATCCAGCATACCTTCTTTGTACTGCGGGCAAGTGTAGCGCGCATAGTACCAGGAGGACTCCATAAAGGTGTCGAAAGTGTCGGTTTCACGCAGTGCTGGCATACCGTTAACGGTAGTTTTCGCCCACTCCGGATCTGCTTTAATCGGGCTGGTAATGCCGTCCATTACCACATCTTCCGGCAGGATCACCGGCAGCTGGTCGTCCGGGGTCGGCATTACGGTACCGTCTTCCAGCGTCACCATCGGAATCGGCGCGCCCCAGTAACGCTGACGGGAAACACCCCAGTCGCGCAGGCGGTAGTTCACTTTACGCTCGCCAACGCCCATCGCAGTCAGTTTATCGGCGATGGCGTTGAAGGCCGCTTCATGGTCAAGACCGTTGAACTCGC  >  NC_000913/673140‑673630

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