Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,521,967	C→A	100%	A171S (GCG→TCG)	xapR ←	transcriptional activator of xapAB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,521,967	0	C	A	100.0%	303.2 / NA	92	A171S (GCG→TCG)	xapR	transcriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (44/48);  total (44/48)

GGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCA  >  NC_000913/2521869‑2522062                                                                                                   |                                                                                                gggCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGcc                                                                                                <  1:2073688/100‑1 (MQ=255)   gcggcATTGTAACAAAATACTCGTAACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAAcgacga                                                                                              >  2:247278/1‑100 (MQ=255)   gcgccATTGTAACAAAATACCCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTCCCGGCATCGCCACCACAAAcgacga                                                                                              >  2:1069352/1‑100 (MQ=255)         ttGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATg                                                                                        >  1:2168916/1‑100 (MQ=255)          tGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTc                                                                                           >  2:34389/1‑96 (MQ=255)          tGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGc                                                                                       <  2:1559761/100‑1 (MQ=255)           gTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCa                                                                                      >  1:804881/1‑100 (MQ=255)            tAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCaa                                                                                     <  2:1813189/100‑1 (MQ=255)              aCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGCTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAAcc                                                                                   <  1:1800927/100‑1 (MQ=255)               cAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCg                                                                                  >  1:1142630/1‑100 (MQ=255)                 aaaTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTa                                                                                >  1:2113015/1‑100 (MQ=255)                 aaaTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCAACAGAAAAGACGATTCATGCAACCGTa                                                                                >  2:849745/1‑100 (MQ=255)                  aaTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTa                                                                                <  1:475802/99‑1 (MQ=255)                   aTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAg                                                                              <  2:1921102/100‑1 (MQ=255)                    tACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGc                                                                             <  1:2481634/100‑1 (MQ=255)                       tCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTgg                                                                          <  2:722825/100‑1 (MQ=255)                       tCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTgg                                                                          <  2:2284128/100‑1 (MQ=255)                        cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                         <  2:1332499/100‑1 (MQ=255)                        cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                         >  2:482041/1‑100 (MQ=255)                        cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                         >  1:586653/1‑100 (MQ=255)                        cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                         <  2:2406011/100‑1 (MQ=255)                        cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGt                                                                         <  1:1765341/100‑1 (MQ=255)                         gTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTa                                                                        <  1:2139410/100‑1 (MQ=255)                         gTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTa                                                                        <  2:1806083/100‑1 (MQ=255)                         gTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTa                                                                        >  1:943904/1‑100 (MQ=255)                          tCACGTAGCGCTTCCAGCGGGACGGTGGAAAATCATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTaa                                                                       >  2:574151/1‑100 (MQ=255)                             cGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAc                                                                    <  1:2458012/100‑1 (MQ=255)                             cGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAAc                                                                    <  1:2488070/100‑1 (MQ=255)                               tAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCa                                                                  <  2:928796/100‑1 (MQ=255)                               tAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCa                                                                  <  1:680575/100‑1 (MQ=255)                               tAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCa                                                                  <  1:680568/100‑1 (MQ=255)                                 gcgcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGt                                                                >  1:99032/1‑100 (MQ=255)                                  cgcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTGAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTc                                                               >  1:511182/1‑100 (MQ=255)                                  cgcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTc                                                               >  2:230422/1‑100 (MQ=255)                                  cgcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAAAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTc                                                               >  2:709916/1‑100 (MQ=255)                                   gcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCg                                                              >  2:898053/1‑100 (MQ=255)                                    cTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCgg                                                             <  2:297252/100‑1 (MQ=255)                                     ttCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggt                                                            <  2:1205990/100‑1 (MQ=255)                                      tCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtg                                                           >  1:666307/1‑100 (MQ=255)                                      tCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtg                                                           >  2:2311051/1‑100 (MQ=255)                                      tCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtg                                                           >  2:2362462/1‑100 (MQ=255)                                      tCCAGCGGGACGGTGGAAAATGATGAGAGACGATGCTCTTAAGGCATCGCCACCAGAAACGACGATTCATGCACCCGTAAGCTGGTAAAACCAGTCgctg                                                           >  2:2004813/1‑100 (MQ=255)                                        cAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtggt                                                         <  1:493279/100‑1 (MQ=255)                                        cAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtggt                                                         <  2:736809/100‑1 (MQ=255)                                        cAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCggtggt                                                         <  2:666300/100‑1 (MQ=255)                                         aGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGtt                                                        >  1:2446348/1‑100 (MQ=255)                                          gCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGtt                                                        >  1:877185/1‑99 (MQ=255)                                          gCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTc                                                       <  2:23929/100‑1 (MQ=255)                                           cGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCt                                                      <  1:1478025/100‑1 (MQ=255)                                            gggACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTg                                                     >  1:32070/1‑100 (MQ=255)                                              gACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTc                                                   <  2:81896/100‑1 (MQ=255)                                               aCGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCg                                                  <  2:2141408/100‑1 (MQ=255)                                                cGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGc                                                 <  1:2557275/100‑1 (MQ=255)                                                    ggAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCAtt                                             >  1:1171500/1‑100 (MQ=255)                                                     gAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTc                                            >  1:932983/1‑100 (MQ=255)                                                     gAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTc                                            <  2:1303064/100‑1 (MQ=255)                                                        aaTGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGcc                                         <  1:1948338/100‑1 (MQ=255)                                                         atgatgAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCa                                        <  1:186780/100‑1 (MQ=255)                                                             tgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATc                                    <  1:598002/100‑1 (MQ=255)                                                             tgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATc                                    >  2:612024/1‑100 (MQ=255)                                                                gagaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAGACCAGTCGGAGGTTCTGTCGCCATTCACCAGATCCCg                                 >  2:393805/1‑100 (MQ=255)                                                                 agaTGATGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCgg                                <  1:31855/100‑1 (MQ=255)                                                                  gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGc                               >  1:1788615/1‑100 (MQ=255)                                                                  gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGc                               <  1:1985969/100‑1 (MQ=255)                                                                  gatgatGCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGc                               >  1:692367/1‑100 (MQ=255)                                                                   atgatgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCa                              <  2:1053089/100‑1 (MQ=255)                                                                   atgatgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCa                              <  1:1339316/100‑1 (MQ=255)                                                                    tgatgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCAt                             >  2:1280462/1‑100 (MQ=255)                                                                     gatgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAAGCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCa                              <  2:516360/98‑1 (MQ=255)                                                                       tgCTCTTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCaa                          >  1:1881521/1‑100 (MQ=255)                                                                         ctctTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGt                        >  1:1093221/1‑100 (MQ=255)                                                                           ctTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTc                      <  1:2540552/100‑1 (MQ=255)                                                                           ctTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTc                      <  2:1704242/100‑1 (MQ=255)                                                                           ctTCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTc                      >  1:1578909/1‑100 (MQ=255)                                                                             tCAGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcgc                    >  1:2419142/1‑100 (MQ=255)                                                                               aGGCATCGCCACCAGAAAGGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcgg                  <  2:828136/100‑1 (MQ=255)                                                                               aGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcgg                  <  2:2101166/100‑1 (MQ=255)                                                                               aGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcgg                  <  1:31973/100‑1 (MQ=255)                                                                               aGGCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcgg                  <  1:1711677/100‑1 (MQ=255)                                                                                ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                 >  2:1616854/1‑100 (MQ=255)                                                                                ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                 <  1:1884725/100‑1 (MQ=255)                                                                                ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                 >  1:508697/1‑100 (MQ=255)                                                                                ggCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggc                 >  1:2235038/1‑100 (MQ=255)                                                                                 gCATCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCgcggcg                >  2:2462002/1‑100 (MQ=255)                                                                                    tCGCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTc             >  2:1700092/1‑100 (MQ=255)                                                                                      gCCACCAGAAACGACGATTCATGCAACGGTAAGCTGGTAAAACCAGTCCGTGGGTCTGTCTCCATCCGCCAGATCCCGGCATCAAGTTAGCGGCGttac            >  2:2218649/1‑97 (MQ=255)                                                                                      gCCACCAGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCa           >  1:2283583/1‑100 (MQ=255)                                                                                       ccaccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAg          <  1:2362477/100‑1 (MQ=255)                                                                                        caccaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGc         <  1:17473/100‑1 (MQ=255)                                                                                         accaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCa        >  2:1181861/1‑100 (MQ=255)                                                                                         accaGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCa        >  2:1254232/1‑100 (MQ=255)                                                                                           caGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAg      >  1:2198426/1‑100 (MQ=255)                                                                                           caGAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAg      <  1:881032/100‑1 (MQ=255)                                                                                             gAAACGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGc    <  1:1355343/100‑1 (MQ=255)                                                                                               aaCGACGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCa  >  2:1593242/1‑100 (MQ=255)                                                                                                   |                                                                                                GGGCGGCATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCA  >  NC_000913/2521869‑2522062

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