Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_pae_tpiA	4,186,457	A→G	45.3%	K398E (AAG→GAG)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	4,186,457	0	A	G	45.3%	16.3 / 63.3	55	K398E (AAG→GAG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (11/19);  new base G (13/12);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTA  >  NC_000913_3_pae_tpiA/4186313‑4186604                                                                                                                                                  |                                                                                                                                                    tCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCgaagaa                                                                                                                                                 >  1:62917/1‑149 (MQ=255)   cGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGgaaa                                                                                                                                               >  1:652635/1‑149 (MQ=255)   cGGTCGTTCCTTAATCACCGTACGTCCCTACCTGCGTCCGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCCTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGaaaa                                                                                                                                               <  2:222399/149‑1 (MQ=255)     gTCGTTCGGTAATCACCGCAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATg                                                                                                                                             <  1:358635/149‑1 (MQ=255)      tCGCTCTGTAATCACCTTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATgg                                                                                                                                            <  2:197234/149‑1 (MQ=255)       cGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGt                                                                                                                                           <  2:20838/149‑1 (MQ=255)        gTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCAGTCGGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAAATGCGTGGTCTTGCTACGCACC‑TTAAAGCTGCGACGaaaaaagtt                                                                                                                                          >  2:855621/1‑144 (MQ=255)            tGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGacc                                                                                                                                      >  2:203022/1‑147 (MQ=255)              tAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAAAAAATGGTTGAgcgc                                                                                                                                    >  1:383916/1‑149 (MQ=255)               caTCACCGTATGTCCATACCTGTGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAAACGTCCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTTAGcgcg                                                                                                                                   <  2:467933/148‑1 (MQ=255)                 tCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCgaa                                                                                                                                 >  1:301769/1‑149 (MQ=255)                    ccGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCgaagaa                                                                                                                              <  1:218679/149‑1 (MQ=255)                     cGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAACCTGCGGAGAAAATGATTGAGCGCaaagagg                                                                                                                             >  1:796881/1‑147 (MQ=255)                         ggTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGt                                                                                                                         <  1:179720/149‑1 (MQ=255)                          gTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCCCCCCGGTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTc                                                                                                                        <  1:222739/149‑1 (MQ=255)                          gTCAATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGACAATGGTTGAGCGCAAAGAAGCTGTc                                                                                                                        >  1:594244/1‑149 (MQ=255)                               tACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTg                                                                                                                   <  2:578078/149‑1 (MQ=255)                                  cTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGCGCTCTTCGGCC‑CACCAGTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGa                                                                                                                >  1:221149/1‑149 (MQ=255)                                  cTGCGGCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCCCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTATAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGa                                                                                                                <  2:173894/149‑1 (MQ=255)                                        cTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCt                                                                                                          <  2:855545/149‑1 (MQ=255)                                         tGCATCACTGCGGGCGGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCAAC‑AACAATAAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTg                                                                                                         <  1:184408/149‑1 (MQ=255)                                                           cGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGaa                                                                                       >  1:850513/1‑149 (MQ=255)                                                             aagaaAATGGCACTTGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAaca                                                                                     >  1:561866/1‑149 (MQ=255)                                                               gaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACAcc                                                                                   <  1:870909/149‑1 (MQ=255)                                                                  aaTGGCACTGGGGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATGCGCGAACACCCgg                                                                                <  1:631860/149‑1 (MQ=255)                                                                  aaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCgg                                                                                <  2:675495/149‑1 (MQ=255)                                                                  aaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGATCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCgg                                                                                >  2:538708/1‑149 (MQ=255)                                                                    tGGCCCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTa                                                                              <  2:537710/149‑1 (MQ=255)                                                                        aCTGTAGCTGTTCAAACCGTTAATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGAActgc                                                                          >  2:349344/1‑149 (MQ=255)                                                                                  ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGc                                                                >  2:660556/1‑149 (MQ=255)                                                                                    cAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                              <  1:660556/149‑1 (MQ=255)                                                                                     aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                             >  1:845917/1‑149 (MQ=255)                                                                                     aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                             >  1:211944/1‑149 (MQ=255)                                                                                     aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACTAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                             <  2:628546/149‑1 (MQ=255)                                                                                      aaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCg                                                            >  1:841387/1‑149 (MQ=255)                                                                                              acgtACGGCAAGTGGGAACTGCGTGATACTTCTAC‑CATCATTAACGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGCGATATCCTGCACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCa                                                    <  2:847210/146‑1 (MQ=255)                                                                                                 tACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCg                                                 >  1:656980/1‑149 (MQ=255)                                                                                                    ggCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                              >  1:707648/1‑149 (MQ=255)                                                                                                      cagtcTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                            <  2:44792/145‑1 (MQ=255)                                                                                                      cAAGTTGGAACTGCGTGGTCTTGCTAC‑CACGTTTAAAGCTTCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                            <  1:871321/149‑1 (MQ=255)                                                                                                      cAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                            <  2:924149/149‑1 (MQ=255)                                                                                                      cAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                            >  1:292631/1‑149 (MQ=255)                                                                                                       aaGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                           <  2:794977/149‑1 (MQ=255)                                                                                                           tGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATc                                       >  2:604231/1‑149 (MQ=255)                                                                                                                 tGCGTGGTCTTGCTAC‑TACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGTCTCTGCACCGTCTGGGTATCCAGGCa                                 <  1:535336/149‑1 (MQ=255)                                                                                                                     tGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGACGCTGTCGTTTGGGATAGCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTg                             <  1:774399/149‑1 (MQ=255)                                                                                                                          ttGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTGCATTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGCCTGGGTATCCAGGCATTAGAACCg                        >  1:645683/1‑149 (MQ=255)                                                                                                                          ttGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                            >  2:59129/1‑129 (MQ=255)                                                                                                                          ttGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                            <  1:59129/129‑1 (MQ=255)                                                                                                                          ttGCTAC‑CACCATTAAAGCTGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                        >  2:734588/1‑149 (MQ=255)                                                                                                                             cTAC‑CACCATTAAAGCTGAGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGTAGTTATCCGCGAACACGCGGTACTGCTGAGCCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTa                     <  2:901093/149‑1 (MQ=255)                                                                                                                                        tAAAGTTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGCCGTGTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg           <  2:128056/149‑1 (MQ=255)                                                                                                                                             cTGCGGAGTAAATGGCTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGTACGAATTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCCCCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa      <  2:106283/149‑1 (MQ=255)                                                                                                                                             cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaaa      <  1:862836/149‑1 (MQ=255)                                                                                                                                              tGCGGAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg     <  2:483310/149‑1 (MQ=255)                                                                                                                                               gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTAGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCCTTTGAACCGGTACTGCTCGAAGGTATAGc    >  2:534851/1‑149 (MQ=255)                                                                                                                                                 gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACGGATCGAAGGTAAAGCTa  >  2:447118/1‑149 (MQ=255)                                                                                                                                                 gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTGATTTTCTAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTa  >  1:271989/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAC‑CACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTA  >  NC_000913_3_pae_tpiA/4186313‑4186604

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