Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,270	G→A	G294S (GGC→AGC)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,270	0	G	A	92.9%	109.5 / NA	42	G294S (GGC→AGC)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/21);  minor base C (2/1);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGG  >  NC_000913/23026‑23508                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                tGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGAGTACTCCATTCTCAGCACg                                                                                                                                                                                                                                            <  2:40864/250‑1 (MQ=255)  gCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAACGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCAAGGt                                                                                                                                                                                                                                          >  2:217639/1‑251 (MQ=255)                            aCTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGCGCTTGAGTtgctgc                                                                                                                                                                                                                >  1:110693/1‑251 (MQ=255)                            aCTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTtgctgc                                                                                                                                                                                                                <  1:114738/251‑1 (MQ=255)                               gAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGTTGCGGAGCTTGAGTTGCTgcg                                                                                                                                                                                                               <  2:227347/249‑1 (MQ=255)                                     aaaaTTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGGGCTTGAGTTGCTTCGCTTTAccc                                                                                                                                                                                                       >  2:8465/1‑251 (MQ=255)                                     aaaaTTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTAccc                                                                                                                                                                                                       <  2:157949/251‑1 (MQ=255)                                                          tAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGAc                                                                                                                                                                                   <  1:130852/250‑1 (MQ=255)                                                                        cGCTGGTGATCTGGACCACGACGCCGTGGACTCTGCCTGCGAACCGCGCAAACTCTATTGCACCTGATTGCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATGCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTc                                                                                                                                                                    <  2:41771/251‑1 (MQ=255)                                                                              tgaTCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGAt                                                                                                                                                              >  1:158707/3‑251 (MQ=255)                                                                                            acGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTg                                                                                                                                                  <  1:194089/249‑1 (MQ=255)                                                                                                 gTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCGGATGTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTCCCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCg                                                                                                                                           <  2:13409/251‑1 (MQ=255)                                                                                                   ggACTCTGCCTGCGAGCCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGCAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCCGCACGGTAAAAGGTGCGGGGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGGCGTTCCGGCAATCCTCGGCGATACCGTTACCCTGGATGCCGGt                                                                                                                                         >  2:43041/1‑251 (MQ=255)                                                                                                        ctgcctgcGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTAAACCATTCTCAGCACGGTGAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCgg                                                                                                                                    >  2:107463/1‑251 (MQ=255)                                                                                                                                   gCACCTGATTTCGACTCTGCGCTGGTGACGATAGACGGTCAGGCCGTGATTCTGGCGAACGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGCTTAACCCATTCTCAGCACGGTAAAAGGTGCGGCGCTTGCGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCCCGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACgg                                                                                                        >  2:16916/1‑251 (MQ=255)                                                                                                                                    catctgattgcggcTATGCGCTGGTGCCGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAATGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTTCCCTGGATGCCGGTACCTG‑TGCCGTTCACACCTCGCGTGGGCACGGc                                                                                                       <  2:41774/238‑1 (MQ=255)                                                                                                                                     acctgATTTCGACTATGCGCTGGTGCAGATCTAGGGGGAGGCCGTGTTTCTGAAGAAAGATCTGGTTGAAAGCGGAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAACGGTGCGGAGCTTAAGTTGCTGCGCCTTACCCATCCGTTTTTGGGCCTCGACAATCCCGCAATCCCCGGCGATCACGTTCACATGGACGCCGGTACCGGGTCCAGTCCACAC‑GCGCCTGGCCCCGGcc                                                                                                      >  2:232770/5‑251 (MQ=255)                                                                                                                                         gATTTCGACTATGCGCTGGGGCAGAGCGACGGGCAGGCCGTGACTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCAGTCTCCGCACTGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGTCATTCCTCGGCGATCCCGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGga                                                                                                  <  2:117698/251‑1 (MQ=255)                                                                                                                                                 cTATGCGCTGGTGCAGCTCGACGGTCAGGCCGTGCTTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCCTGAGTTGCTGCGCTTTACCCAACCGGTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCCCGTTACCCTGGATGCCGGTACCGG‑CTCCGTTACCACCGCGCCTGGTCACGGCCCGGCCGACTAtg                                                                                          >  2:54037/1‑251 (MQ=255)                                                                                                                                                   aTGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAACACGGTAAACGGCGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCCATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGGCTAtgtg                                                                                        >  1:158236/1‑251 (MQ=255)                                                                                                                                                               aGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGCAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGACCACCGCCCGGACGACTATGTGATCGGCCAGaa                                                                             >  1:217128/1‑250 (MQ=255)                                                                                                                                                                     aCGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCCCGGCCCGGCCGACTATGTGATCGGTCAGAAATACGGc                                                                      >  2:86159/1‑251 (MQ=255)                                                                                                                                                                       ggTCAGGCCGTGATTCTGGCGAAAGACCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTTAGTTGCTGCGCTTTCCCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑GGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCt                                                                    <  2:123539/251‑1 (MQ=255)                                                                                                                                                                        gTCAGGCCGTGATTCGGGGGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGGCCGATTACACCATTCTCATCACGGTAAAAGGGGCGGAGCTTGAGTTGCTGCTCGTTACCCATCCGGTTATGGGCTTCGACGTTGCGGCAATCCGCGGCGATCACGTTACCCTGTATGGCTGTGCCTG‑TGCCGTTCACATCGCGCCGGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTg                                                                   <  1:132767/251‑1 (MQ=255)                                                                                                                                                                              ccGTTATTCTGGCGCAAGATCTGGGTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACCCCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCATTCCTCGGCGATCACGTTACCCTGGATGCCGTTACCGG‑TGCCGTTCACACCGCGCCTGGTCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAAcc                                                             <  1:69794/251‑1 (MQ=255)                                                                                                                                                                                gTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGCCGTGCCCGATTACACCATTCTCAGAACGGTAAAAGGTGCGGAGCTTGAGTTGCTGAGCTTTCCCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTAACGG‑TGCCGTTCACACCGCGCCTCGCCACGGCACGGCCGACTATGTGCTCGGTCAGAAATACGCCCTGGAAACCGc                                                           >  1:166124/1‑251 (MQ=255)                                                                                                                                                                                 tGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAACATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCCCGTTACCCTGGATGCCGGGACCGG‑TGCCGTTCACACCGCGCCTGGACACGGCCCGGAAGACTATGTGATCGGTCAGAAATGCGGCCTGGAAACCGCt                                                          >  1:7072/1‑251 (MQ=255)                                                                                                                                                                                        ggCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCg                                                    <  1:8465/250‑1 (MQ=255)                                                                                                                                                                                          cGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACTCCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGtt                                                 <  2:79682/251‑1 (MQ=255)                                                                                                                                                                                          cGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGAGGAGATTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGGTACCCTGGATGCCGGTACCGG‑GGCCGTTCAGACCGCGCCTGGCCACGGCCCGGACGAATATGTGATCGGTCAGACAAACGGCCTGGAAACCGCTAACCCGGtt                                                 >  1:152526/1‑251 (MQ=255)                                                                                                                                                                                                   gggTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGAc                                        <  1:206090/250‑1 (MQ=255)                                                                                                                                                                                                     gTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACgg                                      >  1:175060/1‑251 (MQ=255)                                                                                                                                                                                                        taaaGCGTAATGCAGCGTATCGGCGTGTCCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCa                                    <  2:219497/249‑1 (MQ=255)                                                                                                                                                                                                        gAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCCGCACGGTAAAAGGGGCGGGGCTTGAGTTGCTGCGCTTTACCCACCCGTTTTTGGGCTTCGCCGTTCCGGCAATCCTCGGCAATCCCGTTACCCTGGTTGCCGGGCCCGA‑TGCCGTCAACACCGGGCCTGGCCCCGGCCCGGACGCCTACGTGTTCGGCCAGAACTACGGCCGGAAAACCGCCAACCCGGTTTGCCCGGACGGccc                                   >  2:234027/1‑249 (MQ=255)                                                                                                                                                                                                              gTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCt                             <  1:16916/251‑1 (MQ=255)                                                                                                                                                                                                                aagtcAGCGTATCGGCGGGACCTATTACACCATTCGCAGCACGGTAAAAGGTGCGGAGCTTTATTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCGGGCCACGGCCCGGACGACTATGGGATCGGTCAGAAATACGGCCTGGAAACCGCTATCCCGGTTGGCCCGTACTGCACTTATCTGc                           <  2:256872/247‑1 (MQ=255)                                                                                                                                                                                                                 aTGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGcc                          <  1:86159/251‑1 (MQ=255)                                                                                                                                                                                                                    cAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGATTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATAAGGCCTGGAAACCGCTAACCCGGTTGGCCCGGCCGGCACTTATCTGCCggg                       >  1:136360/1‑251 (MQ=255)                                                                                                                                                                                                                       cGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGGGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCAc                    >  2:86995/1‑251 (MQ=255)                                                                                                                                                                                                                        ggagcGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTCCCCATCCGTTTATGGGCTGCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCGGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACt                   <  2:110693/247‑1 (MQ=255)                                                                                                                                                                                                                              gCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCg             >  1:151956/1‑251 (MQ=255)                                                                                                                                                                                                                                          acacCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTGACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCTCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATgg  <  1:43041/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                TGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGG‑TGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGG  >  NC_000913/23026‑23508

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