Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	1,265,009	C→A	intergenic (+39/‑3)	hemA → / → prfA	glutamyl tRNA reductase/peptide chain release factor RF‑1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	1,265,009	0	C	A	96.9%	162.9 / ‑3.9	66	intergenic (+39/‑3)	hemA/prfA	glutamyl tRNA reductase/peptide chain release factor RF‑1
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (31/32);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

TGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGAT  >  NC_000913_3_pae_tpiA/1264864‑1265152                                                                                                                                                   |                                                                                                                                                tGAGTAACCGCTAGACCCATGCGCAAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTa                                                                                                                                              <  2:681233/149‑1 (MQ=255) tGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGGGACGGGGATAACGAATGCCTGAATATGCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTa                                                                                                                                              <  2:260976/149‑1 (MQ=255)  gACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTAt                                                                                                                                             >  1:1040230/1‑149 (MQ=255)    cTAACCGCTTGATACATGCGCCAAAGAAATCACGACAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑CCGACTATGa                                                                                                                                           >  2:1010583/1‑149 (MQ=255)     tAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACGGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGTGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGaa                                                                                                                                          >  1:619516/1‑149 (MQ=255)        ccGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGGGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCCGTACATCATTTTCTTTTTTTACAGGGTGCATCT‑TCGCCTATGATGcc                                                                                                                                       >  2:336612/1‑149 (MQ=255)           cTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTc                                                                                                                                    >  1:1013324/1‑149 (MQ=255)             tGATCCATGCGGCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTTAATTTTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTa                                                                                                                                  <  1:28339/149‑1 (MQ=255)              gATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGCCTATGAAGCCTTCTAt                                                                                                                                 <  1:980703/149‑1 (MQ=255)                  cATGCGCCTACGAAATCACTTCAACAGGCTGCCCGTGACGGAGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCACTT‑ACGACTATGAAGCCTTCTATCGtt                                                                                                                             <  1:608707/149‑1 (MQ=255)                    tGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGc                                                                                                                           <  1:245817/149‑1 (MQ=255)                        ccAACGAAATCCGTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCaaa                                                                                                                       <  2:1025821/149‑1 (MQ=255)                         cAACGAAATCACTTTAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAAc                                                                                                                      <  1:394761/149‑1 (MQ=255)                         cAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTATATGGTGTATTT‑ACGACTATGAGACCTTCTATCGTTGCCAAAc                                                                                                                      >  2:837403/1‑149 (MQ=255)                          aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAATGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACTACTATGAAGCCTTCTATCGTTGCCAAACt                                                                                                                     >  1:221087/1‑149 (MQ=255)                          aaCGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGGGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTGTGAAGCCTTCTATCGTTGACAACCt                                                                                                                     >  1:36526/1‑149 (MQ=255)                          aaCGAAACCCCTTCAACAGGCTGCCCGTGACGGGGATAAGGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGTGTGCATTT‑ACGACTATGAAGCTTTCTATCGTTGCCAAACt                                                                                                                     <  2:444746/149‑1 (MQ=255)                              aaaTCACGTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTGCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGaa                                                                                                                 <  1:512357/149‑1 (MQ=255)                               aaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAg                                                                                                                <  1:468763/149‑1 (MQ=255)                               aaTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTAAAGCATTTTCTTTTTTTACAGGCTGCATTT‑ACGACTTTGCAGCCCTCTATTGTTACCAAACTGGaat                                                                                                                >  2:313220/1‑148 (MQ=255)                                aTCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGc                                                                                                               >  2:750473/1‑149 (MQ=255)                                  cACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTATATAGTTACCAAACTGAAAGcca                                                                                                             >  1:223433/1‑148 (MQ=255)                                   aCTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCt                                                                                                            >  1:374806/1‑149 (MQ=255)                                     ttCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGCCAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGAGTCCATTT‑ACGACTATAAAGCCTTCTAACGTTGCCAAACTGGAAGCCCTGc                                                                                                          >  1:162770/1‑149 (MQ=255)                                        aaCAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACACCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGAATTT‑ACGACTATAAATCCTTCTATCGGTGCTAAACACGAAGCCCTGCATg                                                                                                       >  2:174934/1‑149 (MQ=255)                                           aGGCCGCCTGTGACGGGGAGAACGAACGGCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTCTTTTACAGGATGCAGTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAAc                                                                                                    <  2:955963/149‑1 (MQ=255)                                           aGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAAc                                                                                                    <  1:673552/149‑1 (MQ=255)                                           aGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTCGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAAc                                                                                                    <  2:910588/149‑1 (MQ=255)                                             gccgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCCCAGGCTGGAGTAGCAGTAGACCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAACACTTCTATCGTTGCCAAACTGGAAGCCCCACATGCACGc                                                                                                  >  2:119079/1‑149 (MQ=255)                                               cgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCa                                                                                                <  2:817627/149‑1 (MQ=255)                                               cgccCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCa                                                                                                <  1:914681/149‑1 (MQ=255)                                                   cGTGACGGGGATAACGAACGCCTGAATATTTTGCGCGACAGGCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATgaa                                                                                            <  2:811633/149‑1 (MQ=255)                                                   cGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACATGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATgaa                                                                                            >  1:529054/1‑149 (MQ=255)                                                        cGGGGATAACGAACGACTGAATATTCTCCGCGACAGCCTCGCGCTGGAGTAGACGTACAGCATTGTCGCGTGTTACAGGGTGCATGG‑ACGACTATGAAGCGTTCGATCTTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGt                                                                                       <  2:1030050/149‑1 (MQ=255)                                                         ggggATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTTCACG‑CTTTGAAACCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGtt                                                                                      >  1:501558/1‑149 (MQ=255)                                                           ggATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAATCTGGACGCCCTGCATGAACGCCTTGAAGAAGTTCa                                                                                    >  2:70863/1‑149 (MQ=255)                                                               aacgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGACTTCTATCGTTGCCAAACTGGAAGCCCTGCATAAACGCAATGAAGACATTCAGGCg                                                                                >  1:394541/1‑149 (MQ=255)                                                                acgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGtt                                                                                                                             <  2:230411/103‑1 (MQ=255)                                                                acgaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGtt                                                                                                                             >  1:230411/1‑103 (MQ=255)                                                                  gaacgCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACACCATGAAGAATTTCAGGCGTTg                                                                             >  2:186678/1‑149 (MQ=255)                                                                  gaacgCCTGAATATTCTGAGCGAAAGCCTCGGGGTGGAGTAGCAGAAAATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCGAAACTGGAAGACCTGCATGAACGACATGAAGAAGTTAAGGCATTg                                                                             >  2:795222/1‑149 (MQ=255)                                                                  gaacgCCTGAATATTCTACGCGACAGCCTCGGGCTGGAGTAGCAGTCCATCATTTTCTTAGTCTAACGGGTGCATTT‑ACTACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCCGCACGAACGCCATAAAGACGTTCAGGctgcg                                                                             >  2:116034/1‑145 (MQ=255)                                                                      ggcTGAATATGCGGCGCGACGGCGCCGGGCTGGAGTAGCAGTACATCATTTTCTTGGTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTgg                                                                         <  1:921774/147‑1 (MQ=255)                                                                        cTGAATATCCTGCGCGACAGCCTCGGGCTGGAGCAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGt                                                                       <  2:925448/149‑1 (MQ=255)                                                                           aaTATTCTGCGCGACAGCCTGGGGCTGGAGTAGCAGTACATCATTTTTTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATAGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGAc                                                                    <  2:171663/149‑1 (MQ=255)                                                                           aaTATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGAc                                                                    >  1:861690/1‑149 (MQ=255)                                                                                  cgctcgACAGCCTCCGGCTGGTGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCaaa                                                             <  1:520929/145‑1 (MQ=255)                                                                                       gACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATc                                                        <  2:529054/149‑1 (MQ=255)                                                                                            ccTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCATTCGATCGTGGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACCACCCCCGa                                                   >  2:995043/1‑149 (MQ=255)                                                                                             cTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTTTTTTTTTTATTT‑ACGACTATGAAGCCTTCTATCTTTTCCATACTTGAAGCCCTGCATTAACGCCTTTAAGAAGTTTCGGCGTCGCTGGGTGACGCGCCAACTTTCTCCgtt                                                  >  2:508597/1‑147 (MQ=255)                                                                                               cGGGCTGTAGTAGCAGTACAGCATTTTCTCTGTTTCCAGGGTACATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCa                                                <  2:283412/149‑1 (MQ=255)                                                                                                      gAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTAAGGCGTTGCTGGGCGACGCGCAAACTATCGCCGACCAGGAACGt                                         >  2:847480/1‑149 (MQ=255)                                                                                                       aGTAGCAGTATATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGtt                                        <  1:705283/149‑1 (MQ=255)                                                                                                          aGCAGTACATTATTTTCTTTTTTTACGGGGTGCATTT‑ACGACTATGAATCCTTCTATCGTTGCCAAACTGGAAGCCCTGAATGAACGACAGGAATATGTTTACGGCTTTCTTGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTc                                     <  1:770617/149‑1 (MQ=255)                                                                                                              gTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGCAGCCCATGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCgcgc                                 <  1:987901/149‑1 (MQ=255)                                                                                                              gTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCaaa                                                             <  2:8856/121‑1 (MQ=255)                                                                                                              gTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCaaa                                                             >  1:8856/1‑121 (MQ=255)                                                                                                               tACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGa                                                                     <  1:588804/112‑1 (MQ=255)                                                                                                               tACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGa                                                                     >  2:588804/1‑112 (MQ=255)                                                                                                                  atcatTTTCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTAGAAGCCCTGAATGAACGCCATGAAGGAGTTAAGGCGTTGCTGGGTGACGCGCAAAATATCGCCGCCCAGGAACGTTTTTGCGCATTa                             >  2:238930/1‑149 (MQ=255)                                                                                                                        ttCTTTTTTTACAGGGTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCAACCAGGAACGTTTTCGCGCATTTTCAcgc                       >  1:76364/1‑149 (MQ=255)                                                                                                                                gtccAGGGTTCATTT‑ACGAATATGAAGCCTTCTATCGTTGCCAAACGGGAAGCCCTGCATGAACGCCATGAAGAATCTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATgc               <  2:374806/146‑1 (MQ=255)                                                                                                                                       gTGCATTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAAAAGGAAGCCCTACAAGCAAACCATGAATGTGATCATGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCACGAACGTTTTCGCGCATTATCACGCGAATATGCCCAGTTa        >  2:119606/1‑149 (MQ=255)                                                                                                                                        cgCATTC‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGCAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCAACCAGGAACGTTTTCGCGCACTATCACGCGAATATGCGCAGTTaa       <  2:566292/148‑1 (MQ=255)                                                                                                                                           aTTT‑ACGACTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGTGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTg    >  2:927655/1‑149 (MQ=255)                                                                                                                                             tt‑ACGACTATGGAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGAt  <  1:508597/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                TGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGCTGGAGTAGCAGTACATCATTTTCTTTTTTTACAGGGTGCATTT‑ACGCCTATGAAGCCTTCTATCGTTGCCAAACTGGAAGCCCTGCATGAACGCCATGAAGAAGTTCAGGCGTTGCTGGGTGACGCGCAAACTATCGCCGACCAGGAACGTTTTCGCGCATTATCACGCGAATATGCGCAGTTAAGTGAT  >  NC_000913_3_pae_tpiA/1264864‑1265152

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