Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	24,646	A→G	L752L (CTA→CTG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	24,646	0	A	G	98.3%	187.5 / NA	60	L752L (CTA→CTG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (28/31);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCAAAAGATATGGTGAAACCGGAAGAGATGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGTGTGGCGCGTCGTAGCTGCCAGACTGCGCTATATCACATCGCAGAAGCGCTGGTGCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTG  >  NC_000913/24404‑24890                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     gCAAAAGATATGGTGAAACCGGAAGAGATGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACAt                                                                                                                                                                                                                                              >  1:157517/1‑251 (MQ=255)        atatGGTTAAACCGGAAGAGATGGGGGTACTGGATCGCTGGGCCGTAGGTTTTGCGAAAGCGGCACAGGAAGTCATCCTCAAGGCGGACGAAGCATACGATTTCCACGAAGTTGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACTTCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGaa                                                                                                                                                                                                                                       <  2:232443/251‑1 (MQ=255)         tatGGTGAAACCGGAAGAGATGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGa                                                                                                                                                                                                                                        >  1:173655/1‑249 (MQ=255)                aaaCCGGAAGAGATGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGt                                                                                                                                                                                                                               >  1:68074/1‑251 (MQ=255)                      gAAGAGATGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTg                                                                                                                                                                                                                         <  1:94199/251‑1 (MQ=255)                       aaGAGATGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCGGAAGCGCTGGATCGCTgg                                                                                                                                                                                                                        >  1:97374/1‑251 (MQ=255)                          agaTGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTgg                                                                                                                                                                                                                                >  1:128620/1‑240 (MQ=255)                          agaTGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACAGCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGGCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTgg                                                                                                                                                                                                                                <  2:128620/240‑1 (MQ=255)                           gatgggggtaatgggtagatGGGCCGTAGATTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGCAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACATCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGatggatgg                                                                                                                                                                                                                    >  1:83265/20‑251 (MQ=255)                            aTGGTGGTACTGGACCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGc                                                                                                                                                                                                                   >  2:216695/1‑251 (MQ=255)                              ggtggtACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAATACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTTCAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCAc                                                                                                                                                                                                                 <  1:245130/251‑1 (MQ=255)                               gtggtACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGTCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCAcc                                                                                                                                                                                                                <  1:5514/251‑1 (MQ=255)                                  gtACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAAt                                                                                                                                                                                                             >  1:187648/1‑251 (MQ=255)                                           cGCTGGGCCGAAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCtt                                                                                                                                                                                                    >  2:235724/1‑251 (MQ=255)                                               gggCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCa                                                                                                                                                                                                  <  1:104513/249‑1 (MQ=255)                                                ggCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCg                                                                                                                                                                                               <  2:217603/251‑1 (MQ=255)                                                        gTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCCCGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGGTCGCTGGATGGCACCAATCCTCTCCTTCACCGCtgatga                                                                                                                                                                                        >  2:82548/1‑250 (MQ=255)                                                             gCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAgtgt                                                                                                                                                                                   <  2:73975/250‑1 (MQ=255)                                                                  aGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGCGATGGGTTCCTTCTACCTCGACATCATCACAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCt                                                                                                                                                                             >  2:229617/1‑251 (MQ=255)                                                                      gCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCTCCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCt                                                                                                                                                                         <  2:220767/251‑1 (MQ=255)                                                                                  aTCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTAACTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACg                                                                                                                                                             >  1:168402/1‑251 (MQ=255)                                                                                       cAAGGCGTACTAAGCTTACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGaa                                                                                                                                                         <  1:119238/250‑1 (MQ=255)                                                                                             gTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACg                                                                                                                                                  >  2:50979/1‑251 (MQ=255)                                                                                               aCGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGt                                                                                                                                                 <  1:251070/250‑1 (MQ=255)                                                                                                cGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCt                                                                                                                                               <  2:187648/251‑1 (MQ=255)                                                                                                 gAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCtt                                                                                                                                              >  1:98410/1‑251 (MQ=255)                                                                                                        aCGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCCAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTACGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCggg                                                                                                                                       >  2:204278/1‑250 (MQ=255)                                                                                                                        ggTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGcc                                                                                                                       >  1:118185/1‑251 (MQ=255)                                                                                                                          tACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTAGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTg                                                                                                                     >  1:128181/1‑251 (MQ=255)                                                                                                                                gTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCCAAGCGGACAGAGTGGCGCGTCGTAGCTGCCAGACTGCGCTATATAACATCGCAGACGCGCTGGGCCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCAACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGAGAGGGGAACGAAGGCCTGTTTggt                                                                                                               >  1:45693/1‑250 (MQ=255)                                                                                                                                    gATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTCCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTggtct                                                                                                             <  1:20987/249‑4 (MQ=255)                                                                                                                                    gATGCGCTTCTGCTCCGTTGAGATGGGTTCCTGCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTg                                                                                                            <  2:73181/250‑1 (MQ=255)                                                                                                                                     aTGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGc                                                                                                          <  1:204329/251‑1 (MQ=255)                                                                                                                                     aTGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCGCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGc                                                                                                          <  1:184302/251‑1 (MQ=255)                                                                                                                                       gcgcTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAg                                                                                                        <  1:198033/251‑1 (MQ=255)                                                                                                                                         gcTTCTGCGCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGAc                                                                                                      <  2:71555/251‑1 (MQ=255)                                                                                                                                                         gaTGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACCCCGCCAAAGCGGACAGCGTCGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGGACGAAGGCCTGTTTGGTCTGGCAGACAGTGAATCAATGAACg                                                                                      <  1:211706/251‑1 (MQ=255)                                                                                                                                                         gaTGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCTCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCTCTGGTGCGCTGGAGGGCACCAATCCTCTCCTTCACCGCTGATGATGTGTGGGGCTACCTGCCGGGCGAACGGGAAAAATACGTCTTCACCGGTGATTGGTACGAAGGCCTGGTTGGTCTGGCAGACAGTGAAGCAATGAACg                                                                                      <  2:253903/251‑1 (MQ=255)                                                                                                                                                             ggTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGc                                                                                  >  1:146545/1‑251 (MQ=255)                                                                                                                                                                tCCTGCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGt                                                                                <  1:115295/250‑1 (MQ=255)                                                                                                                                                                            gACATTATCACAGACCGACAGTCCACCGCCAAAGCTGACCGCGTGTCGCGTCGTAGCTGCCAGACTGCGCTGTATCCCATCGCAGAAGTGCCGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGCAGTGTGGGGCTACCTGCCGGGCGTACGTGCAAAATACGTCTACACCGGTGAGTGGTACTAAGGCCTGATTGGTCTGGCAGACAGTGTAGAAATGAACGTAGCGTGCTGGGGCGAGCt                                                                   >  2:223585/1‑251 (MQ=255)                                                                                                                                                                            gACATCATCAAAGACGGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTACGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCt                                                                   >  2:223592/1‑251 (MQ=255)                                                                                                                                                                             aCATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTg                                                                  <  1:253349/251‑1 (MQ=255)                                                                                                                                                                                tcatcaAAGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGCTGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTg                                                               >  2:120392/1‑251 (MQ=255)                                                                                                                                                                                 catcaAAGACCGTCCGTACACCGCCAAAGCGTACAGCGTGGCGCGTCGTAGCTTCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGa                                                              <  2:142984/251‑1 (MQ=255)                                                                                                                                                                                     aaaGACCGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGt                                                          >  2:229491/1‑251 (MQ=255)                                                                                                                                                                                     aaaGACCGTCAGTAACCCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGt                                                          <  2:223866/251‑1 (MQ=255)                                                                                                                                                                                          ccGTCAGTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTg                                                     <  1:50979/251‑1 (MQ=255)                                                                                                                                                                                           cGTCATTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGGGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTgg                                                    <  1:20974/251‑1 (MQ=255)                                                                                                                                                                                                gTACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGAGCGTGGCGa                                                 >  2:71667/1‑249 (MQ=255)                                                                                                                                                                                                 tACACCGCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGt                                              >  1:80817/1‑251 (MQ=255)                                                                                                                                                                                                       gCCAAAGCGGACAGCGTGGCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACa                                         >  2:77281/1‑250 (MQ=255)                                                                                                                                                                                                                   agcgaggCGCGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGCGCGTGGCGAAGTGAACAAAGTCATTGAGcc                            >  2:102205/6‑250 (MQ=255)                                                                                                                                                                                                                         ggtcgacgtAGCTGCCAGACTGCGCCGCATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCACTCCTCTCCTGCACCGCAGATGAAGTGTGTGGCTACCTGCCGTGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATTCGTTCTGGGTCGATCTGTTGAAAGTGCGTGGCGAGGTGAACAAAGTCTTTGAGCAAGcgcg                      <  2:181359/245‑1 (MQ=255)                                                                                                                                                                                                                         gcgcGTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGcgcg                      <  1:121752/251‑1 (MQ=255)                                                                                                                                                                                                                           gcgTCGTAGCTGCCAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTg                    <  1:215523/251‑1 (MQ=255)                                                                                                                                                                                                                                       ccAGACTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGTCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAg        <  2:71311/251‑1 (MQ=255)                                                                                                                                                                                                                                          gACTGCGCTGTATCACATCGCAGAAGCGCTGCTTCGCTGGATGTCACCACTCCTCTCATTCACCGCTGATGACGTGTGGGGCTACCTGCCGGGTGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGAATGGTCTGGCAGACAGTGAAGCCATGAACGATGAGTTCTGGGACGAGCGGTTGAAAGTGCGTGGCGAGGTGTAAACAGTCAGTAAGCGAGCGCGTGCCAACAAGAAAgtgg     >  2:43672/1‑251 (MQ=255)                                                                                                                                                                                                                                            cTGCGCTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGGGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAgtgggg   >  1:12022/1‑250 (MQ=255)                                                                                                                                                                                                                                              gcgcTGTATCACATCGCAGAAGCGCTGGTTCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAgtgggtg  <  1:229491/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     GCAAAAGATATGGTGAAACCGGAAGAGATGGTGGTACTGGATCGCTGGGCCGTAGGTTGTGCGAAAGCGGCACAGGAAGACATCCTCAAGGCGTACGAAGCATACGATTTCCACGAAGTGGTACAGCGTCTGATGCGCTTCTGCTCCGTTGAGATGGGTTCCTTCTACCTCGACATCATCAAAGACCGTCAGTACACCGCCAAAGCGGACAGTGTGGCGCGTCGTAGCTGCCAGACTGCGCTATATCACATCGCAGAAGCGCTGGTGCGCTGGATGGCACCAATCCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTG  >  NC_000913/24404‑24890

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