Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,204,890	C→T	G348S (GGT→AGT)	purD ←	phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,204,890	0	C	T	96.2%	171.0 / NA	53	G348S (GGT→AGT)	purD	phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (30/20);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTAGAGAAAAC  >  NC_000913/4204643‑4205129                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                tAGTTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GCTATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGAACACGTCACTGGt                                                                                                                                                                                                                                              >  2:339126/1‑251 (MQ=255)     tCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTAACTACCTGATCGTCATCCGCCAGTTTTGCGCCCGAGTGGACCACTTTGCCGCCTGCCACTTCTTCCAGCGGAAGGCCGTGGATCACGTCACTGGTGCgg                                                                                                                                                                                                                                          >  2:44666/1‑251 (MQ=255)        gCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTa                                                                                                                                                                                                                                        >  1:37644/1‑250 (MQ=255)           cgcttgcgataggcggacAGCCGATATTTCCCCGGCAGAAGCAGTTGCCCCGGTG‑GCTATCGGTCATTAAGGCATGGTGGCTTTTCTG‑GCTTTCTGCCACGG‑TTGGCCCCAGGGCGGGGACTCAGAGTCCGCGCCTGCCGTTGGTCATTACCTGCTCTTCATCCGCCAGTTTGGTGCCCGCGTGGAACACTTGGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAAt                                                                                                                                                                                                                                      <  1:207598/232‑1 (MQ=255)            gcgTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCAc                                                                                                                                                                                                                                   >  1:254916/1‑251 (MQ=255)                          gcCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GCTATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGGGCGGTAATCCCCCGGATATCCAccc                                                                                                                                                                                                                     >  2:21382/1‑251 (MQ=255)                            cAGCCGATATCTTCCGGGCAGAAGCAGTCGTCCCAGGG‑AATATCGGTAATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCCGGG‑TATGCCCCAGCGCGGGGCCGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGc                                                                                                                                                                                                                   <  1:154042/251‑1 (MQ=255)                              gCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TCTGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGGCACTGGTGCGGTAATCCCCCGGATTACCACCCGccg                                                                                                                                                                                                                 >  2:155253/1‑249 (MQ=255)                                cGATATCTTCCCGGAAGAGGAGGTCGTCCCAGTG‑GATACCGGTCATAAAGGCATAGGGGCTTTTCTGCCTTT‑CTGCCCCGG‑TATGCCCGGGCGCGGTGACCCACAGTCCGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGcc                                                                                                                                                                                                               <  1:236740/251‑1 (MQ=255)                                 attttCTTCCCGGCAGAAGCGGTCGCCCCGGGG‑GATATCGGTCATTAGGGCATAGGGGCGTTTCTGCCCTT‑CCCCCACGG‑TATGCCCCGGCGCGGTGCCGCCCGGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCa                                                                                                                                                                                                              <  1:268441/247‑1 (MQ=255)                                  atatCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAACGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑CATGACCCAGCGCGGTGACGCACAGTAAGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGCTCACGTCCCCGGTGCGGCAATCACCCGGATATCCACCCGCAGCCAt                                                                                                                                                                                                             >  2:368443/1‑251 (MQ=255)                                   tatCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATc                                                                                                                                                                                                            >  2:351283/1‑251 (MQ=255)                                          cgggaagaagcggtcgcccagggg‑gATATCGGTAATTAAGGCATAGGCGTGTTGGGGCGCTT‑C‑GCCCCGGGTAGGCCCAAGCGGGGGGACGCACTGCCCGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGc                                                                                                                                                                                                     <  1:342726/228‑1 (MQ=255)                                           cGGCAGAAGAAGTCGTCC‑AGTGATATATCGTGCATTAGGGCATGGGCCCGTTCCTGCGCTT‑TGGCCACGG‑TAGGACCCGGAGGGGTAACCCAAAGTTCCCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGcc                                                                                                                                                                                                    <  1:44666/251‑1 (MQ=255)                                              cAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCg                                                                                                                                                                                                   >  1:302850/1‑249 (MQ=255)                                               aGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCgagg                                                                                                                                                                                                >  1:326589/1‑250 (MQ=255)                                                  aGCAGACGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑ATTGACCCAGCCCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTCCCTGCTCGACATCCGCCAGTTTTGTGCACGCGCGGTACACTTTGCCGCCTGCCACTTCTTCCATCGGCAGTCACTGGAGCACGTCACACGTCCGGTAATTAACCCGATTCCCACCCACGGCCATCACCACGCCgcgaga                                                                                                                                                                                              >  2:335652/1‑250 (MQ=255)                                                          ccccAGTG‑GATATCGGTCATTAAGGCATGGGCGTGTTTCTGCGCTT‑CTGCCCCGT‑TATGCCCCAGCGCGGTGCCGCCCAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCAGCACCACGCCGAGAGAAGCGCGTTc                                                                                                                                                                                     <  1:154819/250‑1 (MQ=255)                                                                      atCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGCGGATCACGTCACTGGGGCGGTAATAACCCGGTTATCCACCAGCAGCCATCACCACCCCGCAAGAAGCGCGTTCCTCCCCCTCgg                                                                                                                                                                          >  2:329769/1‑251 (MQ=255)                                                                          gTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTCCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGAACACGTCACTGGGGCGGTAATCACCCGGAAATCCACCCGCAGCCATCCCCACGCCGAGAGAAGCGCCTTCATCCCACTCGGACGt                                                                                                                                                                      >  2:238075/1‑251 (MQ=255)                                                                             aTTAAGGAATAGGAGCGTTTCTGCGCTT‑CTGCAACGG‑TATGACCCAGAGCGGTAACGCACCGTACGCGACCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGGACACTTTGCCGCCTGCCACTTCTTACAGCGGCGGGCCGTGGATCACGTCACTGGGGCGGTAATCACCAGGCAATCCACCCGAAGACATCACCAGGCCGAGAGAGGCGAGTTCAtacccagcgcacgtttt                                                                                                                                                                   >  2:31383/1‑235 (MQ=255)                                                                                aaGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTtctc                                                                                                                                                                >  2:126867/1‑251 (MQ=255)                                                                                       aGGCGCGTTTCTGCGCTC‑CCGC‑ACGGGTATGCCCCGGCGGGGGGACGCACAGTCCCCTCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGt                                                                                                                                                         <  1:433422/251‑1 (MQ=255)                                                                                        ggCGCGTTTCTGCCCTT‑CTGCCACAG‑TCTGACCCAGCGCGGTGACACAACGCACGCGCCCGACGTTGGTAACTCCCTGCTCGTCATCCGCCAGTTTTCTGCCCGCGTGGAACACTTTGCCGACTGCCACTTCTTCCAGCAGCAGGCCGAGGCTCAAGCCACTGGTCCGGTAATCCCCCGGATATCCACCCGGAGCCATCACCACGCCGAGAAAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCactt                                                                                                                                                        >  2:424093/1‑248 (MQ=255)                                                                                             gTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCtt                                                                                                                                                   >  1:236887/1‑251 (MQ=255)                                                                                                     gcTT‑CTGCCACGG‑TATGACACAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGCTCACGTCACCGGTGCGGTAATCACCCGGCTATCCACCCGCAGCCATCCCCACGCCGAGAAAAGCCCGTTCATCCCACCCGGACGTTTTCTAGTCCAGTTTGCTTTCACAGgc                                                                                                                                           >  2:44185/1‑251 (MQ=255)                                                                                                           tGCCACGG‑TCTGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCCCTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTACTCACCCGGCTATCCACCCGCAGCCATCACCACGCCGCGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTGTGCTTTCACCGgccgccc                                                                                                                                      >  2:285672/1‑250 (MQ=255)                                                                                                                cGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCAcc                                                                                                                                                                                                                      >  1:113660/1‑166 (MQ=255)                                                                                                                      ggccccAGCGGGGTGACGCAGAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTc                                                                                                                            <  1:90291/248‑1 (MQ=255)                                                                                                                                  gtgccgcccattACGCCCCCGCCGTTGGTCACTCCCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGa                                                                                                                <  1:431974/240‑1 (MQ=255)                                                                                                                                       gCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTc                                                                                                            >  1:253844/1‑250 (MQ=255)                                                                                                                                            gTACGCGCCAGCCGTTGGTCACTACATGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCCCCCGCAGCCATCACCACGCCGAGAGCAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCCGATCCCAACCAGATCGGACTTCCTgcgc                                                                                                      >  2:63263/1‑251 (MQ=255)                                                                                                                                               cgcccccgccGTGGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCgcagc                                                                                                   <  1:351283/247‑1 (MQ=255)                                                                                                                                                   cccgccGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATaa                                                                                               <  1:326988/251‑1 (MQ=255)                                                                                                                                                         gTTGGTAACTACCTGCTCGTCATCCGACAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAACGCGTTCATCCCACCCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGAACGGACTTCATGCGCAGCATAATCGGcg                                                                                         >  2:178411/1‑250 (MQ=255)                                                                                                                                                               cACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCCGTTTGCTTTCACAGGCCGCCAGGCAGCGCTCAACCAGACCGGCCCTCACGCGCCGCATAATCGGCTGGGttt                                                                                   >  2:9300/1‑251 (MQ=255)                                                                                                                                                                aCTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTc                                                                                  >  2:304983/1‑251 (MQ=255)                                                                                                                                                                   accgccTCGCCATCCGCCGGTTTTGTGCCCGGGTGGAACGCTTTGCCGCCGGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGGGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTGCACAGGCCGCCAGGCAGAGCGCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTcc                                                                                 <  2:189835/244‑1 (MQ=255)                                                                                                                                                                    ccTGCTCGTAATCCGCAAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGCCCCTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGCTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCCGGCAGGGCTCAACCAGATCGGACTTCATGCGCAGCCTAATCGGCTGGGTTTcccgg                                                                              >  2:27078/1‑248 (MQ=255)                                                                                                                                                                               cccGCCAGTTGTGTGCCCGCGTGGAACACTTAGCCGCCTGCCACTTCTTCCAGCTGCAGGCCGTGGATCACGTCACTGGTGCGGTAAGCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTGCAGCCCACTCGGACGTTGTCGCGGCCAGTTTGCTGTCACAGGCCGCCAGGCAGAGCTCAACCAGATCTGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCg                                                                   <  1:424093/250‑1 (MQ=255)                                                                                                                                                                                ccGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCgg                                                                  <  1:285672/251‑1 (MQ=255)                                                                                                                                                                                  gCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGc                                                                 >  1:23170/1‑250 (MQ=255)                                                                                                                                                                                  gCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGCTTCCGGATCGCCAAAGCGGc                                                                 >  1:417234/1‑250 (MQ=255)                                                                                                                                                                                     aGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGtt                                                             <  1:21382/251‑1 (MQ=255)                                                                                                                                                                                     aGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTAAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGtt                                                             <  1:126867/251‑1 (MQ=255)                                                                                                                                                                                                                     tGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGt                             >  1:342088/1‑251 (MQ=255)                                                                                                                                                                                                                      gCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTACTCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAACGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGtt                            >  1:225867/1‑250 (MQ=255)                                                                                                                                                                                                                      gCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTc                            >  2:49921/1‑251 (MQ=255)                                                                                                                                                                                                                                  aGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATGCAATAACCTTCGTATTGCCCTGTTTGTCGATCATCAGGcc                <  1:367215/251‑1 (MQ=255)                                                                                                                                                                                                                                  aGCGGCAGGCCGTGGATCACGTCACTGGTGCGGTAAACACCCGGATATCCACCCGCAGCCATCACCCCGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCCAAGCGGCAGTTAAATTCAATACCCTTCGGATTGCCCCGTTTGTCGATCATCAGGcc                >  2:162185/1‑251 (MQ=255)                                                                                                                                                                                                                                    cGGCAGGCCGTGGATCACGACACTGGGGCGGTAATCACCCGGATATCCACCCGCAGCCATCCCCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGGCGATCATCAGGCCcg              <  2:207766/251‑1 (MQ=255)                                                                                                                                                                                                                                            cGTGGATCACTTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTAgaga      <  2:137515/251‑1 (MQ=255)                                                                                                                                                                                                                                                gATCACGTCACTGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCCAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTAGAGAAAAc  >  1:362252/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                TAGTTCTGCTCGCGTTCGATAGCGCGCCAGCCGATATCTTTCCGGCAGAAGCAGTCGTCCCAGTG‑GATATCGGTCATTAAGGCATAGGCGCGTTTCTGCGCTT‑CTGCCACGG‑TATGACCCAGCGCGGTGACGCACAGTACGCGCCCGCCGTTGGTCACTACCTGCTCGTCATCCGCCAGTTTTGTGCCCGCGTGGAACACTTTGCCGCCTGCCACTTCTTCCAGCGGCAGGCCGTGGATCACGTCACCGGTGCGGTAATCACCCGGATATCCACCCGCAGCCATCACCACGCCGAGAGAAGCGCGTTCATCCCACTCGGACGTTTTCTCGTCCAGTTTGCTTTCACAGGCCGCCAGGCAGAGCTCAACCAGATCGGACTTCATGCGCAGCATAATCGGCTGGGTTTCCGGATCGCCAAAGCGGCAGTTAAATTCAATAACCTTCGGATTGCCCTGTTTGTCGATCATCAGGCCCGCGTAGAGAAAAC  >  NC_000913/4204643‑4205129

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