Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	613,412	C→T	P226S (CCC→TCC)	fes →	enterobactin/ferrienterobactin esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	613,412	0	C	T	100.0%	244.9 / NA	64	P226S (CCC→TCC)	fes	enterobactin/ferrienterobactin esterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/32);  total (32/32)

AAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGC  >  NC_000913/613264‑613558                                                                                                                                                     |                                                                                                                                                   aaaGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACttccttcc                                                                                                                                                  <  2:527157/151‑1 (MQ=255)  aaGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTccg                                                                                                                                                 <  1:583854/151‑1 (MQ=255)      cACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTccgcc                                                                                                                                               >  2:365782/1‑149 (MQ=255)        aggCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCCCGATGGCGAATTTTGGGCGCAAAGAATGCCCGCCTGGCCAGTGCTGAATTCGCTGACCCATCGTCAGCAACTTCCTTccgccg                                                                                                                                              >  2:657033/2‑148 (MQ=255)          gCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTAt                                                                                                                                         <  2:318169/151‑1 (MQ=255)            gtgtATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTAtgt                                                                                                                                       >  1:902503/1‑151 (MQ=255)              gtATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATgtgt                                                                                                                                     <  2:871354/151‑1 (MQ=255)                aTGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTg                                                                                                                                   >  1:219343/1‑151 (MQ=255)                    aTTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCg                                                                                                                               >  1:586786/1‑151 (MQ=255)                     tttttACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCg                                                                                                                               >  1:162991/1‑150 (MQ=255)                     tttttACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGCGTTGATCg                                                                                                                               >  1:331376/1‑150 (MQ=255)                       tttACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACg                                                                                                                            >  2:140783/1‑151 (MQ=255)                           ccaccGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTAt                                                                                                                        <  2:600386/151‑1 (MQ=255)                                    aTGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGAcaccac                                                                                                               >  2:754617/1‑151 (MQ=255)                                    aTGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGAcaccac                                                                                                               <  1:234628/151‑1 (MQ=255)                                       cAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCa                                                                                                            <  2:168471/151‑1 (MQ=255)                                           agcagAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgc                                                                                                        >  2:808699/1‑151 (MQ=255)                                            gcagAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg                                                                                                       >  2:605263/1‑151 (MQ=255)                                              agaagaACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGcc                                                                                                     >  2:246418/1‑151 (MQ=255)                                                  gaaCGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACg                                                                                                 <  2:166109/151‑1 (MQ=255)                                                   aaCGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACg                                                                                                 >  2:206649/1‑150 (MQ=255)                                                     cGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAAc                                                                                              <  2:148792/151‑1 (MQ=255)                                                     cGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAAc                                                                                              <  1:260968/151‑1 (MQ=255)                                                       cccGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTg                                                                                            <  2:512877/151‑1 (MQ=255)                                                        ccGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGc                                                                                           >  1:723673/1‑151 (MQ=255)                                                         cGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGcc                                                                                          <  2:586736/151‑1 (MQ=255)                                                           cTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGCACTGcc                                                                                          >  2:230156/1‑149 (MQ=255)                                                                aGTTTTGCTCGATGACGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAAt                                                                                   >  1:856254/1‑151 (MQ=255)                                                                   tttGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCg                                                                                >  1:216420/1‑151 (MQ=255)                                                                      gCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGa                                                                              <  1:787552/150‑1 (MQ=255)                                                                      gCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGAt                                                                             >  2:868521/1‑151 (MQ=255)                                                                      gCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGAt                                                                             >  1:57524/1‑151 (MQ=255)                                                                         cGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTc                                                                          >  1:540328/1‑151 (MQ=255)                                                                          gATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACGGCCGTGTAATGCGGATTTCt                                                                         <  2:466739/151‑1 (MQ=255)                                                                              gCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGc                                                                      >  2:698838/1‑150 (MQ=255)                                                                                       gggCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACa                                                            <  1:873000/151‑1 (MQ=255)                                                                                        ggCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAg                                                           <  1:266381/151‑1 (MQ=255)                                                                                         gcgcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGc                                                          <  1:129603/151‑1 (MQ=255)                                                                                         gcgcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACGGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGc                                                          <  2:398672/151‑1 (MQ=255)                                                                                          cgcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCa                                                         <  2:892662/151‑1 (MQ=255)                                                                                           gcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCa                                                         >  2:619171/1‑150 (MQ=255)                                                                                            cAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAag                                                       <  1:221110/151‑1 (MQ=255)                                                                                               aGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGt                                                    <  1:48441/151‑1 (MQ=255)                                                                                                 tATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGtt                                                   >  1:289972/1‑150 (MQ=255)                                                                                                    gCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGttatta                                               <  2:112304/151‑1 (MQ=255)                                                                                                    gCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGACAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACGCCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGttatta                                               <  2:872824/151‑1 (MQ=255)                                                                                                       cGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTAccc                                            <  1:451765/151‑1 (MQ=255)                                                                                                        gTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTAcccc                                           <  1:870882/151‑1 (MQ=255)                                                                                                        gTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTAcccc                                           >  2:617920/1‑151 (MQ=255)                                                                                                           tGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTgg                                        >  1:596016/1‑151 (MQ=255)                                                                                                            ggCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTg                                         <  1:273804/149‑1 (MQ=255)                                                                                                                aGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGaaa                                   <  2:241083/151‑1 (MQ=255)                                                                                                                   gCTGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCt                                <  1:49393/151‑1 (MQ=255)                                                                                                                     tGACTTCGCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTAt                              >  2:135439/1‑151 (MQ=255)                                                                                                                            gCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAg                                                                <  2:2085/110‑1 (MQ=255)                                                                                                                            gCTGACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAg                                                                >  1:2085/1‑110 (MQ=255)                                                                                                                               gACCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCttt                     <  2:339403/150‑1 (MQ=255)                                                                                                                                aCCCATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTa                   <  2:349850/151‑1 (MQ=255)                                                                                                                                   cATCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCg                >  2:522337/1‑151 (MQ=255)                                                                                                                                    aTCGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCg                >  1:846470/1‑150 (MQ=255)                                                                                                                                      cGTCAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATc             >  2:281341/1‑151 (MQ=255)                                                                                                                                         cAGCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGt           <  1:49540/150‑1 (MQ=255)                                                                                                                                           gCAACTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGc         <  2:296333/150‑1 (MQ=255)                                                                                                                                             aaCTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCg       >  1:635841/1‑150 (MQ=255)                                                                                                                                             aaCTTCCTTCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCg       >  1:129584/1‑150 (MQ=255)                                                                                                                                                ttccttccGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCg   >  2:514838/1‑151 (MQ=255)                                                                                                                                                 tccttcCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGc  >  2:747265/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AAAGTCACGGCGTGTATGGATTTTTACCACCGGCGATGCAACAGCAGAAGAACGCCCGCTGGCAGTTTTGCTCGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAACTGCCGTGTAATGCGGATTTCTGGCTCGCAGTACAGCAAGAGTTATTACCCCTGGTGAAAGCTATTGCCCCTTTTAGCGATCGTGCCGATCGC  >  NC_000913/613264‑613558

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