BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	817,654	+G	coding (377/513 nt)	moaB →	molybdopterin biosynthesis protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	817,651	1	.	G	100.0%	330.6 / NA	101	A125G (GCG→GGG)	moaB	molybdopterin biosynthesis protein B
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (60/41); total (60/41)

Only 100 of 128 total aligned reads displayed.

GTTCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑G‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCA  >  NC_000913/817559‑817738

gttCGACCGTGAAGTTTAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGc                                                                                           <  2:2344159/93‑1 (MQ=255)

gttCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGc                                                                                           >  2:3855432/1‑93 (MQ=255)

gttCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGc                                                                                           <  1:821937/93‑1 (MQ=255)

 ttCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑g                                                                                         >  2:823794/1‑93 (MQ=255)

  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTGCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑g‑g                                                                                       <  1:3301715/93‑1 (MQ=255)

  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑g‑g                                                                                       <  1:3372148/93‑1 (MQ=255)

  tCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑g‑g                                                                                       <  1:1763895/93‑1 (MQ=255)

     aCCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑g‑gggc                                                                                    >  2:4172562/1‑91 (MQ=255)

      ccGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑g‑gggcg                                                                                   >  2:2255498/1‑90 (MQ=255)

      ccGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑g‑gggcg                                                                                   >  1:796936/1‑90 (MQ=255)

       cGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑g‑gggcgg                                                                                  >  1:223217/1‑89 (MQ=255)

         tGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑GCGGCGTCg                                                                                 >  1:389783/1‑93 (MQ=255)

           aaGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGcc                                                                               <  2:225090/93‑1 (MQ=255)

           aaGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGcc                                                                               >  1:3142808/1‑93 (MQ=255)

           aaGTTGAAGGTTTTGGGGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGcc                                                                               >  2:2019965/1‑93 (MQ=255)

             gTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCcaa                                                                             >  2:3880096/1‑93 (MQ=255)

             gTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCcaa                                                                             <  1:3340961/93‑1 (MQ=255)

             gTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCcaa                                                                             >  2:750552/1‑93 (MQ=255)

               tGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCcgaca                                                                           >  2:942341/1‑93 (MQ=255)

                gAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCcaacaa                                                                          >  1:1266279/1‑93 (MQ=255)

                  aGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACaaaa                                                                        <  2:4386211/93‑1 (MQ=255)

                  aGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACaaaa                                                                        >  2:2830522/1‑93 (MQ=255)

                     ttttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGc                                                                     <  1:2662158/93‑1 (MQ=255)

                     ttttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGc                                                                     <  2:492496/93‑1 (MQ=255)

                     ttttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGc                                                                     <  1:364505/93‑1 (MQ=255)

                     ttttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGc                                                                     >  2:712310/1‑93 (MQ=255)

                      tttGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCt                                                                    <  1:718693/93‑1 (MQ=255)

                            gAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGAtttt                                                              >  2:2124710/1‑93 (MQ=255)

                            gAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGAtttt                                                              <  1:2836684/93‑1 (MQ=255)

                             aaGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTc                                                             >  1:4577200/1‑93 (MQ=255)

                             aaGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTc                                                             >  1:3881240/1‑93 (MQ=255)

                              aGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCg                                                            <  1:248979/93‑1 (MQ=255)

                              aGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCg                                                            >  2:2639840/1‑93 (MQ=255)

                                tgtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGcc                                                          <  1:2432829/93‑1 (MQ=255)

                                 tttCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCa                                                         <  2:503579/92‑1 (MQ=255)

                                 gtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCa                                                         <  2:674094/93‑1 (MQ=255)

                                 gtTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCa                                                         >  2:1311844/1‑93 (MQ=255)

                                    tcGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGc                                                      <  2:2644286/92‑1 (MQ=255)

                                    ccGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGc                                                      <  1:3331646/93‑1 (MQ=255)

                                     cGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGcc                                                     >  1:2292187/1‑93 (MQ=255)

                                     cGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGcc                                                     >  2:2327497/1‑93 (MQ=255)

                                     cGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGcc                                                     >  2:2327498/1‑93 (MQ=255)

                                     cGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGcc                                                     >  2:449448/1‑93 (MQ=255)

                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    >  1:3228886/1‑93 (MQ=255)

                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    >  2:287060/1‑93 (MQ=255)

                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    <  2:3509171/93‑1 (MQ=255)

                                      ggATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCg                                                    <  2:3509170/93‑1 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  1:1912754/1‑93 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  1:4140989/1‑93 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  2:3979759/1‑93 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  1:1355443/1‑93 (MQ=255)

                                       gATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCgg                                                   >  2:1413886/1‑93 (MQ=255)

                                        aTGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCACCAAAACGCTGTTTTTGGCCATGCCggt                                                  >  2:3048999/1‑92 (MQ=255)

                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 >  1:953774/1‑93 (MQ=255)

                                         tgttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGt                                                 >  1:328244/1‑93 (MQ=255)

                                          gttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCTGGGtt                                                >  1:3173455/1‑93 (MQ=255)

                                          gttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGtt                                                >  1:2279353/1‑93 (MQ=255)

                                           ttgtCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTc                                               <  1:2840437/93‑1 (MQ=255)

                                                gTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCa                                          >  2:2799458/1‑93 (MQ=255)

                                                 tttGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGAccca                                         >  1:1773565/1‑91 (MQ=255)

                                                  ttGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCaaa                                        <  2:103952/93‑1 (MQ=255)

                                                     aaGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCg                                     <  2:3352872/93‑1 (MQ=255)

                                                     aaGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCg                                     <  1:3249154/93‑1 (MQ=255)

                                                      agagATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGt                                    >  1:2582870/1‑93 (MQ=255)

                                                       gagaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTg                                   >  1:753875/1‑93 (MQ=255)

                                                       gagaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTg                                   >  2:3156025/1‑93 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 <  1:213362/93‑1 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 >  1:1455409/1‑93 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 <  2:1730078/93‑1 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 >  2:2685149/1‑93 (MQ=255)

                                                         gaTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGcc                                 >  1:1452531/1‑93 (MQ=255)

                                                          aTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCg                                >  2:2612369/1‑93 (MQ=255)

                                                          aTTGGCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGGGCCg                                >  2:1874287/1‑93 (MQ=255)

                                                             ggCACTTCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTAc                             >  2:151111/1‑93 (MQ=255)

                                                                  ttCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCAt                        >  2:580313/1‑93 (MQ=255)

                                                                  ttCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCAt                        <  1:1156804/93‑1 (MQ=255)

                                                                   tCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATg                       <  1:1805655/93‑1 (MQ=255)

                                                                   tCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATg                       >  2:1593383/1‑93 (MQ=255)

                                                                   tCCACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATg                       <  1:3738913/93‑1 (MQ=255)

                                                                    tcACGTTGCAATCCCGTGCGGGAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATgg                      <  2:1710088/92‑1 (MQ=255)

                                                                    ccACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATgg                      >  2:3300953/1‑93 (MQ=255)

                                                                     cACGTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATggg                     >  2:1577773/1‑93 (MQ=255)

                                                                        gTTGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGACGGGTTCGACCAAAGCGTGCCGTACCGCATGGGgaa                  >  2:195595/1‑93 (MQ=255)

                                                                         ttGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaaa                 >  2:90878/1‑93 (MQ=255)

                                                                         ttGCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGaaaa                 >  1:4050019/1‑93 (MQ=255)

                                                                           gCAATCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAata               >  1:3133718/1‑93 (MQ=255)

                                                                             aaTCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatc             >  2:4041581/1‑93 (MQ=255)

                                                                              aTCTCGTGCGGTAGCGG‑GGCGTCGCCACCAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAatataa            >  1:1277027/1‑91 (MQ=255)

                                                                              aTCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatca            >  1:1277026/1‑93 (MQ=255)

                                                                              aTCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatca            >  2:4325573/1‑93 (MQ=255)

                                                                              aTCTCGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatca            >  1:910509/1‑93 (MQ=255)

                                                                               tctcGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           <  1:28086/93‑1 (MQ=255)

                                                                               tctcGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcat           >  2:4373946/1‑93 (MQ=255)

                                                                                ctcGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATatcatc          <  1:199189/93‑1 (MQ=255)

                                                                                  cGTGCGGTAGCGG‑GGCGTCGCCAACAACACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcgc        <  1:3381730/93‑1 (MQ=255)

                                                                                  cGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcgc        <  1:1876095/93‑1 (MQ=255)

                                                                                  cGTGCGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATcgc        <  1:1456924/93‑1 (MQ=255)

                                                                                      cGGTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgccg    <  2:3991390/93‑1 (MQ=255)

                                                                                       ggTAGCGG‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGcgccgc   >  2:1379102/1‑93 (MQ=255)

                                                                                        gTAGCGG‑GGCCTCGCCAACTAAACGCTGATTTTCGCCATGCCGGGATCGACAAAACAGTGCCGAACCGCTTGGTAAAATATCATAGCGCCGCa  >  2:2250875/1‑93 (MQ=255)

GTTCGACCGTGAAGTTGAAGGTTTTGGTGAAGTGTTCCGGATGTTGTCGTTTGAAGAGATTGGCACTTCCACGTTGCAATCTCGTGCGGTAGC‑G‑GGCGTCGCCAACAAAACGCTGATTTTCGCCATGCCGGGTTCGACCAAAGCGTGCCGTACCGCATGGGAAAATATCATCGCGCCGCA  >  NC_000913/817559‑817738

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