Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	24,935	T→G	S849A (TCG→GCG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	24,935	0	T	G	98.1%	168.6 / ‑4.5	53	coding (2545/2817 nt)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (28/24);  minor base . (0/1);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCC  >  NC_000913/24688‑25177                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    ccTCTCCTTCACAGCTGATGAAGTGAGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCgg                                                                                                                                                                                                                                                  >  2:90616/1‑251 (MQ=255)                    gaAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTgg                                                                                                                                                                                                                                <  1:204278/250‑1 (MQ=255)                                 aCCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACg                                                                                                                                                                                                                   >  1:244445/1‑250 (MQ=255)                                         ggCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTGGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTTCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGTTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTctt                                                                                                                                                                                                          <  2:97374/251‑3 (MQ=255)                                           cGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTGGTCTTGt                                                                                                                                                                                                        <  2:157517/251‑1 (MQ=255)                                              aCGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGa                                                                                                                                                                                                    <  1:77281/251‑1 (MQ=255)                                                gTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGAcc                                                                                                                                                                                                  <  1:82548/251‑1 (MQ=255)                                                 tgcaaaaTACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGCTTTGTCTTGT‑TGGCCt                                                                                                                                                                                                 >  2:3093/4‑251 (MQ=255)                                                 tGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAACGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCt                                                                                                                                                                                                 >  2:17340/1‑251 (MQ=255)                                                   aaaaaTACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCCGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTACCCTTGTATGCAGAACCGGAAC‑TGGCGGCGCAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TTACCTcc                                                                                                                                                                                               >  1:130622/1‑251 (MQ=255)                                                     aaaTACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCgg                                                                                                                                                                                             >  2:92386/1‑251 (MQ=255)                                                      actACGTCTTCACCGGTGAGTGGGACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGAGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACAGCGCTGGGCGATGAATTACGATTTGTC‑TGTGTGTCCTCCGgc                                                                                                                                                                                            >  2:6859/3‑251 (MQ=255)                                                      aaTACGTCTTCACCGGTGAGTGGTACGAAGGACTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGAGATGAATTACGATTTGTCTTGT‑TGACCTCCGgc                                                                                                                                                                                            >  1:74642/1‑251 (MQ=255)                                                      aaTACGTCTGCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGgc                                                                                                                                                                                            <  2:98410/251‑1 (MQ=255)                                                        tACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCGGGGCGATGAATTACGATTTGTCTTGT‑TGCCCTCCGgcgc                                                                                                                                                                                          >  1:166630/1‑251 (MQ=255)                                                              ttCACCGGTGAGTGGTACGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGt                                                                                                                                                                                    >  1:78463/1‑251 (MQ=255)                                                                    ggTGAGTGGTACGAAGGCCTGTTTGGTCTGGCCGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGCAAGTGCGCGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGc                                                                                                                                                                                 >  2:194445/1‑248 (MQ=255)                                                                               cGAAGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGGCCTCCGGCGCTACCGTTGGAGAGTATAacgacg                                                                                                                                                                   <  2:128181/251‑1 (MQ=255)                                                                                  aGGCCTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCAc                                                                                                                                                                <  2:218829/251‑1 (MQ=255)                                                                                     ccTGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACctg                                                                                                                                                             <  1:120392/251‑1 (MQ=255)                                                                                       tGTTTGGTCTGGCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGACGTGAACAAAGTCATTGAGCACGCGCGTGCCGACAAGACAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGCTGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACctgct                                                                                                                                                           >  1:35913/1‑251 (MQ=255)                                                                                                 ggCAGACAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCAGTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcag                                                                                                                                                  <  2:195948/250‑1 (MQ=255)                                                                                                      aCAGTGAAGCAATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCAAAGTGGACAAAGTAATTGAGCCAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGc                                                                                                                                            >  1:39156/1‑251 (MQ=255)                                                                                                                aaTGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTc                                                                                                                                   >  1:5508/2‑250 (MQ=255)                                                                                                                    aaCGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                                                              <  2:68074/251‑1 (MQ=255)                                                                                                                      cgatgctgtCTGGGACGAGCTGTGGAAAGTGCGTGGCGAAGTGAACAAAGTCAGTGAGCAAGGGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGc                                                                                                                            <  2:173655/243‑1 (MQ=255)                                                                                                                            gTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGaa                                                                                                                        >  1:215880/1‑249 (MQ=255)                                                                                                                                 gggACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTcgcg                                                                                                                 >  1:81710/1‑251 (MQ=255)                                                                                                                                        gCTGTTGAAATTGCGGGGCGAAGTAAACAAAGTCATTGAGCAAGCGCGTGCCTGCAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTa                                                                                                          <  2:80817/251‑1 (MQ=255)                                                                                                                                           gTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaa                                                                                                         >  2:56455/1‑249 (MQ=255)                                                                                                                                                    ggggggggaaggtgaCAAAGTCAATGGGCAAGCGGGTGCCCGCAAGAAAGGGGGGGGGTCGGTGGGAGCGGCAGTAACCTTGTATGCAGGACCGGAACCTG‑CTGCGAAACTGGCCGCGCTGGGCCATGAATTACGGTTTGTCCTTT‑TGACCCCCGGCGCTTCCGTTGCAGACCATAACGACGCACCTTCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGCCGCGTTGAGTAAAGCCGAAGGTg                                                                                              <  2:168402/237‑1 (MQ=255)                                                                                                                                                     cggggCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTga                                                                                             <  2:173175/248‑1 (MQ=255)                                                                                                                                                      gTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACctgct                                                                                                                                                           >  1:59170/1‑188 (MQ=255)                                                                                                                                                      gTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACctgct                                                                                                                                                           <  2:59170/188‑1 (MQ=255)                                                                                                                                                      gTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTg                                                                                              >  1:243560/1‑249 (MQ=255)                                                                                                                                                           gAAGAGAACAAAGACATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAAACGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGCTGCATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAAGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTg                                                                                       >  2:138854/1‑251 (MQ=255)                                                                                                                                                            aaGTGAACAAAGTCAGTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGc                                                                                      <  2:118185/251‑1 (MQ=255)                                                                                                                                                             aGTGAACAAAGGCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCGGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGcc                                                                                     <  2:123695/251‑1 (MQ=255)                                                                                                                                                              gTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGccc                                                                                    >  1:48375/1‑251 (MQ=255)                                                                                                                                                                   aaaaGGCATTGAGCAAGCGCGGGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTGGTATGCAGAACCGGAAC‑GGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCGTGT‑GGACCTCCGGCGAGAACGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACGCAAAGGGCTGAAAGGCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgct                                                                               <  2:35913/250‑1 (MQ=255)                                                                                                                                                                        tCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTgg                                                                          >  2:70067/1‑251 (MQ=255)                                                                                                                                                                                 caGCGCATGCCGAATAGCAAGTGGGTGGCTGGCTGGACGCGGCAGTCACCTTGTATGCAGAACCGGCAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATCACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCCGAGCGACGTACTCCCACCGCTGAAAGTCGCGTTGAGTACCGCCGAAGTGGGGACGTGCCCAAGCTGCTGGCGCTacccc                                                                 >  2:44687/2‑248 (MQ=255)                                                                                                                                                                                         gCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATg                                                          <  1:92386/250‑1 (MQ=255)                                                                                                                                                                                                                  tGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGa                                 >  2:76020/1‑250 (MQ=255)                                                                                                                                                                                                                   ggAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGaa                                >  1:40707/1‑250 (MQ=255)                                                                                                                                                                                                                      aGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCt                            <  1:235724/251‑1 (MQ=255)                                                                                                                                                                                                                         ggCAGTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCg                         <  2:5508/251‑1 (MQ=255)                                                                                                                                                                                                                             gTAACCTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCg                     <  1:46458/251‑1 (MQ=255)                                                                                                                                                                                                                                  cTTGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTACAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGAGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACACAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtgtg                >  1:163003/1‑251 (MQ=255)                                                                                                                                                                                                                                    tGTATGCAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGCTGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTc              >  2:205037/1‑251 (MQ=255)                                                                                                                                                                                                                                          cAGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAAc        <  1:102205/251‑1 (MQ=255)                                                                                                                                                                                                                                           aGAACCGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACg       <  1:56455/251‑1 (MQ=255)                                                                                                                                                                                                                                                cGGAAC‑TGGCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCTTGT‑TGACCTCCGGCGCTACCGTTGAAGACTATAACGACGCACCTGCTGCTGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGACGCGTTGAGTAAAGCCGGAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGGAGGCAAGGAGGCGGAACACGAAGAAATCTGCGGCCGCTGTGGCAGCAACGTCgac  >  1:267743/1‑249 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    CCTCTCCTTCACCGCTGATGAAGTGTGGGGCTACCTGCCGGGCGAACGTGAAAAATACGTCTTCACCGGTGAGTGGTACGAAGGCCTGTTTGGCCTGGCAGACAGTGAAGCGATGAACGATGCGTTCTGGGACGAGCTGTTGAAAGTGCGTGGCGAAGTGAACAAAGTCATTGAGCAAGCGCGTGCCGACAAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAAC‑TGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGT‑TGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCC  >  NC_000913/24688‑25177

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