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%	376.5 / NA	100	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 . (56/44);  total (56/44)

AGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGA  >  NC_000913/1208403‑1208667                                                                                                                                    |                                                                                                                                      aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGa                                                                                                                               >  1:1276405/1‑139 (MQ=255)       cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACAt                                                                                                                         >  2:273842/1‑139 (MQ=255)       cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACAt                                                                                                                         >  1:232388/1‑139 (MQ=255)             gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  1:1523008/1‑139 (MQ=255)                 aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGt                                                                                                               <  2:1255207/139‑1 (MQ=255)                 aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGt                                                                                                               >  2:735332/1‑139 (MQ=255)                  cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                              <  1:273842/139‑1 (MQ=255)                           gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                     >  2:546847/1‑139 (MQ=255)                           gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                     >  1:1766574/1‑139 (MQ=255)                               aGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  2:1495242/121‑1 (MQ=255)                               aGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  1:1495242/1‑121 (MQ=255)                                gCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGc                                                                                                >  1:1236970/1‑139 (MQ=255)                                   gAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGt                                                                                             >  1:844564/1‑139 (MQ=255)                                          tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                                      >  2:1747757/1‑139 (MQ=255)                                          tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                                      >  1:1268869/1‑139 (MQ=255)                                          tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                                      >  1:529696/1‑139 (MQ=255)                                           ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                                     <  2:1523008/139‑1 (MQ=255)                                              aTTTGTGGCTGGATTATCTGGAAGCACCGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCaa                                                                                  >  1:1212665/1‑139 (MQ=255)                                               tttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAAt                                                                                 >  2:1664242/1‑139 (MQ=255)                                               tttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAAt                                                                                 >  1:1185756/1‑139 (MQ=255)                                                ttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                <  2:1236970/139‑1 (MQ=255)                                                    ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                                            >  2:1074611/1‑139 (MQ=255)                                                     gCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  2:561906/1‑99 (MQ=255)                                                     gCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  1:561906/99‑1 (MQ=255)                                                     gCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATc                                                                           >  1:688881/1‑139 (MQ=255)                                                      cTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGc                                                                                                             >  1:703432/1‑104 (MQ=255)                                                      cTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGc                                                                                                             <  2:703432/104‑1 (MQ=255)                                                         gATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGc                                                                       >  2:1466719/1‑139 (MQ=255)                                                          aTTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                                                      <  2:1185756/139‑1 (MQ=255)                                                          aTTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                                                      <  1:735332/139‑1 (MQ=255)                                                           ttATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                                                     <  1:769916/139‑1 (MQ=255)                                                              tCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  2:428996/1‑90 (MQ=18)                                                              tCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  1:428996/90‑1 (MQ=18)                                                                   aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTaa                                                             >  1:209964/1‑139 (MQ=255)                                                                       aCTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                <  1:610471/116‑1 (MQ=255)                                                                       aCTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                >  2:610471/1‑116 (MQ=255)                                                                        cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  2:133532/80‑1 (MQ=17)                                                                        cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  2:1517529/1‑80 (MQ=17)                                                                        cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  1:1517529/80‑1 (MQ=17)                                                                        cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  1:133532/1‑80 (MQ=17)                                                                        cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                                 >  2:798477/1‑82 (MQ=17)                                                                        cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGcgcg                                                                                                                 <  1:798477/82‑1 (MQ=17)                                                                        cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTg                                                        >  1:1354209/1‑139 (MQ=255)                                                                         tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                       <  1:1307541/139‑1 (MQ=255)                                                                         tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                       >  2:1699940/1‑139 (MQ=255)                                                                                ggTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTc                                                >  2:393815/1‑139 (MQ=255)                                                                                 gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                              <  1:496656/76‑1 (MQ=14)                                                                                 gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTg                                                                                                              >  2:496656/1‑76 (MQ=14)                                                                                 gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                       <  2:317117/131‑1 (MQ=255)                                                                                 gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                       >  1:317117/1‑131 (MQ=255)                                                                                 gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTct                                               >  1:485799/1‑139 (MQ=255)                                                                                    gATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACGTCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                            <  2:864514/139‑1 (MQ=255)                                                                                    gATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                            <  2:1212665/139‑1 (MQ=255)                                                                                     aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                <  2:329487/102‑1 (MQ=255)                                                                                     aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                                                >  1:329487/1‑102 (MQ=255)                                                                                     aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcct                                           >  2:507/1‑139 (MQ=255)                                                                                         aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                       <  2:529696/139‑1 (MQ=255)                                                                                         aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                       <  1:507/139‑1 (MQ=255)                                                                                         aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                       <  1:225517/139‑1 (MQ=255)                                                                                         aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                                       <  1:1466719/139‑1 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTaa                                                                             <  2:1414343/99‑1 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTaa                                                                             >  1:1414343/1‑99 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                                               >  2:1087283/1‑113 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                                               <  1:1087283/113‑1 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                            >  1:71603/1‑132 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                            <  2:71603/132‑1 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                     >  2:430424/1‑139 (MQ=255)                                                                                           ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                     >  2:501007/1‑139 (MQ=255)                                                                                            cAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                                                     >  2:1316149/1‑106 (MQ=255)                                                                                            cAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGcac                                                                     <  1:1316149/106‑1 (MQ=255)                                                                                               tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATg                                                                                                                             >  1:1079207/1‑47 (MQ=14)                                                                                               tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATg                                                                                                                             <  2:1079207/47‑1 (MQ=14)                                                                                               tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   <  1:1054347/57‑1 (MQ=11)                                                                                               tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGgcg                                                                                                                   >  2:1054347/1‑57 (MQ=11)                                                                                               tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt                                 >  2:215987/1‑139 (MQ=255)                                                                                                 cgcCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                                                       >  2:97789/1‑67 (MQ=11)                                                                                                 cgcCAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                                                       <  1:97789/67‑1 (MQ=11)                                                                                                   ccAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGt                             >  1:1441062/1‑139 (MQ=255)                                                                                                    cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc                            <  2:1441062/139‑1 (MQ=255)                                                                                                    cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc                            >  1:97654/1‑139 (MQ=255)                                                                                                    cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc                            >  2:135243/1‑139 (MQ=255)                                                                                                    cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc                            >  2:1699207/1‑139 (MQ=255)                                                                                                    cAGATATTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc                            >  2:135244/1‑139 (MQ=255)                                                                                                        tatTGAATGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGc                        <  1:1664242/139‑1 (MQ=255)                                                                                                               tGGCTTACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAt                 <  2:97654/139‑1 (MQ=255)                                                                                                                 gcttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                       <  1:911973/96‑1 (MQ=255)                                                                                                                 gcttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                                       >  2:911973/4‑99 (MQ=255)                                                                                                                  cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                     <  2:363338/50‑1 (MQ=14)                                                                                                                  cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                                                                                                     >  1:363338/3‑52 (MQ=11)                                                                                                                  cttACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc              >  2:791322/3‑139 (MQ=255)                                                                                                                    tACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAg            <  1:393815/139‑1 (MQ=255)                                                                                                                    tACGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAg            <  1:1699207/139‑1 (MQ=255)                                                                                                                     aCGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGt           <  1:430424/139‑1 (MQ=255)                                                                                                                     aCGCCTCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGt           <  1:1699940/139‑1 (MQ=255)                                                                                                                          tCCGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAAt      <  2:844564/139‑1 (MQ=255)                                                                                                                           ccGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                                             >  1:387192/1‑99 (MQ=255)                                                                                                                           ccGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                                             <  2:387192/99‑1 (MQ=255)                                                                                                                           ccGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATc     >  2:101803/1‑139 (MQ=255)                                                                                                                            cGGCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCa    >  1:405222/1‑139 (MQ=255)                                                                                                                              gCAGTT‑AGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGa  <  1:546847/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      AGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGA  >  NC_000913/1208403‑1208667

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