Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,316,893	C→T	E212K (GAA→AAA)	nusA ←	transcription termination/antitermination L factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,316,893	0	C	T	100.0%	340.9 / NA	94	E212K (GAA→AAA)	nusA	transcription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (46/48);  total (46/48)

CTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATAT  >  NC_000913_3_bme_pgi/3316747‑3317030                                                                                                                                                   |                                                                                                                                          cTGAACACGCGCGCCACGCATACCTGCGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAAt                                                                                                                                        >  4:11788/1‑150 (MQ=255)  tGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATa                                                                                                                                       <  1:300973/150‑1 (MQ=255)     acacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACgg                                                                                                                                    >  3:432533/1‑150 (MQ=255)       acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGaa                                                                                                                                  <  3:373702/150‑1 (MQ=255)           gcgcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGt                                                                                                                              >  2:1043486/1‑150 (MQ=255)              ccACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCg                                                                                                                           <  2:701166/150‑1 (MQ=255)                aCGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTCTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGAt                                                                                                                         <  1:723464/150‑1 (MQ=255)                 cGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATc                                                                                                                        <  1:444821/150‑1 (MQ=255)                 cGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGAGATTTTAGAACGAGAAAAAGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATc                                                                                                                        <  3:622842/150‑1 (MQ=255)                  gCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAGCCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCa                                                                                                                       >  1:657316/1‑150 (MQ=255)                        cTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                     >  1:1031560/1‑146 (MQ=255)                        cTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                     <  2:1031560/146‑1 (MQ=255)                         tACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTc                                                                                                                <  2:846969/150‑1 (MQ=255)                              aaGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCt                                                                                                           >  2:522095/1‑150 (MQ=255)                               aGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                          <  1:119518/150‑1 (MQ=255)                                   ccTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGaa                                                                                                      <  2:657316/150‑1 (MQ=255)                                    cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                                     >  4:440381/1‑150 (MQ=255)                                    cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                                     >  4:567071/1‑150 (MQ=255)                                     tACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACg                                                                                                    >  1:894362/1‑150 (MQ=255)                                      aCCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGa                                                                                                   <  1:1017300/150‑1 (MQ=255)                                       ccGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                  >  1:377890/1‑150 (MQ=255)                                       ccGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                  >  2:352010/1‑150 (MQ=255)                                        cGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGGTCAGCATTTCCGGCTTGGAACGAGt                                                                                                 >  4:653577/1‑150 (MQ=255)                                        cGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGt                                                                                                 <  2:2284/150‑1 (MQ=255)                                         ggATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTg                                                                                                >  1:338439/1‑150 (MQ=255)                                         ggATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTg                                                                                                >  3:467799/1‑150 (MQ=255)                                          gATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGa                                                                                               <  2:894362/150‑1 (MQ=255)                                           aTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGAc                                                                                              >  1:1026363/1‑150 (MQ=255)                                                   gTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATGAGCATTTCCGGCTTGGAACGAGTGACGAACAGtt                                                                                      >  1:780142/1‑150 (MQ=255)                                                       aTCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgcgc                                                                                  <  1:522095/150‑1 (MQ=255)                                                          gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                     >  4:217758/1‑112 (MQ=255)                                                          gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                     <  3:217758/112‑1 (MQ=255)                                                           ttGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCa                                                                              >  2:82471/1‑150 (MQ=255)                                                           ttGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCa                                                                              >  4:467103/1‑150 (MQ=255)                                                            tGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAc                                                                             >  2:65502/1‑150 (MQ=255)                                                            tGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAc                                                                             <  2:424547/150‑1 (MQ=255)                                                               ttttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCAttt                                                                                                                 >  3:291469/1‑111 (MQ=255)                                                               ttttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCAttt                                                                                                                 <  4:291469/111‑1 (MQ=255)                                                                 ttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCt                                                                        >  1:13077/1‑150 (MQ=255)                                                                    aCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTTACGAACAGTTGCGCGCCACGCGCTTcc                                                                     >  4:593635/1‑150 (MQ=255)                                                                    aCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTcc                                                                     >  1:316049/1‑150 (MQ=255)                                                                    aCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTcc                                                                     >  1:747551/1‑150 (MQ=255)                                                                      cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                   >  1:1035882/1‑150 (MQ=255)                                                                      cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                   <  2:338439/150‑1 (MQ=255)                                                                      cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                   <  2:377890/150‑1 (MQ=255)                                                                      cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                   >  4:164683/1‑150 (MQ=255)                                                                      cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                   >  4:178586/1‑150 (MQ=255)                                                                      cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                   <  1:82471/150‑1 (MQ=255)                                                                      cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                   >  2:762455/1‑150 (MQ=255)                                                                        cgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGc                                                                 <  1:689980/150‑1 (MQ=255)                                                                            tttCGCACGGGAAGGCGTATCGGTAGCGGCTCCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGATCAGTTCGATCAGCTTTTCCGGCTTGGAACGAGTGACCAACATTTGCGCGCCACTCTCTTCCGGGCGAAc                                                             <  2:534154/150‑1 (MQ=255)                                                                            tttCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAAc                                                             <  2:316049/150‑1 (MQ=255)                                                                              tCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACgg                                                           <  3:467103/150‑1 (MQ=255)                                                                                  aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTg                                                                                                >  1:722954/1‑109 (MQ=255)                                                                                  aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTg                                                                                                <  2:722954/109‑1 (MQ=255)                                                                                  aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                       <  4:584635/150‑1 (MQ=255)                                                                                  aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                       <  3:185890/150‑1 (MQ=255)                                                                                    gagaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATaga                                                     <  3:178586/150‑1 (MQ=255)                                                                                      gaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                  >  2:646/1‑135 (MQ=255)                                                                                      gaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                  <  1:646/135‑1 (MQ=255)                                                                                         cccGGATCGCGAGCCGCTGCTTTAATTTCGATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAAc                                                                                          >  2:63505/1‑108 (MQ=255)                                                                                         cccGGATCGCGAGCCGCTGCTTTAATTTCGATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAAc                                                                                          <  1:63505/108‑1 (MQ=255)                                                                                         cccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACg                                                >  4:569396/1‑150 (MQ=255)                                                                                            ggATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGaa                                                                                                                                  <  2:86075/65‑1 (MQ=255)                                                                                            ggATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGaa                                                                                                                                  >  1:86075/1‑65 (MQ=255)                                                                                            ggATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCa                                             <  1:634705/150‑1 (MQ=255)                                                                                             gATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCAc                                            <  4:500913/150‑1 (MQ=255)                                                                                                  cgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTg                                                                                                         <  2:429981/84‑1 (MQ=255)                                                                                                  cgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTg                                                                                                         >  1:429981/1‑84 (MQ=255)                                                                                                  cgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcg                                       <  2:747551/150‑1 (MQ=255)                                                                                                   gAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcgc                                      >  1:825154/1‑150 (MQ=255)                                                                                                   gAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcgc                                      <  2:437892/150‑1 (MQ=255)                                                                                                    aGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcgcg                                     <  3:11788/150‑1 (MQ=255)                                                                                                    aGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTGCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcgcg                                     <  1:551015/150‑1 (MQ=255)                                                                                                      ccGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGt                                   >  1:218407/1‑150 (MQ=255)                                                                                                           gctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCa                              <  1:352010/150‑1 (MQ=255)                                                                                                            ctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAg                             >  1:606002/1‑150 (MQ=255)                                                                                                             tttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAgg                            <  1:762455/150‑1 (MQ=255)                                                                                                              ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAggg                           >  1:429891/1‑150 (MQ=255)                                                                                                                  tttCAATCACTTCTTCGCCGGTTTCTGGCACTTTAATACGGAACAGtt                                                                                                                             <  1:769561/48‑1 (MQ=255)                                                                                                                  tttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGtt                                                                                                                             >  2:769561/1‑48 (MQ=255)                                                                                                                  tttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGa                       >  1:1066354/1‑150 (MQ=255)                                                                                                                   ttCAATCCCTTCTTCGCCGATTTCTGGAACTTTAATACGGAACAGTCCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGaa                      <  4:389489/150‑1 (MQ=255)                                                                                                                   ttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCa                              <  1:89930/142‑1 (MQ=255)                                                                                                                   ttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCa                              >  2:89930/1‑142 (MQ=255)                                                                                                                       aTCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGttt                  <  2:80575/150‑1 (MQ=255)                                                                                                                            ttcttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcgcg                                     <  4:307315/126‑1 (MQ=255)                                                                                                                            ttcttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcgcg                                     >  3:307315/1‑126 (MQ=255)                                                                                                                            ttcttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAcg             <  3:567071/150‑1 (MQ=255)                                                                                                                            ttcttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAcg             <  3:440381/150‑1 (MQ=255)                                                                                                                            ttcttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAcg             >  1:312967/1‑150 (MQ=255)                                                                                                                               ttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCgg          >  2:460388/1‑150 (MQ=255)                                                                                                                                gggCCGATGGCTGGCACTTTAATACGTAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCGCGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgc         <  4:520125/148‑1 (MQ=255)                                                                                                                                       tttCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCatat  >  3:527219/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                          CTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATAT  >  NC_000913_3_bme_pgi/3316747‑3317030

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