Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,208,534	Δ1 bp	coding (403/414 nt);coding (586/603 nt)	tfaP → tfaE ←	e14 prophage, uncharacterized protein e14 prophage, putative tail fiber assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,208,534	0	C	.	100.0%	300.1 / NA	77	coding (403/414 nt) coding (586/603 nt)	tfaP tfaE	e14 prophage, uncharacterized protein e14 prophage, putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (39/38);  total (39/38)

TAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGA  >  NC_000913/1208401‑1208667                                                                                                                                      |                                                                                                                                      tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCa                                                                                                                                 >  2:1775596/1‑139 (MQ=255)   aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGa                                                                                                                               >  2:315925/1‑139 (MQ=255)              ggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgc                                                                                                                    >  1:1499660/1‑139 (MQ=255)               gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  2:2013542/139‑1 (MQ=255)               gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  2:1499660/139‑1 (MQ=255)                    cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                              <  1:1760001/139‑1 (MQ=255)                              aTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTg                                                                                                    >  2:1476965/1‑139 (MQ=255)                               tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGt                                                                                                   <  2:2032510/139‑1 (MQ=255)                                            tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                                      >  2:66383/1‑139 (MQ=255)                                            tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                                      >  2:14484/1‑139 (MQ=255)                                             ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                     <  1:1838666/139‑1 (MQ=255)                                             ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                     <  2:731319/139‑1 (MQ=255)                                             ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCACCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                     <  2:1095209/139‑1 (MQ=255)                                                aTTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCaa                                                                                  >  1:845138/1‑139 (MQ=255)                                                  ttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                <  2:845138/139‑1 (MQ=255)                                                    gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTa                                                                              >  1:1094155/1‑139 (MQ=255)                                                      ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                                            >  1:209008/1‑139 (MQ=255)                                                        cTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  1:1961327/1‑98 (MQ=255)                                                        cTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  2:1961327/98‑1 (MQ=255)                                                           gATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                                       >  1:1431282/1‑139 (MQ=255)                                                            aTTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                                                      <  1:1476965/139‑1 (MQ=255)                                                             ttATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                                                     <  2:89513/139‑1 (MQ=255)                                                             ttATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                                                     <  2:1094155/139‑1 (MQ=255)                                                                     aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                                               <  2:573587/137‑1 (MQ=255)                                                                     aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                                               >  1:573587/1‑137 (MQ=255)                                                                           tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                       <  1:458272/139‑1 (MQ=255)                                                                           tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                       >  2:847292/1‑139 (MQ=255)                                                                               aCTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGgttgt                                                   <  1:14484/139‑1 (MQ=255)                                                                                   gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcaca                                                                    >  2:1812147/1‑118 (MQ=255)                                                                                   gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcaca                                                                    <  1:1812147/118‑1 (MQ=255)                                                                                   gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTct                                               >  2:284292/1‑139 (MQ=255)                                                                                       aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                                       >  1:47272/1‑111 (MQ=255)                                                                                       aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                                       <  2:47272/111‑1 (MQ=255)                                                                                        tACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  1:461845/66‑1 (MQ=14)                                                                                        tACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  2:461845/1‑66 (MQ=14)                                                                                           aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcg                                                                                        <  1:649165/90‑1 (MQ=255)                                                                                           aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcg                                                                                        >  2:649165/1‑90 (MQ=255)                                                                                            tCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTc                                      <  1:284292/139‑1 (MQ=255)                                                                                             ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGt                                                           <  1:676833/117‑1 (MQ=255)                                                                                             ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGt                                                           >  2:676833/1‑117 (MQ=255)                                                                                             ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                     >  2:2041002/1‑139 (MQ=255)                                                                                             ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                     >  1:1915015/1‑139 (MQ=255)                                                                                                gTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGt                                   <  2:1381694/138‑1 (MQ=255)                                                                                                gTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGt                                   >  1:1381694/1‑138 (MQ=255)                                                                                                 tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                       >  2:1422069/1‑133 (MQ=255)                                                                                                 tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                       <  1:1422069/133‑1 (MQ=255)                                                                                                 tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt                                 >  2:1481645/1‑139 (MQ=255)                                                                                                 tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt                                 >  1:1319966/1‑139 (MQ=255)                                                                                                  gcgcCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                              >  2:1491428/1‑61 (MQ=11)                                                                                                  gcgcCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                              <  1:1491428/61‑1 (MQ=11)                                                                                                  gcgcCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTa                                >  2:1094887/1‑139 (MQ=255)                                                                                                      cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc                            <  1:2041002/139‑1 (MQ=255)                                                                                                         atatTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCgg                         <  2:1915015/139‑1 (MQ=255)                                                                                                           atTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                <  2:1437710/82‑1 (MQ=255)                                                                                                           atTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                >  1:1437710/1‑82 (MQ=255)                                                                                                             tGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTg                     <  2:178784/139‑1 (MQ=255)                                                                                                                aTGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                >  1:49945/1‑77 (MQ=11)                                                                                                                aTGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                <  2:49945/77‑1 (MQ=11)                                                                                                                   gcttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  2:1274595/36‑1 (MQ=14)                                                                                                                   gcttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  1:1274595/4‑39 (MQ=14)                                                                                                                    cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc              >  2:73359/3‑139 (MQ=255)                                                                                                                    cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc              >  2:1744446/3‑139 (MQ=255)                                                                                                                    cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc              >  2:150628/3‑139 (MQ=255)                                                                                                                     ttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATc                                                                           <  2:255138/76‑1 (MQ=255)                                                                                                                     ttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATc                                                                           >  1:255138/2‑77 (MQ=255)                                                                                                                      tACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAg            >  2:183042/1‑139 (MQ=255)                                                                                                                      tACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAg            <  1:183042/139‑1 (MQ=255)                                                                                                                       aCGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGt           <  1:1775596/139‑1 (MQ=255)                                                                                                                            tCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAAt      <  1:1744446/139‑1 (MQ=255)                                                                                                                             ccGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATc     >  2:947641/1‑139 (MQ=255)                                                                                                                             ccGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATc     >  2:1864061/1‑139 (MQ=255)                                                                                                                              cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTGATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                                             <  2:877797/98‑1 (MQ=255)                                                                                                                              cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCt                                                         >  2:638791/1‑86 (MQ=255)                                                                                                                              cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCt                                                         <  1:638791/86‑1 (MQ=255)                                                                                                                              cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                                             >  1:877797/1‑98 (MQ=255)                                                                                                                              cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCa    <  1:847292/139‑1 (MQ=255)                                                                                                                                gCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGCTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGa  <  2:1319966/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      TAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGA  >  NC_000913/1208401‑1208667

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