Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	199,912	G→A	G662D (GGC→GAC)	bamA →	BamABCDE complex OM biogenesis outer membrane pore‑forming assembly factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	199,912	0	G	A	100.0%	261.6 / NA	74	G662D (GGC→GAC)	bamA	BamABCDE complex OM biogenesis outer membrane pore‑forming assembly factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/38);  total (36/38)

CTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATCGGATGATGCTGT  >  NC_000913/199777‑200044                                                                                                                                        |                                                                                                                                     cTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTAAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACtt                                                                                                                                   >  1:951140/1‑139 (MQ=255)     tgtgCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAg                                                                                                                               <  1:1081956/139‑1 (MQ=255)            atcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATa                                                                                                                        >  1:779913/1‑139 (MQ=255)            atcgatGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATa                                                                                                                        >  1:21767/1‑139 (MQ=255)                   aCGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAg                                                                                                                               >  2:867057/1‑125 (MQ=255)                   aCGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAg                                                                                                                               <  1:867057/125‑1 (MQ=255)                           aaaTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAg                                                                                                         >  1:58054/1‑139 (MQ=255)                             aTGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCa                                                                                                       <  2:295872/139‑1 (MQ=255)                                ggTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGtt                                                                                                    >  2:336960/1‑139 (MQ=255)                                ggTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGtt                                                                                                    >  1:1104665/1‑139 (MQ=255)                                  ttgttCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTa                                                                                                  <  1:1038195/139‑1 (MQ=255)                                            gCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCAt                                                                                        <  2:21767/139‑1 (MQ=255)                                            gCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCAt                                                                                        <  1:347381/139‑1 (MQ=255)                                            gCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCAt                                                                                        <  1:350226/139‑1 (MQ=255)                                               tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAg                                                                                                                               >  1:138997/1‑97 (MQ=255)                                               tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAg                                                                                                                               <  2:138997/97‑1 (MQ=255)                                               tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAg                                                                                     >  2:311279/1‑139 (MQ=255)                                               tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAg                                                                                     >  2:637279/1‑139 (MQ=255)                                               tACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAg                                                                                     >  2:96394/1‑139 (MQ=255)                                                 cccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGc                                                                                   <  2:1104665/139‑1 (MQ=255)                                                 cccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGc                                                                                   <  1:774865/139‑1 (MQ=255)                                                 cccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGc                                                                                   <  1:914578/139‑1 (MQ=255)                                                  ccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGGCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGcc                                                                                  <  1:792015/139‑1 (MQ=255)                                                  ccGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGcc                                                                                  >  1:1085248/1‑139 (MQ=255)                                                    gCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAg                                                                                >  1:423184/1‑139 (MQ=255)                                                     cTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGt                                                                               <  1:637279/139‑1 (MQ=255)                                                             aTGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGa                                                                       <  2:603677/139‑1 (MQ=255)                                                              tGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGAt                                                                      <  2:59619/139‑1 (MQ=255)                                                               ggTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATc                                                                     <  2:1085248/139‑1 (MQ=255)                                                                 tGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCg                                                                   <  2:779913/139‑1 (MQ=255)                                                                  gATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATc                                                                     >  1:491001/1‑136 (MQ=255)                                                                  gATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATc                                                                     <  2:491001/136‑1 (MQ=255)                                                                     ggTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACt                                                               <  1:493290/139‑1 (MQ=255)                                                                       tttAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTAt                                                             >  2:396794/1‑139 (MQ=255)                                                                         tAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGa                                                           >  1:42261/1‑139 (MQ=255)                                                                          aGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGAt                                                          >  1:70445/1‑139 (MQ=255)                                                                             cggcAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTAc                                                       <  2:50877/139‑1 (MQ=255)                                                                               gcAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGa                                                     >  1:608693/1‑139 (MQ=255)                                                                                cAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGc                                                                                   >  1:833341/1‑108 (MQ=255)                                                                                cAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGc                                                                                   <  2:833341/108‑1 (MQ=255)                                                                                 aaaGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTAt                                                             <  1:91339/129‑1 (MQ=255)                                                                                 aaaGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTAt                                                             >  2:91339/1‑129 (MQ=255)                                                                                    gagaTGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAAtgtg                                                <  2:517871/139‑1 (MQ=255)                                                                                    gagaTGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAAtgtg                                                <  2:245379/139‑1 (MQ=255)                                                                                    gagaTGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAAtgtg                                                <  2:461828/139‑1 (MQ=255)                                                                                                tACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACg                                    <  2:58054/139‑1 (MQ=255)                                                                                                  cGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGgc                                  >  1:189709/1‑139 (MQ=255)                                                                                                    agaACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGa                                                                 >  1:903322/1‑106 (MQ=255)                                                                                                    agaACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGa                                                                 <  2:903322/106‑1 (MQ=255)                                                                                                      aaCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCaa                                                                                                                          <  1:101485/47‑1 (MQ=255)                                                                                                      aaCTTCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCaa                                                                                                                          >  2:101485/1‑47 (MQ=255)                                                                                                          tCTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGa                          >  2:234689/1‑139 (MQ=255)                                                                                                           cTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGAc                         <  1:2337/139‑1 (MQ=255)                                                                                                           cTATGCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGAc                         <  2:423184/139‑1 (MQ=255)                                                                                                               gCCGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCtgt                     >  1:261541/1‑139 (MQ=255)                                                                                                                ccGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTAt                                                             >  2:71457/1‑98 (MQ=255)                                                                                                                ccGGTGGTTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTAt                                                             <  1:71457/98‑1 (MQ=255)                                                                                                                      gtTCCAGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATc              <  2:172445/139‑1 (MQ=255)                                                                                                                           aGCACCGTGCGTGACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATCGgatg         <  2:42261/139‑1 (MQ=255)                                                                                                                                cgtgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATCGGATGATGCt    >  2:649518/1‑139 (MQ=255)                                                                                                                                 gtgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATCGGATGATGCTg   >  1:990959/1‑139 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAAtt                                                                           >  1:190005/1‑66 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAAtt                                                                           <  2:190005/66‑1 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAAtgtg                                                >  1:808004/1‑93 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAAtgtg                                                <  2:808004/93‑1 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACgg                                   <  1:261420/106‑1 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACgg                                   >  2:261420/1‑106 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGAc                         >  1:540353/1‑116 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGAc                         <  2:540353/116‑1 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATCGGATGATGCTGt  >  1:472618/1‑139 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATCGGATGATGCTGt  >  2:877946/1‑139 (MQ=255)                                                                                                                                  tgcgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATCGGATGATGCTGt  >  1:137794/1‑139 (MQ=255)                                                                                                                                    cgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGc                                                                                   <  2:761042/56‑1 (MQ=255)                                                                                                                                    cgtgACTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGc                                                                                   >  1:761042/1‑56 (MQ=255)                                                                                                                                        |                                                                                                                                     CTTATGTGCCGATCGATGACGATCACAAATGGGTTGTTCTGGGGCGTACCCGCTGGGGTTATGGTGATGGTTTAGGCGGCAAAGAGATGCCGTTCTACGAGAACTTCTATGCCGGTGGTTCCAGCACCGTGCGTGGCTTCCAGTCCAATACCATTGGTCCGAAAGCAGTTTACTTCCCGCATCAGGCCAGTAATTATGATCCGGACTATGATTACGAATGTGCGACTCAGGACGGCGCGAAAGACCTGTGTAAATCGGATGATGCTGT  >  NC_000913/199777‑200044

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