Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_pae_tpiA	4,032,422	G→A	5.7%	R421H (CGT→CAT)	pepQ →	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	4,032,422	0	G	A	5.7%	256.4 / 4.9	88	R421H (CGT→CAT)	pepQ	proline dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (47/35);  new base A (3/2);  total (50/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

CCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTT  >  NC_000913_3_pae_tpiA/4032275‑4032568                                                                                                                                                    |                                                                                                                                                   cccGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGt                                                                                                                                                   >  1:787299/1‑149 (MQ=255)  ccGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTGGTa                                                                                                                                                  >  2:770102/1‑149 (MQ=255)  ccGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGACAGAAAATTGAAGCACTGAAACCGTTCGCCGGCATTCGTa                                                                                                                                                  >  2:334785/1‑149 (MQ=255)  ccGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCACTACAAGCACCGTGGGCTGAAGGGCAGTTCAGCAAGCACTTCAACGGGCAGAAAATTAAAGCACTGAAACCGGTCGGAGGCATACGTa                                                                                                                                                  >  2:203018/1‑149 (MQ=255)   cGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTAt                                                                                                                                                 <  2:765427/149‑1 (MQ=255)       ttCTCCAGCCGCGCATGGTGTTAACCATCGAACCGGGTATCTACTTCCTCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAACCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTAACGaa                                                                                                                                             >  2:355335/1‑149 (MQ=255)       tcctcCAGCCGGGCATGGTCTTAACCATCGAACCGGGTATCTACTCCATCGAATCGCGACTGGCACCGTGGCGTGAAAGGCAGTTCAGCAAGCATTTCAACTGTCATAATATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGaa                                                                                                                                             <  1:279621/147‑1 (MQ=255)         ctcCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTGCATCGAATCGCTACGGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGa                                                                                                                                           <  1:81987/149‑1 (MQ=255)          tcCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGGAGCACTGAAACGGTTCGGCGGCCTTCGTATCGacaac                                                                                                                                          >  1:197827/1‑145 (MQ=255)          tcCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTACCGAAggc                                                                                                                                          >  1:793771/1‑147 (MQ=255)                cGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACgtg                                                                                                                                    <  2:781176/149‑1 (MQ=255)                  gtacatgtgtTACCCTGCGAACCGGTTATTTTCTTTCTCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCCGCAAGCTCTTCAACTGGCAGAAAATTGAAGCACTGAACCCGTTCGGCGGCATTCGTATCGAAGACAACgtggt                                                                                                                                  <  2:123539/143‑1 (MQ=255)                  ggCATGGTGTTAACCATCGAACCGGGTATCTACCTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACgtggt                                                                                                                                  <  1:132361/149‑1 (MQ=255)                      tGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAGGACAACACGGAGATc                                                                                                                              >  2:337384/1‑149 (MQ=255)                      tGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATc                                                                                                                              >  2:818146/1‑149 (MQ=255)                        gtgtTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCa                                                                                                                            <  2:769134/149‑1 (MQ=255)                        gtgtTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCa                                                                                                                            <  1:774907/149‑1 (MQ=255)                         tgtTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCCTTCGGCTGCATTCGTATCGAAGACAACGTGGTGATCCAc                                                                                                                           <  2:902208/149‑1 (MQ=255)                          gtTAACCCTGGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGCGCAGTTCATCAACGACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACg                                                                                                                          <  2:104019/149‑1 (MQ=255)                            tAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAATTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGaa                                                                                                                        <  2:83634/149‑1 (MQ=255)                             aaCCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGaaa                                                                                                                       >  2:68103/1‑149 (MQ=255)                              aCCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATGCACGaaaa                                                                                                                      >  1:497551/1‑149 (MQ=255)                                cATCGAATCGTGTATCTCCCTAATCTAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCATAAAATTGAAGCACTGAAACCGTTCGGCGGCACTCGTATCGAAGACAACGTGGTGATCCACGAAaaca                                                                                                                    <  2:448240/149‑1 (MQ=255)                                cATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGCAACCGTTCGGCGGCATTCGTATCGACGACAACGTGGTGACCCACGaaaaaa                                                                                                                    >  2:373180/1‑147 (MQ=255)                                cATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAaaca                                                                                                                    >  1:199929/1‑149 (MQ=255)                                cATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAaaca                                                                                                                    >  1:642518/1‑149 (MQ=255)                                 aTCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAgacaa                                                                                                                   >  1:786071/1‑149 (MQ=255)                                 aTCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAaacaa                                                                                                                   >  2:485301/1‑149 (MQ=255)                                     aaCCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTg                                                                                                               >  2:739501/1‑149 (MQ=255)                                       ccGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGa                                                                                                             >  2:656836/1‑149 (MQ=255)                                        cGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGaa                                                                                                            <  1:110724/149‑1 (MQ=255)                                         gggTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTAGGCGGCATTAGTATAGAAGACAACGTGGATATCCACGAAAACAACGGGGaaa                                                                                                           >  2:137842/1‑149 (MQ=255)                                            tATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACa                                                                                                        <  1:166340/149‑1 (MQ=255)                                              tCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATg                                                                                                      <  1:270856/149‑1 (MQ=255)                                               cTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGa                                                                                                     <  1:573395/149‑1 (MQ=255)                                                tACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAAAGTGGTAATCCACGAAAAAAACGTGGAAAACATGaa                                                                                                    >  2:619205/1‑148 (MQ=255)                                                 aCTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGAcc                                                                                                   <  1:893214/149‑1 (MQ=255)                                                   ttCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATTGCAGACAACGGGGTGAACAAACAAAACAACGTGGAAAACATGACCCg                                                                                                 >  2:99196/1‑149 (MQ=255)                                                   ttCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACGACGCGGAAAACATGACCCg                                                                                                 >  2:177385/1‑149 (MQ=255)                                                   ttCATCGAATCGCGACTGGCACCGTGGCGTGAAGGGCAGTGCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCg                                                                                                 <  1:890846/149‑1 (MQ=255)                                                     cATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCATATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCggg                                                                                               >  2:923014/1‑149 (MQ=255)                                                     cATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCATATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCggg                                                                                               >  2:923013/1‑149 (MQ=255)                                                     cATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCATATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCggg                                                                                               >  2:922996/1‑149 (MQ=255)                                                     cATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTCTAGCAAGCACTTCGACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCggg                                                                                               <  2:660977/149‑1 (MQ=255)                                                       tCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGAt                                                                                             <  2:248499/149‑1 (MQ=255)                                                         gAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCt                                                                                           >  1:576272/1‑149 (MQ=255)                                                                 aCTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCGGAAACTgg                                                                                   >  2:91787/1‑149 (MQ=255)                                                                   tGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGCAAACAACGTGGAAAACATGACCCGGGATCTGAACCTggac                                                                                 >  2:328152/1‑147 (MQ=255)                                                                   tGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACAAGACCCGGGATCTGAAACTGGCg                                                                                 >  1:650133/1‑149 (MQ=255)                                                                     gCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTg                                                                               >  1:227044/1‑149 (MQ=255)                                                                     acACCGTGGCGTGAAGGGCCGGTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACGGGTCGGCGGCATTCGTATCGAAGCCAACGTGGTGATCCAGGGGAACAACGTGGAAAACATGACCCGGGACGTGACATTGGCggg                                                                               >  2:83755/2‑147 (MQ=255)                                                                       agcGTGTCGTGAAGGGTAGTTCATCAAGCTCTTCAACTGGCAGAAAATTGAAGCACTGAAACCTTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTCAt                                                                             <  2:715008/147‑1 (MQ=255)                                                                              cgggAAGGGCGGTTCAGCAAGCACGTCAACTGGCAGAAAATTGGAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGc                                                                      <  2:88021/146‑1 (MQ=255)                                                                               gTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCCTTCGTATCGAAGACAACGTGGTGATCCACGAACGCTACGTGGCAAACATGACCCGGGATCTGAAACTGGCATCATGGACAGCt                                                                     >  1:22785/1‑149 (MQ=255)                                                                                tGAAGGGCGGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGCCGTGTCTAACGAAAACAACGTGGAAAACATGACCCGGGATCTGATACTGGCGTGATGGACAGCTg                                                                    >  1:681564/1‑149 (MQ=255)                                                                                     ggCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTaa                                                               >  1:657106/1‑149 (MQ=255)                                                                                      gCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAAt                                                              >  1:79113/1‑149 (MQ=255)                                                                                          ttCAGCAAGCACTTCAACTGGCAGAAAATGGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCt                                                          >  1:623442/1‑149 (MQ=255)                                                                                              gcaagcaCTTCCACGGGCAGAAAATTGAAGGAATGAACCCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCgg                                                      <  1:443512/149‑1 (MQ=255)                                                                                                aaccaCTTCAATCTGCACAAAATTGAATCACTAAAACCATTCGTCGCCATTCGTATCGAAGACAACGTGGTAATCCACGAAAACAACTTGCAAAACATGACCCTCTATCTGAAACTGACTTAATTAAAACCTCCTTAATTCCTGCGGCa                                                    <  2:487801/146‑1 (MQ=255)                                                                                                 agcaCTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCAc                                                   <  2:720881/149‑1 (MQ=255)                                                                                                     cTTCAACTGGCAGAAAATTGAAGCCCTGAAACCGTTCGGCGGCATTCGTACCGCAGACAACGTGGTGAACCAAGAAGACAACGTGGAAAGCATGACCCGGGATTTGGAACAGGCGTGATGGGAAGCGGGTTGATGCCTGTGGCACAggg                                               >  2:73811/1‑148 (MQ=255)                                                                                                      ttCAACTGGCAGAAAATTGAAGCACTGAAACCGTACGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACACCGTGGAAAACATGACCCGGGACCTGAAACTGCCGTGAAGGAAAGCTCGTTAATTCCTGCAGCACCGGTc                                              >  1:397245/1‑149 (MQ=255)                                                                                                       tCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGGGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCAGGCTAATTCCTGCGGCACCGGTCa                                             >  1:107531/1‑149 (MQ=255)                                                                                                        cAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCAc                                            >  1:252015/1‑149 (MQ=255)                                                                                                            tGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACggac                                        >  1:131013/1‑147 (MQ=255)                                                                                                                aggaaaTTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTg                                    <  1:449624/146‑1 (MQ=255)                                                                                                                 gAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGa                                   >  1:768462/1‑149 (MQ=255)                                                                                                                  aaaaTTGAAGCACTGAAACCGTTCGGCGGCCTTCGCTCCAGATACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGaa                                  >  1:468349/1‑149 (MQ=255)                                                                                                                  aaaaTTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGGCCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGaa                                  <  2:533177/149‑1 (MQ=255)                                                                                                                  aaaaTTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGaa                                  <  2:681062/149‑1 (MQ=255)                                                                                                                  aaaaTTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATATGAATCTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGaa                                  <  2:442218/149‑1 (MQ=255)                                                                                                                  aaaaTTGAAGCACTGAAACAGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGaa                                  <  1:33705/149‑1 (MQ=255)                                                                                                                      ttGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAgaga                              <  2:576272/149‑1 (MQ=255)                                                                                                                      taaaaccACTGAAACCGTTCGGTGGCTTTCCGAACCAAGACAACGAAGTGAACAAAGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAgaga                              <  2:894769/143‑1 (MQ=255)                                                                                                                            cACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGGTCaaaa                        >  1:742000/1‑149 (MQ=255)                                                                                                                             aCTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCaaaaa                       <  2:484430/149‑1 (MQ=255)                                                                                                                              cTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGGCAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAag                      >  2:570627/1‑149 (MQ=255)                                                                                                                               tGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAaga                     <  2:306227/149‑1 (MQ=255)                                                                                                                               tGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCGGCGGCACCGGTCACGGTCGTTAAAGAGATCAAAAaga                     >  1:84810/1‑149 (MQ=255)                                                                                                                                      gTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCtt              >  2:3972/1‑149 (MQ=255)                                                                                                                                        tCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCa            <  1:818146/149‑1 (MQ=255)                                                                                                                                          ggcggcATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAACAGAGCCGCTTCATa          >  1:491588/1‑149 (MQ=255)                                                                                                                                            cggcATTAGAATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAAc        <  1:1586/149‑1 (MQ=255)                                                                                                                                             ggcATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGGGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACg       >  2:295623/1‑149 (MQ=255)                                                                                                                                              gcATCCATATGGAAGACAACGTGGTGATACCCGAAACCCCTGTGATAAACATGACCCGGGATCTGCATCTTTTCTGATGGAACGCAAAACTTTCTTTTCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGa      <  1:923013/149‑1 (MQ=255)                                                                                                                                              gaaTTCATATCGAAGACAACGTGGTGATCCACGAACACCCCGTGCGCGTCATGACCCGGGATCCGCCGCCAGGGTGACGGAGTAGTGTCGAGTGGCTGCGGCACCGGGCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGa      <  1:923014/147‑1 (MQ=255)                                                                                                                                                  tCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGAAGtt  >  2:341920/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAAAACATGACCCGGGATCTGAAACTGGCGTGATGGAAAGCTGGTTAATTCCTGCGGCACCGGTCACGGTCGTTGAAGAGATCAAAAAGAGCCGCTTCATAACGATGTT  >  NC_000913_3_pae_tpiA/4032275‑4032568

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