Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,403,176	C→T	P380L (CCG→CTG)	hflK →	modulator for HflB protease specific for phage lambda cII repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,403,176	0	C	T	100.0%	213.8 / NA	56	P380L (CCG→CTG)	hflK	modulator for HflB protease specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/26);  total (30/26)

ATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTATC  >  NC_000913/4403028‑4403317                                                                                                                                                     |                                                                                                                                              aTGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCtgctgc                                                                                                                                             >  1:882595/1‑151 (MQ=255)  tGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgctgcc                                                                                                                                            >  2:290919/1‑151 (MQ=255)   ggAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgctgcc                                                                                                                                            <  1:302675/150‑1 (MQ=255)    gAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgctgcc                                                                                                                                            >  1:279309/1‑149 (MQ=255)      aaaaGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGcc                                                                                                                                        >  1:653292/1‑151 (MQ=255)       aaaGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCct                                                                                                                                       >  1:443332/1‑151 (MQ=255)           tgtTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTtc                                                                                                                                   >  1:645852/1‑151 (MQ=255)                 gTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCcac                                                                                                                             <  2:101855/151‑1 (MQ=255)                     cacCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAAcc                                                                                                                         >  2:481743/1‑151 (MQ=255)                      acCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAAcc                                                                                                                         >  2:332654/1‑150 (MQ=255)                      acCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAAcc                                                                                                                         >  2:584207/1‑150 (MQ=255)                       cccGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCACTTCCTCCACAACCa                                                                                                                        <  2:621109/150‑1 (MQ=255)                          gCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTgg                                                                                                                    <  2:120073/151‑1 (MQ=255)                          gCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTgg                                                                                                                    >  2:203287/1‑151 (MQ=255)                              aGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGagca                                                                                                                <  1:573984/151‑1 (MQ=255)                                  cTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                                                                                             >  2:80048/1‑150 (MQ=255)                                   tGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGagcaagca                                                                                                            >  1:401972/1‑150 (MQ=255)                                     gTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAATGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAac                                                                                                          >  1:656558/1‑150 (MQ=255)                                        aaCGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgt                                                                                                      >  2:336972/1‑151 (MQ=255)                                         aCGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAAc                                                                                                                          <  2:23052/130‑1 (MQ=255)                                         aCGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAAc                                                                                                                          >  1:23074/1‑130 (MQ=255)                                             tAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcc                                                                                                 >  2:749936/1‑151 (MQ=255)                                              aaaGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcc                                                                                                 >  1:575242/1‑150 (MQ=255)                                              aaaGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcc                                                                                                 <  1:598184/150‑1 (MQ=255)                                               aaGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTccac                                                                                               <  2:527118/151‑1 (MQ=255)                                                   tGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGt                                                                                           >  1:535758/1‑151 (MQ=255)                                                      cAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTc                                                                                          >  2:112692/1‑149 (MQ=255)                                                      cAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTc                                                                                          >  1:194161/1‑149 (MQ=255)                                                       aaCCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAgg                                                                                       <  1:435788/151‑1 (MQ=255)                                                       aaCCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAgg                                                                                       <  2:585787/151‑1 (MQ=255)                                                          cTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCg                                                                                    >  2:402935/1‑151 (MQ=255)                                                             aTGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGat                                                                                  >  2:108417/1‑150 (MQ=255)                                                              tGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGatat                                                                                >  2:832309/1‑151 (MQ=255)                                                                          gaGAACAGATGCTGACAGGGGGTAACCCCCCTGCGGCGAAGAGAGATAACGGTGCCAGCAATCTGCGGCGTCTGCTGTGAGCGTCTTCATACACAAAAAGTGGAGCAAGCAACACGACGTACACCAGCCAGGCGGATATAATGGACCa                                                                       >  2:549853/2‑148 (MQ=255)                                                                           aGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAAcgc                                                                   <  2:235306/151‑1 (MQ=255)                                                                            gACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAAcgcc                                                                  <  1:286822/151‑1 (MQ=255)                                                                               cAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCcgcg                                                               >  1:728325/1‑151 (MQ=255)                                                                                    gCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAAc                                                          <  1:342536/151‑1 (MQ=255)                                                                                        aaaGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACgcgc                                                      >  1:478076/1‑151 (MQ=255)                                                                                         aaGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCa                                                     <  1:687425/151‑1 (MQ=255)                                                                                              ggtAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGt                                                 >  1:24539/1‑150 (MQ=255)                                                                                                        cTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTAcc                                       >  2:636213/1‑150 (MQ=255)                                                                                                          gcggcgAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCa                                      <  2:350683/149‑1 (MQ=255)                                                                                                          gcggcgAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGc                                    >  2:704806/1‑151 (MQ=255)                                                                                                            ggcgAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGt                                  >  1:4107/1‑151 (MQ=255)                                                                                                             gcgAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTc                                 >  1:385486/1‑151 (MQ=255)                                                                                                                 agagCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAggg                             >  1:863270/1‑151 (MQ=255)                                                                                                                  gagCGATAACGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAgggg                            <  1:805552/151‑1 (MQ=255)                                                                                                                         aaCGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAAc                     <  1:247498/151‑1 (MQ=255)                                                                                                                         aaCGGTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAAc                     <  2:138959/151‑1 (MQ=255)                                                                                                                             gTGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATg                 <  1:723826/151‑1 (MQ=255)                                                                                                                              tGCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGc                >  2:533108/1‑151 (MQ=255)                                                                                                                               gCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCg               <  2:401946/151‑1 (MQ=255)                                                                                                                               gCCAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCg               <  1:753808/151‑1 (MQ=255)                                                                                                                                 cAGCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTa             <  2:342966/151‑1 (MQ=255)                                                                                                                                   gCAATCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAg           <  1:189613/151‑1 (MQ=255)                                                                                                                                     aaTCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTc         <  1:112714/151‑1 (MQ=255)                                                                                                                                     aaTCTGCTGCGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTc         <  2:194136/151‑1 (MQ=255)                                                                                                                                          gctgcGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTa    <  1:834135/151‑1 (MQ=255)                                                                                                                                             gcGTCTGCTGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCTCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTATc  <  2:777821/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              ATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCGCCAACGCGCAGCGTAACGACTACCAGCGTCAGGGGGAATAACGATGCGTAAGTCAGTTATC  >  NC_000913/4403028‑4403317

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