Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,725,672	(C)5→4	coding (75/1299 nt)	kgtP ←	alpha‑ketoglutarate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,725,668	0	C	.	100.0%	272.7 / NA	57	coding (79/1299 nt)	kgtP	alpha‑ketoglutarate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (31/26);  total (31/26)

GCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTAA  >  NC_000913/2725521‑2725810                                                                                                                                                    |                                                                                                                                               gcATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGcccc                                                                                                                                             <  1:980371/151‑1 (MQ=255) gcATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGcccc                                                                                                                                             >  1:660674/1‑151 (MQ=255) gaaTCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGcc                                                                                                                                               >  1:112175/3‑149 (MQ=255)   aTCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGccccac                                                                                                                                           >  2:890378/1‑149 (MQ=255)     cAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGccccaca                                                                                                                                          >  2:834182/1‑147 (MQ=255)      aGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAAt                                                                                                                                       >  2:275168/1‑151 (MQ=255)       ggAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATc                                                                                                                                      <  1:763797/151‑1 (MQ=255)             ccGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCa                                                                                                                                 >  1:426383/1‑150 (MQ=255)               gcagccAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAAt                                                                                                                              >  1:875936/1‑151 (MQ=255)                      aaaCACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGcga                                                                                                                       <  2:1045844/151‑1 (MQ=255)                       aaCACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGcgac                                                                                                                      <  2:565122/151‑1 (MQ=255)                        acacCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGcgacg                                                                                                                     <  2:173207/151‑1 (MQ=255)                         cacCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGcgacg                                                                                                                     >  1:164541/1‑150 (MQ=255)                             tgctgTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTAt                                                                                                                >  1:1063303/1‑151 (MQ=255)                                tgTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCac                                                                                                             >  2:1057165/1‑151 (MQ=255)                                    tGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCactact                                                                                                         >  2:625949/1‑151 (MQ=255)                                        gtagtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTc                                                                                                     >  2:415542/1‑151 (MQ=255)                                        gtagtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTc                                                                                                     >  1:780240/1‑151 (MQ=255)                                             tGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGtt                                                                                                 <  1:674361/150‑1 (MQ=255)                                                  cgtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCTACTAGTATCACTACTTGTCAGTTTGCTg                                                                                            <  1:550650/150‑1 (MQ=255)                                                   gtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTc                                                                                          >  1:951157/1‑151 (MQ=255)                                                     cgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTg                                                                                        >  2:305384/1‑151 (MQ=255)                                                      gtgtTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGc                                                                                       <  2:977419/151‑1 (MQ=255)                                                      gtgtTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTATCAGTTTGCTGTCTGc                                                                                       >  1:552594/1‑151 (MQ=255)                                                            ccGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTAc                                                                                 >  2:851320/1‑151 (MQ=255)                                                            ccGGAAGGAACGAAGATGTGAGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAAACGAACCACAAGACCAGATTACCTGCAGAAG‑AACAACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTAAGTTTGCTGTCTGCCGTTAc                                                                                 <  2:646655/151‑1 (MQ=255)                                                              ggaagggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTAc                                                                                 >  1:198154/1‑149 (MQ=255)                                                               gaagggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGt                                                                              >  1:192851/1‑151 (MQ=255)                                                                 agggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTAc                                                                            >  1:497880/1‑151 (MQ=255)                                                                  gggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACt                                                                           <  1:694563/151‑1 (MQ=255)                                                                     aagaagATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTc                                                                        >  2:714010/1‑151 (MQ=255)                                                                        aagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgc                                                                     >  2:463/1‑151 (MQ=255)                                                                         agaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgcc                                                                    >  2:1011541/1‑151 (MQ=255)                                                                          gaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgcca                                                                   <  1:291183/151‑1 (MQ=255)                                                                            tgtgGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccatg                                                                 >  2:733708/1‑151 (MQ=255)                                                                             gtgGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccatgc                                                                >  1:690580/1‑151 (MQ=255)                                                                             gtgGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccatg                                                                 <  2:1063480/150‑1 (MQ=255)                                                                               gggCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccatgcca                                                              <  2:1039844/151‑1 (MQ=255)                                                                                  cAAGGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCAtt                                                            >  2:603074/1‑150 (MQ=255)                                                                                        agagTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGtctc                                                     >  2:25636/1‑151 (MQ=255)                                                                                          agTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCt                                                   >  1:569486/1‑151 (MQ=255)                                                                                          agTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCt                                                   <  1:508584/151‑1 (MQ=255)                                                                                           gTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTg                                                  <  2:746115/151‑1 (MQ=255)                                                                                                 aGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGCACTTTCAGCCATGCCATTATGTCTCCTGCCGTa                                             <  2:297894/150‑1 (MQ=255)                                                                                                        gTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATg                                     <  2:913208/151‑1 (MQ=255)                                                                                                          aGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCt                                   <  1:275100/151‑1 (MQ=255)                                                                                                             cATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCtttt                                >  2:949218/1‑151 (MQ=255)                                                                                                             cATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCtttt                                >  1:672007/1‑151 (MQ=255)                                                                                                                aGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgt                              <  2:692475/150‑1 (MQ=255)                                                                                                                        aaCCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCtt                     <  2:521734/151‑1 (MQ=255)                                                                                                                        aaCCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCtt                     >  2:619676/1‑151 (MQ=255)                                                                                                                          ccACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCttt                    <  1:280375/150‑1 (MQ=255)                                                                                                                           cccTCGACCAGATTCCCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTCCTGTCTGCGGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTg                  <  1:890229/149‑1 (MQ=255)                                                                                                                           cACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTg                  >  2:730543/1‑151 (MQ=255)                                                                                                                                gACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCTGTCGCTTTTGTTTAt             >  2:455540/1‑151 (MQ=255)                                                                                                                                    aGATTCCCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAttttt         <  1:947943/151‑1 (MQ=255)                                                                                                                                         aCCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTg     <  2:376682/150‑1 (MQ=255)                                                                                                                                          ccTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTa   <  2:646305/151‑1 (MQ=255)                                                                                                                                          ccTAAAGAGG‑CCCCACATTCGCCCAAATGCGGCGACGAGTATCACTCCTTGTCAGTTTGCTGTCTTCCGTTACAGTCCTTTCAGCCTTGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAtttttttttt    <  1:1064565/150‑3 (MQ=255)                                                                                                                                           cTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTa   >  1:1004709/1‑150 (MQ=255)                                                                                                                                            tGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTaa  >  2:68265/1‑150 (MQ=255)                                                                                                                                             gAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTaa  <  2:102142/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               GCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTAA  >  NC_000913/2725521‑2725810

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