BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	966,368	A→G	A683A (GCA→GCG)	ycaI →	ComEC family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	966,368	0	A	G	100.0%	440.3 / NA	125	A683A (GCA→GCG)	ycaI	ComEC family inner membrane protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (67/58); total (67/58)

Only 100 of 125 total aligned reads displayed.

GCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCC  >  NC_000913/966272‑966463

gcagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTc                                                                                              >  2:2739470/1‑100 (MQ=255)

 cagcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCg                                                                                             <  2:152825/100‑1 (MQ=255)

  agcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgc                                                                                            >  1:2030370/1‑100 (MQ=255)

  agcaAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgc                                                                                            >  2:1492231/1‑100 (MQ=255)

   gcaATTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCAGCGGTTCCATTAAGTCAGCGAGTGAATGGAAAAGTGGCGCTCgca                                                                                           <  1:2429532/100‑1 (MQ=255)

      aGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcg                                                                                        >  2:3109833/1‑100 (MQ=255)

       gTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcgg                                                                                       <  1:2630902/100‑1 (MQ=255)

         gCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggca                                                                                     <  1:361747/100‑1 (MQ=255)

           aGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcatc                                                                                   <  1:2832496/100‑1 (MQ=255)

           aGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcatc                                                                                   <  1:87495/100‑1 (MQ=255)

            ggCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCgcatcggcatcg                                                                                  >  2:2624597/1‑100 (MQ=255)

                aCATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCt                                                                              >  2:1727898/1‑100 (MQ=255)

                 cATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCta                                                                             >  2:2416319/1‑100 (MQ=255)

                  aTTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCtat                                                                            <  2:2131409/100‑1 (MQ=255)

                     gCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAAc                                                                         >  1:2152450/1‑100 (MQ=255)

                      cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                        >  2:1227332/1‑100 (MQ=255)

                      cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                        >  1:96171/1‑100 (MQ=255)

                      cTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACg                                                                        >  1:2456655/1‑100 (MQ=255)

                        tCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCa                                                                      >  2:1767530/1‑100 (MQ=255)

                        tCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCa                                                                      >  1:747/1‑100 (MQ=255)

                        tCAGGTACCTCACCATGGCAGTAACACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCa                                                                      <  2:1347855/100‑1 (MQ=255)

                         cAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCAt                                                                     <  1:155584/100‑1 (MQ=255)

                          aGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTCATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACCCATg                                                                    >  1:1576853/1‑100 (MQ=255)

                          aGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATg                                                                    <  2:3360292/100‑1 (MQ=255)

                           ggTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATgg                                                                   >  1:825621/1‑100 (MQ=255)

                             tACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCg                                                                 <  2:906818/100‑1 (MQ=255)

                               ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                               >  2:2137915/1‑100 (MQ=255)

                               ccTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGAc                                                               >  2:3304076/1‑100 (MQ=255)

                                cTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACt                                                              >  1:2516198/1‑100 (MQ=255)

                                 tCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTg                                                             >  1:1010352/1‑100 (MQ=255)

                                  cACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGc                                                            >  1:3009889/1‑100 (MQ=255)

                                    ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                          >  1:2641624/1‑100 (MQ=255)

                                    ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                          >  1:2667582/1‑100 (MQ=255)

                                    ccATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGccc                                                          <  1:2033081/100‑1 (MQ=255)

                                       tGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCctct                                                       >  1:1830377/1‑100 (MQ=255)

                                          cAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAAc                                                    >  2:3396264/1‑100 (MQ=255)

                                            gTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACaa                                                  >  2:3149168/1‑100 (MQ=255)

                                              aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                <  1:1518730/100‑1 (MQ=255)

                                              aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                <  1:1627495/100‑1 (MQ=255)

                                              aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                <  2:2169697/100‑1 (MQ=255)

                                              aaTACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAg                                                <  1:2657074/100‑1 (MQ=255)

                                                tACCTCATCATCGTTGCCATTAATTCGGCGAGTGAATGGAAAAGTGGCGCTCGCAGCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTGACAAAGt                                               <  1:3284964/99‑1 (MQ=255)

                                                 aCCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTa                                             <  1:3292330/100‑1 (MQ=255)

                                                   cTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                           <  2:1510010/100‑1 (MQ=255)

                                                   cTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                           <  2:441326/100‑1 (MQ=255)

                                                   cTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGAGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAg                                           >  1:1903323/1‑100 (MQ=255)

                                                    tcatccTCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGtt                                              <  1:2905376/96‑1 (MQ=255)

                                                     catcatCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCa                                         >  1:132136/1‑100 (MQ=255)

                                                     catcatCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCa                                         >  2:3101382/1‑100 (MQ=255)

                                                         atcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAAGGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                                     >  2:3305597/1‑100 (MQ=255)

                                                         atcGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGc                                     <  1:1712397/100‑1 (MQ=255)

                                                           cGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTa                                   <  2:3146989/100‑1 (MQ=255)

                                                           cGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTa                                   <  1:1727875/100‑1 (MQ=255)

                                                             ttGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATc                                 <  1:416989/100‑1 (MQ=255)

                                                              tGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCa                                <  1:1509251/100‑1 (MQ=255)

                                                                 cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                             >  1:1160342/1‑100 (MQ=255)

                                                                 cATTAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACt                             >  1:421922/1‑100 (MQ=255)

                                                                    tAATTCAGCGAGTTAATGGAAAAGTGGCGCTCGCATCGACATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCa                          >  1:2147271/1‑100 (MQ=255)

                                                                    tAATTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCa                          <  1:3295527/100‑1 (MQ=255)

                                                                     aaTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                         >  2:2554107/1‑100 (MQ=255)

                                                                     aaTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                         >  1:2592896/1‑100 (MQ=255)

                                                                     aaTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                         >  1:845606/1‑100 (MQ=255)

                                                                     aaTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCaa                         >  2:1186530/1‑100 (MQ=255)

                                                                      aTTCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAg                        >  2:688757/1‑100 (MQ=255)

                                                                       ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                       <  2:3196799/100‑1 (MQ=255)

                                                                       ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                       <  1:2137884/100‑1 (MQ=255)

                                                                       ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                       <  2:2978274/100‑1 (MQ=255)

                                                                       ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAgg                       <  2:2168322/100‑1 (MQ=255)

                                                                       ttCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGAATCGCGCCCTAACGCATGGCGACTGCCGTCTACCAAAGGTAAGCATCGCTATCAACTGCAAgg                       >  1:1258260/1‑100 (MQ=255)

                                                                        tCAGCGAGTGAATGGAAAAGTGGCGCTCGCCTCCGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCGACTGCAAGGa                      >  1:2281538/1‑100 (MQ=255)

                                                                        tCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGa                      >  1:3406066/1‑100 (MQ=255)

                                                                        tCAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGa                      <  1:779264/100‑1 (MQ=255)

                                                                         cAGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGat                     >  1:1380755/1‑100 (MQ=255)

                                                                          aGCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGata                    >  2:675760/1‑100 (MQ=255)

                                                                           gCGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGatat                   <  2:1136694/100‑1 (MQ=255)

                                                                            cGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATc                  <  2:708070/100‑1 (MQ=255)

                                                                            cGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATc                  <  2:751440/100‑1 (MQ=255)

                                                                            cGAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCCCTAACAAAGCTAAGCATCGCTATCAACTGCAAGGATATc                  <  2:3173371/100‑1 (MQ=255)

                                                                             gAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCa                 >  2:2437788/1‑100 (MQ=255)

                                                                             gAGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCa                 >  2:2713464/1‑100 (MQ=255)

                                                                              aGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCaa                >  1:2788711/1‑100 (MQ=255)

                                                                              aGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCaa                >  1:710508/1‑100 (MQ=255)

                                                                              aGTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCaa                >  1:1973727/1‑100 (MQ=255)

                                                                               gTGAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAAt               >  2:1714395/1‑100 (MQ=255)

                                                                                tGAATGGAAAAGTGGCGCTCGCATCGGCATCGCTCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATg              <  1:1059333/100‑1 (MQ=255)

                                                                                 gAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg             >  2:651701/1‑100 (MQ=255)

                                                                                 gAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATgg             >  2:759058/1‑100 (MQ=255)

                                                                                 gAATGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGACGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATg              <  1:736904/99‑1 (MQ=255)

                                                                                  aaTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGa            <  1:3305546/100‑1 (MQ=255)

                                                                                  aaTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGa            <  2:250193/100‑1 (MQ=255)

                                                                                  aaTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGa            >  2:1467045/1‑100 (MQ=255)

                                                                                   aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt           >  2:456210/1‑100 (MQ=255)

                                                                                   aTGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAt           <  2:488170/100‑1 (MQ=255)

                                                                                    tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt          <  2:2909398/100‑1 (MQ=255)

                                                                                    tGGAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGAtt          <  1:1767505/100‑1 (MQ=255)

                                                                                      gAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa        <  2:3009938/100‑1 (MQ=255)

                                                                                      gAAAAGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGa        <  2:626967/100‑1 (MQ=255)

                                                                                          aGTGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACt    <  2:474004/100‑1 (MQ=255)

                                                                                            tGGCGCTCGCATCGGCATCGCGCTCTAACGCATGGCGACTGCCCTATACCAAAGTTCAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTcc  <  2:1897132/100‑1 (MQ=255)

                                                                                            tGGCGCTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTcc  <  1:1624966/100‑1 (MQ=255)

GCAGCAAGTGCAGGCAACATTGCTTCAGGTACCTCACCATGGCAGTAATACCTCATCATCGTTGCCATTAATTCAGCGAGTGAATGGAAAAGTGGCACTCGCATCGGCATCGCGCTATAACGCATGGCGACTGCCCTCTAACAAAGTTAAGCATCGCTATCAACTGCAAGGATATCAATGGATTGATACTCC  >  NC_000913/966272‑966463

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