Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,381	T→C	intergenic (‑325/‑45)	yigE ← / → corA	DUF2233 family protein/magnesium/nickel/cobalt transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,001,381	0	T	C	100.0%	303.8 / NA	82	intergenic (‑325/‑45)	yigE/corA	DUF2233 family protein/magnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/50);  total (32/50)

ACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCG  >  NC_000913/4001237‑4001527                                                                                                                                                 |                                                                                                                                                   aCTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCAc                                                                                                                                              <  2:522202/151‑1 (MQ=255)   tGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCg                                                                                                                                            >  1:887952/1‑151 (MQ=255)    ggCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCAcc                                                                                                                                             <  1:730461/149‑1 (MQ=255)    ggCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCg                                                                                                                                            <  2:1209297/150‑1 (MQ=255)    ggCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGa                                                                                                                                           <  1:1226535/151‑1 (MQ=255)     gCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGaa                                                                                                                                          <  2:1121408/151‑1 (MQ=255)     gCTCACTATCCTGTCATCCAGGATAAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGATGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGa                                                                                                                                           >  1:1247288/1‑150 (MQ=255)        cACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACt                                                                                                                                        <  1:694165/150‑1 (MQ=255)        cACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTCTCTGTGCGCCCACCGAACTg                                                                                                                                       <  1:1196585/151‑1 (MQ=255)            aTCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCg                                                                                                                                   <  1:1107907/151‑1 (MQ=255)              ccTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGat                                                                                                                                 <  1:418429/151‑1 (MQ=255)              ccTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGat                                                                                                                                 >  1:96461/1‑151 (MQ=255)                   cATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATAtttt                                                                                                                            >  1:510059/1‑151 (MQ=255)                   cATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATAtttt                                                                                                                            <  2:984074/151‑1 (MQ=255)                    aTCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTa                                                                                                                           <  2:265040/151‑1 (MQ=255)                      ccAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTaa                                                                                                                          <  1:567035/150‑1 (MQ=255)                      ccAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAg                                                                                                                         <  2:788402/151‑1 (MQ=255)                      ccAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAg                                                                                                                         <  1:924715/151‑1 (MQ=255)                        aGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGc                                                                                                                        >  1:683743/1‑150 (MQ=255)                        aGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCa                                                                                                                       <  2:223277/151‑1 (MQ=255)                           aTCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTg                                                                                                                    >  2:155783/1‑151 (MQ=255)                           aTCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCAGTCGCGCTGAGTCAGGCTGTTTAATGGTCTAAAACCCAATTTGTTTCTGTGCGCCCGCCGAACTGTCCGATATCTTAAGCAtt                                                                                                                     <  1:572682/150‑1 (MQ=255)                            tCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTgg                                                                                                                   >  2:115152/1‑151 (MQ=255)                             cAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTggg                                                                                                                  >  2:559986/1‑151 (MQ=255)                              aaCTCCTAAGGCTTTCCCTTTTCGCTGATAGCCATAGCGGTTGTCTGCGACCTCAATTTTTCCCGTCTCGCTGAGTCAGTCTTTTTAATGGTCTGAAACACAATTTGTTTCTGTGCGGCCACCGAACGGTCCGATATTTTGCGCATTgaga                                                                                                                 <  1:755911/151‑1 (MQ=255)                              aaCTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTggg                                                                                                                  >  1:992997/1‑150 (MQ=255)                                cTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGt                                                                                                               >  1:24423/1‑151 (MQ=255)                                  ccTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTc                                                                                                              <  2:1282986/150‑1 (MQ=255)                                    tAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCg                                                                                                           >  2:431777/1‑151 (MQ=255)                                     aaGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCg                                                                                                           <  2:641000/150‑1 (MQ=255)                                     aaGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCg                                                                                                           <  2:241276/150‑1 (MQ=255)                                      aGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGCGTCCCgg                                                                                                          >  1:160582/1‑150 (MQ=255)                                      aGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCgg                                                                                                          >  1:740394/1‑150 (MQ=255)                                           aTCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCAt                                                                                                      <  2:756852/149‑1 (MQ=255)                                           aTCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGc                                                                                                    >  2:451458/1‑151 (MQ=255)                                            tCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCt                                                                                                   >  2:247031/1‑151 (MQ=255)                                              ccTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGa                                                                                                 <  2:1272873/151‑1 (MQ=255)                                              ccTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGa                                                                                                 <  2:550543/151‑1 (MQ=255)                                              ccTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCACGCTGa                                                                                                 <  2:136973/151‑1 (MQ=255)                                                tttatGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCCGGCTGTTTAATGGTCTGAAACCCCATTTGTTTCTGGGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAg                                                                                                >  1:645956/1‑150 (MQ=255)                                                 ttttGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAgc                                                                                               <  2:46880/150‑1 (MQ=255)                                                    tGCTGATAGCCTTAGCGGTTGTCAGTGACCTCAATTTTTCCCGTCACGCTGAGTCATGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCCCCTAACTGTTCGATATTTTACGCATTGGGAGTCCCGGTCATGCTGAGCGCa                                                                                            <  1:7208/150‑1 (MQ=255)                                                            gCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTg                                                                                   <  2:641491/151‑1 (MQ=255)                                                                 aGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGaaa                                                                               >  1:1247986/1‑150 (MQ=255)                                                                    ggTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATaa                                                                           <  2:166539/151‑1 (MQ=255)                                                                     gTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAAc                                                                          <  2:735861/151‑1 (MQ=255)                                                                       tGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAAcc                                                                         <  2:569370/150‑1 (MQ=255)                                                                           aGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTg                                                                    >  1:590516/1‑151 (MQ=255)                                                                                   cAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTCCAACTGCCAACTCACCGCCTTACCCGGCt                                                            <  1:621722/151‑1 (MQ=255)                                                                                    aaTTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTg                                                           <  2:246397/151‑1 (MQ=255)                                                                                       ttttCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGa                                                         >  2:511889/1‑150 (MQ=255)                                                                                       ttttCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGa                                                         <  1:989536/150‑1 (MQ=255)                                                                                         ttCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAg                                                       <  2:1041132/150‑1 (MQ=255)                                                                                                cgcgCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAaga                                                <  2:241879/150‑1 (MQ=255)                                                                                                cgcgCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAag                                                 >  1:961829/1‑149 (MQ=255)                                                                                                 gcgcTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGt                                              <  2:1248023/151‑1 (MQ=255)                                                                                                   gcTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTc                                             >  2:1248266/1‑150 (MQ=255)                                                                                                    cTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTcac                                           >  1:154813/1‑151 (MQ=255)                                                                                                      gAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTcaca                                          >  1:91007/1‑150 (MQ=255)                                                                                                      gAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACaa                                         <  1:108645/151‑1 (MQ=255)                                                                                                        gTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAAcc                                       <  1:191581/151‑1 (MQ=255)                                                                                                            ggCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCtt                                   <  1:1268590/151‑1 (MQ=255)                                                                                                            ggCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCt                                    <  2:17224/150‑1 (MQ=255)                                                                                                             gCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTg                                  <  2:505234/151‑1 (MQ=255)                                                                                                                gTTTAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTa                                <  1:1074491/150‑1 (MQ=255)                                                                                                                   tAATGGTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATg                            >  2:1277532/1‑151 (MQ=255)                                                                                                                        gTCTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGt                        >  1:259112/1‑150 (MQ=255)                                                                                                                          cTGAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATg                     >  1:856395/1‑151 (MQ=255)                                                                                                                            gAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATgg                    >  1:1283860/1‑150 (MQ=255)                                                                                                                            gAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATgg                    >  1:510428/1‑150 (MQ=255)                                                                                                                            gAAACCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATg                     <  2:1294846/149‑1 (MQ=255)                                                                                                                               aCCCAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTg                <  1:15342/151‑1 (MQ=255)                                                                                                                                 ccAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGAt              <  1:1284735/151‑1 (MQ=255)                                                                                                                                  cAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATc             <  1:611993/151‑1 (MQ=255)                                                                                                                                  cAATTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATc             >  2:416330/1‑151 (MQ=255)                                                                                                                                   aaTTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCt            >  2:906086/1‑151 (MQ=255)                                                                                                                                    aTTTGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCtt           <  1:411489/151‑1 (MQ=255)                                                                                                                                     tttGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTg          >  2:279926/1‑151 (MQ=255)                                                                                                                                       tGTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGt         <  1:600403/150‑1 (MQ=255)                                                                                                                                        gTTTCTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTc        <  1:799045/150‑1 (MQ=255)                                                                                                                                            cTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCTGTATGAATTGATCTTGTCGaa     <  1:158575/149‑1 (MQ=255)                                                                                                                                            cTGTGCGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAAcc   >  2:215220/1‑151 (MQ=255)                                                                                                                                             tgtgcgCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCg  >  2:1060926/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                   ACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCAATTTGTTTCTGTGTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCG  >  NC_000913/4001237‑4001527

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