Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,221,366	T→C	K208K (AAA→AAG)	bglX ←	beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,221,366	0	T	C	100.0%	311.5 / NA	79	K208K (AAA→AAG)	bglX	beta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (42/37);  total (42/37)

GCGAGTTCAGCGCCACCATCACCGCGCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGTTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCGGATCGCGCGAGA  >  NC_000913/2221218‑2221508                                                                                                                                                     |                                                                                                                                               gCGAGTTCAGCGCCACCATCACCGCGCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGctt                                                                                                                                              <  2:1045627/151‑1 (MQ=255)            gcCACCATCACCGCGCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGtc                                                                                                                                    >  1:713974/1‑150 (MQ=255)             ccaccaTCACCGCGCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGtcat                                                                                                                                  <  2:87895/151‑1 (MQ=255)              caccaTCACCGCGCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGtcatc                                                                                                                                 <  1:567285/151‑1 (MQ=255)                    cACCGCGCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCg                                                                                                                             >  2:446969/1‑149 (MQ=255)                      ccGCGCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTa                                                                                                                         <  2:330330/151‑1 (MQ=255)                       cgcgCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGt                                                                                                                          >  1:359760/1‑149 (MQ=255)                       cgcgCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAg                                                                                                                        <  1:591275/151‑1 (MQ=255)                        gcgcCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGc                                                                                                                       >  1:184078/1‑151 (MQ=255)                          gccgcTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCg                                                                                                                      >  2:660485/1‑150 (MQ=255)                            cgcTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATc                                                                                                                   >  1:892137/1‑151 (MQ=255)                             gcTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCt                                                                                                                  <  2:686657/151‑1 (MQ=255)                                gcctgcGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGcc                                                                                                               >  1:417812/1‑151 (MQ=255)                                gcctgcGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGcc                                                                                                               <  1:40027/151‑1 (MQ=255)                                  ctgcGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCgg                                                                                                             >  2:1024557/1‑151 (MQ=255)                                   tgcGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCggg                                                                                                            >  1:578928/1‑151 (MQ=255)                                    gcGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCggg                                                                                                            <  2:548725/150‑1 (MQ=255)                                       tCCAGCCACGCTTTGTACGTCGGCACATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGCGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGAttt                                                                                                        >  2:111148/1‑151 (MQ=255)                                        ccAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCtttt                                                                                                       <  2:471577/151‑1 (MQ=255)                                           gCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTAcc                                                                                                    <  2:840164/151‑1 (MQ=255)                                           gCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTAcc                                                                                                    >  1:145886/1‑151 (MQ=255)                                             cccGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCt                                                                                                  >  2:183039/1‑151 (MQ=255)                                              ccGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTg                                                                                                 <  1:1038221/151‑1 (MQ=255)                                               cGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGc                                                                                                >  1:597136/1‑151 (MQ=255)                                                gCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGc                                                                                                >  1:967878/1‑150 (MQ=255)                                                  tttGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATc                                                                                             >  2:323467/1‑151 (MQ=255)                                                   ttGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATc                                                                                             <  1:98302/150‑1 (MQ=255)                                                    tGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCg                                                                                            <  1:326574/150‑1 (MQ=255)                                                    tGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGc                                                                                           >  2:666669/1‑151 (MQ=255)                                                      tACGGCGGCTTATAATCATTAAACAGGCGCGGCGGACTCATATCGACGGGGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGGGCTTGACGCCGGTCATCACCGAGTACCGATCGGCCGGGCTTTTACCCGCCATCCCt                                                                                          >  2:499067/1‑150 (MQ=255)                                                        cggcggCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTcc                                                                                       <  1:769216/151‑1 (MQ=255)                                                           cggcATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTccacc                                                                                    >  1:1035825/1‑151 (MQ=255)                                                               atatAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATgg                                                                                <  1:782474/151‑1 (MQ=255)                                                                 ataATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGtt                                                                              <  1:200794/151‑1 (MQ=255)                                                                   aaTCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGttt                                                                             <  1:858829/150‑1 (MQ=255)                                                                       aTTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTAccc                                                                        <  2:971615/151‑1 (MQ=255)                                                                         tAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCa                                                                       <  2:1038894/150‑1 (MQ=255)                                                                          aaaCAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCAtt                                                                     <  1:217170/151‑1 (MQ=255)                                                                            aCAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCAttgt                                                                   <  1:293951/151‑1 (MQ=255)                                                                              aGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCAttgttg                                                                 >  2:688093/1‑151 (MQ=255)                                                                              aGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCAttgttg                                                                 >  2:839377/1‑151 (MQ=255)                                                                                gcgcTGCGGACTCATATCGACGGTGATGTACTCGATACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCAttgttg                                                                 >  1:623488/1‑149 (MQ=255)                                                                                 cgcTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAgg                                                              >  2:7594/1‑151 (MQ=255)                                                                                   cTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTg                                                            >  2:1040512/1‑151 (MQ=255)                                                                                     gCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTg                                                            >  1:120026/1‑149 (MQ=255)                                                                                       ggACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTg                                                            >  1:160471/1‑147 (MQ=255)                                                                                         aCTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATAc                                                      >  1:820322/1‑151 (MQ=255)                                                                                         aCTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATAc                                                      >  2:357598/1‑151 (MQ=255)                                                                                         aCTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATAc                                                      >  2:1035224/1‑151 (MQ=255)                                                                                         aCTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATAc                                                      <  2:196338/151‑1 (MQ=255)                                                                                            cATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACg                                                     <  2:393787/149‑1 (MQ=255)                                                                                            cATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTa                                                   <  2:759162/151‑1 (MQ=255)                                                                                            cATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGG‑TTTACCCATTGTTGAGGTGAGATACGt                                                    <  1:1061490/149‑1 (MQ=255)                                                                                              tatCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATc                                                 >  1:682133/1‑151 (MQ=255)                                                                                                 cGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTc                                              <  2:593/151‑1 (MQ=255)                                                                                                 cGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTc                                              >  1:875979/1‑151 (MQ=255)                                                                                                   aCGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGc                                            >  2:845132/1‑151 (MQ=255)                                                                                                     ggTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGcc                                           >  2:640683/1‑150 (MQ=255)                                                                                                        gttgtACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCaaa                                        <  2:1015231/150‑1 (MQ=255)                                                                                                           gtACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCt                                    <  1:188795/151‑1 (MQ=255)                                                                                                           gtACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCt                                    <  1:208448/151‑1 (MQ=255)                                                                                                             aCTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCtt                                   <  2:970653/150‑1 (MQ=255)                                                                                                             aCTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTc                                  <  1:766280/151‑1 (MQ=255)                                                                                                              ctctTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCg                                 >  2:326435/1‑151 (MQ=255)                                                                                                                 tttACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGa                               >  1:541088/1‑150 (MQ=255)                                                                                                                  ttACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCgg                                <  1:795474/148‑1 (MQ=255)                                                                                                                  ttACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTCACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCgg                                <  1:370070/148‑1 (MQ=255)                                                                                                                           tCTACCGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTcc                      <  1:247154/149‑1 (MQ=255)                                                                                                                               ccGCGCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAgcg                <  1:314240/151‑1 (MQ=255)                                                                                                                                cgcgCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAgcgc               >  2:980659/1‑151 (MQ=255)                                                                                                                                cgcgCCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAgcgc               <  2:840241/151‑1 (MQ=255)                                                                                                                                 gcgcCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGcgcg              >  1:354701/1‑151 (MQ=255)                                                                                                                                 gcgcCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGcgcg              >  1:777119/1‑151 (MQ=255)                                                                                                                                  cgcCGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCgg             >  2:1018268/1‑151 (MQ=255)                                                                                                                                    ccGTATGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCGGAt           <  1:726255/151‑1 (MQ=255)                                                                                                                                        aTGCGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCGGATcgcg       >  1:312193/1‑151 (MQ=255)                                                                                                                                           cGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCGGATcgcgcg     >  1:211818/1‑150 (MQ=255)                                                                                                                                           cGGCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCGGATcgcgcg     >  1:332659/1‑150 (MQ=255)                                                                                                                                            ggCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCGGATcgcgcg     >  1:815583/1‑149 (MQ=255)                                                                                                                                             gCAAAGTGCTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCGGATCGCGCgaga  <  2:208523/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               GCGAGTTCAGCGCCACCATCACCGCGCCGCTGCCTGCGTCCAGCCCCGCTTTGTACGGCGGCATATAATCATTAAACAGGCGCTGCGGACTCATATCGACGGTGTTGTACTCTTTACCGCCTTCTACCGCGCCGTATGCGGCAAAGTGTTTGACGCTGGTCATCACCGAGTAGCGATCTGCCGGGCTTTTACCCTGCATCGCTTCCACCATGGTTTTACCCATTGTTGAGGTGAGATACGTATCTTCGCCAAAACCTTCGGAAGCACGTCCCCAGCGCGGATCGCGCGAGA  >  NC_000913/2221218‑2221508

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