Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	612,932	C→A	Q66K (CAG→AAG)	fes →	enterochelin esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	612,932	0	C	A	100.0%	263.3 / NA	74	Q66K (CAG→AAG)	fes	enterochelin esterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (37/37);  total (37/37)

GCAGCGTCTGAATGACGAAATGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCCAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTC  >  NC_000913/612796‑613066                                                                                                                                         |                                                                                                                                       gCAGCGTCTGAATGACGAAATGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGt                                                                                                                                     <  2:666593/140‑1 (MQ=255)         tGAATGACGAAATGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAg                                                                                                                             <  2:664637/140‑1 (MQ=255)          gAATGACGAAATGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGc                                                                                                                            <  2:1240630/140‑1 (MQ=255)           aaTGACGAAATGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGt                                                                                                                                     <  1:1248313/130‑1 (MQ=255)           aaTGACGAAATGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGt                                                                                                                                     >  2:1248313/1‑130 (MQ=255)                 gAAATGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGc                                                                                                                     <  1:1141235/140‑1 (MQ=255)                   aaTGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAg                                                                                                                   >  1:525885/1‑140 (MQ=255)                         tGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGc                                                                                                                            >  1:144928/1‑125 (MQ=255)                         tGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGc                                                                                                                            <  2:144928/125‑1 (MQ=255)                             gTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGt                                                                                                         <  1:699924/140‑1 (MQ=255)                                 cTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTgg                                                                                                     >  2:1311853/1‑140 (MQ=255)                                        gtggCGTGTTCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGga                                                                                              >  2:1073594/1‑140 (MQ=255)                                        gtggCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGga                                                                                              >  2:1231043/1‑140 (MQ=255)                                        gtggCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGga                                                                                              <  1:1231043/140‑1 (MQ=255)                                           gCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGgacga                                                                                           <  1:1029095/140‑1 (MQ=255)                                            cGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCa                                                                                                   <  2:806498/131‑1 (MQ=255)                                            cGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCa                                                                                                   >  1:806498/1‑131 (MQ=255)                                            cGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGgacgac                                                                                          >  2:5428/1‑140 (MQ=255)                                            cGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCG‑‑TGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGacac                                                                                        >  1:823444/1‑140 (MQ=255)                                             gTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGaca                                                                                         <  2:1354111/140‑1 (MQ=255)                                              tGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGacac                                                                                        >  1:643878/1‑140 (MQ=255)                                                aTCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCG‑‑TGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACt                                                                                    <  2:823444/140‑1 (MQ=255)                                                  ccccAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACt                                                                                    >  2:754244/1‑140 (MQ=255)                                                  ccccAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACt                                                                                    <  2:790036/140‑1 (MQ=255)                                                    ccAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCa                                                                                  <  2:654258/140‑1 (MQ=255)                                                          ttCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTTGACGACACAACTCAATGCCa                                                                            >  2:1062696/1‑140 (MQ=255)                                                            cTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAAc                                                                          >  2:670752/1‑140 (MQ=255)                                                            cTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAAc                                                                          >  2:122532/1‑140 (MQ=255)                                                             tGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACt                                                                         >  2:933947/1‑140 (MQ=255)                                                                agaaTACTCGACGACAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGgc                                                                      <  1:1311853/140‑1 (MQ=255)                                                                    tACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCgg                                                                  >  2:10291/1‑140 (MQ=255)                                                                      cTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCa                                                                >  1:126966/1‑140 (MQ=255)                                                                      cTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCa                                                                >  1:336171/1‑140 (MQ=255)                                                                         gacgaTAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCt                                                             >  2:1255726/1‑140 (MQ=255)                                                                          acgaTAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTa                                                            <  1:10291/140‑1 (MQ=255)                                                                          acgaTAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTa                                                            <  1:1073594/140‑1 (MQ=255)                                                                            gaTAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACt                                                          <  2:313163/140‑1 (MQ=255)                                                                            gaTAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACt                                                          <  1:933947/140‑1 (MQ=255)                                                                              tAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGc                                                        >  1:43367/1‑140 (MQ=255)                                                                               aaaGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCt                                                       <  1:1255726/140‑1 (MQ=255)                                                                                   cgcgTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACt                                                          >  2:327834/1‑133 (MQ=255)                                                                                   cgcgTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACt                                                          <  1:327834/133‑1 (MQ=255)                                                                                       tATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTa                                                            >  2:1058733/1‑127 (MQ=255)                                                                                       tATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTa                                                            <  1:1058733/127‑1 (MQ=255)                                                                                       tATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTccc                                               >  2:848407/1‑140 (MQ=255)                                                                                       tATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTccc                                               >  1:839362/1‑140 (MQ=255)                                                                                        aTGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCa                                              <  1:5428/140‑1 (MQ=255)                                                                                                      cTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGAc                                >  1:385292/1‑140 (MQ=255)                                                                                                      cTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGAc                                >  2:383087/1‑140 (MQ=255)                                                                                                      cTGGTGTGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGAc                                >  1:850547/1‑140 (MQ=255)                                                                                                         gtgtGACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACAtt                             >  2:874878/1‑140 (MQ=255)                                                                                                           gtgACCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACAtttt                           >  2:231184/1‑140 (MQ=255)                                                                                                              aCCGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTc                        >  1:927409/1‑140 (MQ=255)                                                                                                               ccGATCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCt                       <  2:839362/140‑1 (MQ=255)                                                                                                                  aTCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTa                    <  2:126966/140‑1 (MQ=255)                                                                                                                  aTCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTa                    <  2:336171/140‑1 (MQ=255)                                                                                                                   tCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAc                   <  2:385292/140‑1 (MQ=255)                                                                                                                   tCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAc                   <  1:874878/140‑1 (MQ=255)                                                                                                                   tCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAc                   <  2:525885/140‑1 (MQ=255)                                                                                                                   tCACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAc                   <  2:107165/140‑1 (MQ=255)                                                                                                                    cACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAcc                  <  2:643878/140‑1 (MQ=255)                                                                                                                    cACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAcc                  <  2:850547/140‑1 (MQ=255)                                                                                                                    cACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAcc                  <  1:383087/140‑1 (MQ=255)                                                                                                                    cACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAcc                  <  1:754244/140‑1 (MQ=255)                                                                                                                    cACCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAcc                  <  1:231184/140‑1 (MQ=255)                                                                                                                     aCCATCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTTGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCa                 <  1:1062696/140‑1 (MQ=255)                                                                                                                         tCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACAtttttt                         >  1:1003865/1‑128 (MQ=255)                                                                                                                         tCAGAACAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACAtttttt                         <  2:1003865/128‑1 (MQ=255)                                                                                                                             aaCAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGa         >  2:234368/1‑140 (MQ=255)                                                                                                                             aaCAGCCAGCCCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGa         >  1:1297408/1‑140 (MQ=255)                                                                                                                                 gccagccCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGc     >  2:736627/1‑140 (MQ=255)                                                                                                                                 gccagccCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGc     >  1:173950/1‑140 (MQ=255)                                                                                                                                  ccagccCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGcc    >  2:693141/1‑140 (MQ=255)                                                                                                                                    agccCAAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTc  >  1:751716/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                       GCAGCGTCTGAATGACGAAATGTTTGAGGTCACTTTCTGGTGGCGTGATCCCCAAGGTTCTGAAGAATACTCGACGATAAAGCGCGTATGGGTCTACATCACTGGTGTGACCGATCACCATCAGAACAGCCAGCCCCAGTCGATGCAGCGAATTGCAGGCACTAACGTCTGGCAGTGGACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTC  >  NC_000913/612796‑613066

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