Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	219,876	T→G	K40N (AAA→AAC)	rcsF ←	putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	219,876	0	T	G	100.0%	213.2 / NA	60	K40N (AAA→AAC)	rcsF	putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/30);  total (30/30)

ACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCAAAAC  >  NC_000913/219729‑220022                                                                                                                                                    |                                                                                                                                                   aCGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGgttt                                                                                                                                                 <  1:711676/151‑1 (MQ=255)     gCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAg                                                                                                                                             <  1:877970/151‑1 (MQ=255)     gCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAg                                                                                                                                             <  2:281737/151‑1 (MQ=255)         ttGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTc                                                                                                                                         >  1:585040/1‑151 (MQ=255)            gAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGc                                                                                                                                      >  2:754999/1‑151 (MQ=255)             aaTGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCt                                                                                                                                     >  2:796931/1‑151 (MQ=255)               tGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCttt                                                                                                                                   >  1:406252/1‑151 (MQ=255)                 cTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCCGCGCTTTCCGGTTTGCAGGCTCCGCttttg                                                                                                                                 >  1:764737/1‑149 (MQ=255)                      cggAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTgg                                                                                                                            >  2:805986/1‑151 (MQ=255)                                tGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGGTTGCAGGCTCCGCTTTCGGCTGGGGGGCAGTg                                                                                                                   >  1:794499/1‑150 (MQ=255)                                    tAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCttt                                                                                                               <  2:736945/150‑1 (MQ=255)                                    tAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTg                                                                                                              >  2:59433/1‑151 (MQ=255)                                          ccTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACggg                                                                                                        >  2:146488/1‑151 (MQ=255)                                             ggCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTc                                                                                                     <  1:329349/151‑1 (MQ=255)                                                   agTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAgg                                                                                               >  2:379994/1‑151 (MQ=255)                                                      cGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGa                                                                                            <  2:740056/151‑1 (MQ=255)                                                      cGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGa                                                                                            <  1:904338/151‑1 (MQ=255)                                                       gCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGTGGAt                                                                                           <  1:505446/151‑1 (MQ=255)                                                              aCTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGcct                                                                                     >  1:251688/1‑148 (MQ=255)                                                              aCTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTa                                                                                    <  1:430788/151‑1 (MQ=255)                                                                     cGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATgg                                                                              >  2:738129/1‑150 (MQ=255)                                                                        gaTCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACa                                                                          <  2:483068/151‑1 (MQ=255)                                                                             cgGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGc                                                                     <  1:901039/151‑1 (MQ=255)                                                                                    gTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCAt                                                              >  1:599035/1‑151 (MQ=255)                                                                                       tGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATg                                                             >  1:207723/1‑149 (MQ=255)                                                                                        gCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAg                                                           <  1:18975/150‑1 (MQ=255)                                                                                            aCTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGct                                                       >  1:287332/1‑150 (MQ=255)                                                                                             cTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGcta                                                      <  1:630169/150‑1 (MQ=255)                                                                                              tAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGGAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTTACGTGGAACAGCCGCTTAGCATGAGTGctact                                                    >  1:703912/1‑151 (MQ=255)                                                                                                aTTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGGGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGctacta                                                   >  1:431532/1‑150 (MQ=255)                                                                                                aTTCTTCTGCATTGGTATAAATTCGGACCGGCGTCGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCCGTGCCTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGctccta                                                   >  1:212752/1‑150 (MQ=255)                                                                                                        gCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATc                                           <  2:289308/150‑1 (MQ=255)                                                                                                            tGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTa                                       <  2:909465/150‑1 (MQ=255)                                                                                                              gTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTaa                                      >  1:412852/1‑149 (MQ=255)                                                                                                                 taAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCa                                  >  2:569679/1‑150 (MQ=255)                                                                                                                  aaaTTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACg                                <  2:448115/151‑1 (MQ=255)                                                                                                                   aaTTCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGc                               <  2:73861/151‑1 (MQ=255)                                                                                                                     ttCGGACCGGCGTGGCGTGCGGCGCTTTCGGGTTTGCAGGGTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTTACATGGGACAGCCGCTTCGCATGAGTGCTACTAAACAGATCGGTAAAGCACGGAt                             >  1:839583/1‑151 (MQ=255)                                                                                                                     ttCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCa                              >  1:713694/1‑150 (MQ=255)                                                                                                                      tCGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCCTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCAt                             >  1:472128/1‑150 (MQ=255)                                                                                                                       cGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATa                            <  1:404093/150‑1 (MQ=255)                                                                                                                       cGGACCGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATa                            <  2:167873/150‑1 (MQ=255)                                                                                                                            cGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTAttt                      >  1:814223/1‑151 (MQ=255)                                                                                                                            cGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTAttt                      <  1:460963/151‑1 (MQ=255)                                                                                                                            cGGCGTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTAttt                      <  1:46900/151‑1 (MQ=255)                                                                                                                                gTGGCGCGCGGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTc                  <  1:609786/151‑1 (MQ=255)                                                                                                                                     gcgcgGCGCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGt             <  1:291488/151‑1 (MQ=255)                                                                                                                                       gcggcgCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTAt           >  2:22253/1‑151 (MQ=255)                                                                                                                                        cggcgCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATg          >  1:45905/1‑151 (MQ=255)                                                                                                                                        cggcgCTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATg          <  2:862904/151‑1 (MQ=255)                                                                                                                                         ggcgCTTTCGGGTTTGGAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAGCGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATg          <  2:799843/150‑1 (MQ=255)                                                                                                                                            gcTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGGCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGa         >  2:208426/1‑148 (MQ=255)                                                                                                                                            gcTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCa      <  2:534372/151‑1 (MQ=255)                                                                                                                                            gcTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCa      <  2:585058/151‑1 (MQ=255)                                                                                                                                            gcTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCa      <  2:197740/151‑1 (MQ=255)                                                                                                                                            gcTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGCGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATc       >  2:477742/1‑150 (MQ=255)                                                                                                                                             cTTTCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGAt        >  1:288826/1‑148 (MQ=255)                                                                                                                                               ttCGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCaaa    >  2:84100/1‑150 (MQ=255)                                                                                                                                                 cGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCAAAAc  <  1:661901/150‑1 (MQ=255)                                                                                                                                                 cGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCAAAAc  <  1:79593/150‑1 (MQ=255)                                                                                                                                                 cGGGTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCAAAAc  <  2:152417/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGTGCTTTGAACGGGTTCGACAGGGGATCTGCTTAACATGGAACAGCCGCTTAGCATGAGTGCTACTAAACAGATCGGTAAAGCACGCATAGTATTTCCTCAATGTATGATCAAAAC  >  NC_000913/219729‑220022

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