Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	284,450	A→G	E417G (GAG→GGG)	yagF →	CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	284,450	0	A	G	100.0%	263.1 / NA	75	E417G (GAG→GGG)	yagF	CP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/41);  total (34/41)

GTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGC  >  NC_000913/284302‑284595                                                                                                                                                     |                                                                                                                                                  gTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGc                                                                                                                                                 >  2:675548/1‑151 (MQ=255)     tGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCgg                                                                                                                                              >  2:700962/1‑150 (MQ=255)         ccACCTGCGGGTCCCCGGCCTGCTGCATCTGGCCGCCAGGGCCGTGACCGGCCAGCCGGTAGGCGAGAGCCTTGAATGGTGGCAGGCGTCGGAGCGCCGCAGGGGCTTCCGCCAGTGCCTGCGCGAGGAGGACGGCGTTGGGCCGGATGa                                                                                                                                          <  2:18158/150‑1 (MQ=255)          ccccTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCATAGGGCCGGATGACg                                                                                                                                        <  2:1152341/149‑1 (MQ=255)            ccTGCTCGACCTCGGCCTGCTGCATCTGGACGCCAAGGCCGTGGCCAGCCAGACGGTGGGCGAGACCCTTGAATGGAGGGAGGCGTCCGAACGCCGGGGACGCTTCCGCCAGTGCCTGCGCGGGCAGGACGGCGTAGGCGCGGATGACGTg                                                                                                                                      <  2:378694/151‑1 (MQ=255)              tGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGAt                                                                                                                                    >  1:1056239/1‑151 (MQ=255)              tGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGAt                                                                                                                                    >  2:494767/1‑151 (MQ=255)                  cgACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTg                                                                                                                                >  2:707794/1‑151 (MQ=255)                   gACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTg                                                                                                                                <  2:979737/150‑1 (MQ=255)                        cGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTgccgcc                                                                                                                           <  1:668665/150‑1 (MQ=255)                        cGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGccgccg                                                                                                                          >  1:787791/1‑151 (MQ=255)                         ggCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCgg                                                                                                                         <  1:415279/151‑1 (MQ=255)                           ccTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGgag                                                                                                                       <  1:95964/151‑1 (MQ=255)                             tgctgcATCTGGACGCCATGACCGTGACCGGCCAGAAGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCCGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTCCCGCCGgag                                                                                                                       >  2:660406/1‑149 (MQ=255)                              gctgcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAg                                                                                                                    <  1:1136337/151‑1 (MQ=255)                                 gcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCa                                                                                                                 <  1:432921/151‑1 (MQ=255)                                 gcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCa                                                                                                                 <  2:384894/151‑1 (MQ=255)                                     cTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCaca                                                                                                               >  2:183041/1‑147 (MQ=255)                                       ggACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCaaaa                                                                                                              <  2:1138669/148‑1 (MQ=255)                                       ggACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGc                                                                                                            <  2:772712/150‑1 (MQ=255)                                       ggACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCg                                                                                                           <  1:1241615/151‑1 (MQ=255)                                       ggACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCg                                                                                                           >  1:475106/1‑151 (MQ=255)                                               tgaccgtgaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCt                                                                                                   >  2:41101/1‑151 (MQ=255)                                               tgaccgtgaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCt                                                                                                   >  2:1069471/1‑151 (MQ=255)                                                  ccgtgaccgGCCAGGCGGTGGGCGAGAACCGGGAAGGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGa                                                                                                 <  2:85332/150‑1 (MQ=255)                                                   cgtgaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGAcc                                                                                               <  2:1212310/151‑1 (MQ=255)                                                        ccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGAc                                                                                          <  2:345802/151‑1 (MQ=255)                                                          ggCCAGGCGGTGGGCCGGAACCTTGAAGGGTGGCAGGCGTCGGAGCGCCGGGCGCGATTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACg                                                                                         <  2:46676/150‑1 (MQ=255)                                                               gACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTg                                                                                    <  2:154139/150‑1 (MQ=255)                                                                aCGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCt                                                                                  <  2:1056989/151‑1 (MQ=255)                                                                 cGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCtt                                                                                 >  1:968021/1‑151 (MQ=255)                                                                  ggTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCtt                                                                                 <  1:439166/150‑1 (MQ=255)                                                                    tGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTccc                                                                              >  1:653872/1‑151 (MQ=255)                                                                    tGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTccc                                                                              >  2:810259/1‑151 (MQ=255)                                                                     gggCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTccct                                                                             <  1:235424/151‑2 (MQ=255)                                                                            aaCCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGc                                                                       >  1:250930/1‑150 (MQ=255)                                                                                 tGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATc                                                                 >  2:1092337/1‑151 (MQ=255)                                                                                  gAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATc                                                                 <  2:1182156/150‑1 (MQ=255)                                                                                  gAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATc                                                                 >  1:823748/1‑150 (MQ=255)                                                                                   aaTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGc                                                               <  1:675675/151‑1 (MQ=255)                                                                                    aTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGc                                                               >  2:812220/1‑150 (MQ=255)                                                                                    aTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCt                                                              <  1:450765/151‑1 (MQ=255)                                                                                     tgggggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGAGGGCGTAGGGCCGGATGACGTTATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCt                                                              <  2:956261/150‑1 (MQ=255)                                                                                      ggtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTcc                                                            >  1:174717/1‑151 (MQ=255)                                                                                      ggtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTc                                                             <  2:534659/150‑1 (MQ=255)                                                                                       gtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCg                                                           <  1:1242715/151‑1 (MQ=255)                                                                                        tggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCgg                                                          >  1:958120/1‑151 (MQ=255)                                                                                         ggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGa                                                         <  2:49626/151‑1 (MQ=255)                                                                                           cAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGaa                                                        <  2:1056244/150‑1 (MQ=255)                                                                                                 tCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCgg                                                 >  2:451176/1‑151 (MQ=255)                                                                                                   cGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                               >  1:570029/1‑151 (MQ=255)                                                                                                   cGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                               <  2:1189107/151‑1 (MQ=255)                                                                                                    gAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                               >  1:374515/1‑150 (MQ=255)                                                                                                       cgcCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATc                                            >  2:1202839/1‑150 (MQ=255)                                                                                                        gcCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCa                                           >  1:958991/1‑150 (MQ=255)                                                                                                         ccGGGCGCGCTTCCGCCAGTGCCGGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGGCGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGAt                                             <  2:725925/147‑1 (MQ=255)                                                                                                          cGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAgg                                        >  1:1140502/1‑151 (MQ=255)                                                                                                           gggCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGGGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCCTCCCGACGGGCAACAACGCTCCGGAAGGTTCGGTGATCa                                           <  2:529546/147‑1 (MQ=255)                                                                                                           gggCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGc                                       >  1:160429/1‑151 (MQ=255)                                                                                                               gcgcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACg                                   >  1:226117/1‑151 (MQ=255)                                                                                                                   ttCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                >  2:550478/1‑150 (MQ=255)                                                                                                                   ttCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                >  2:729085/1‑150 (MQ=255)                                                                                                                    tCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCTCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                <  1:296776/149‑1 (MQ=255)                                                                                                                     ccGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATc                             <  1:1060217/151‑1 (MQ=255)                                                                                                                      cGCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTACGGAAGGTTCGGTGATCAAGGCCACGGCGa                               <  1:763986/148‑1 (MQ=255)                                                                                                                       gCCAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGa                           <  2:63128/151‑1 (MQ=255)                                                                                                                         cAGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAcc                         <  1:626121/151‑1 (MQ=255)                                                                                                                          aGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAccc                        >  2:525354/1‑151 (MQ=255)                                                                                                                          aGTGCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAccc                        <  1:1196145/151‑1 (MQ=255)                                                                                                                             gCCTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGt                      <  1:29701/150‑1 (MQ=255)                                                                                                                               cTGCGCGAGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCg                    <  1:183029/150‑1 (MQ=255)                                                                                                                                      aGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTggg            >  2:696677/1‑151 (MQ=255)                                                                                                                                      aGCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTggg            >  2:556438/1‑151 (MQ=255)                                                                                                                                       gCAGGACGGCGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGc           <  2:804158/151‑1 (MQ=255)                                                                                                                                                cGTAGGGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGc  >  1:137635/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGC  >  NC_000913/284302‑284595

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